Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,169,338 G→C 28.7% G84G (GGG→GGC tufA → protein chain elongation factor EF‑Tu

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077794,169,3380G→C28.7% 165.6 ‑0.3 ‑1.4167G84G (GGG→GGCtufAprotein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand):  new base (20/28):  ref base (54/65):  total (74/93)
Fisher's exact test strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.37e-02

Only 100 of 179 total aligned reads displayed.

ACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTG  >  NC_007779/4169277‑4169438
                                                              |                                                                                                    
acacTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTgctgct                                                               <  1:8969986/101‑1 (MQ=11)
acacTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTgctgct                                                               >  1:1183981/1‑101 (MQ=11)
acacTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTgctgct                                                               <  1:5916857/101‑1 (MQ=11)
 cgtctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTc                                                              <  1:10044258/94‑1 (MQ=2)
  acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                                                             >  1:210590/1‑101 (MQ=32)
  acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                                                             >  1:10121127/1‑101 (MQ=32)
  acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                                                             >  1:5219838/1‑101 (MQ=32)
  acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                                                             <  1:5905181/101‑1 (MQ=32)
   tctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAg                                                            <  1:7199432/96‑1 (MQ=2)
    ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                                                           <  1:14000799/101‑1 (MQ=32)
    ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                                                           <  1:11780201/101‑1 (MQ=32)
    ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                                                           <  1:9747902/101‑1 (MQ=32)
     tctcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGAt                                                          <  1:6112270/101‑1 (MQ=32)
      ctcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATg                                                         <  1:11286160/101‑1 (MQ=32)
      ctcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATg                                                         <  1:7090720/101‑1 (MQ=32)
      cgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATg                                                         <  1:6368232/99‑1 (MQ=2)
      cgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATg                                                         >  1:9878484/3‑101 (MQ=2)
       tcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATgg                                                        <  1:4860957/101‑1 (MQ=32)
         aCGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGAc                                                      <  1:6647398/101‑1 (MQ=2)
         aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc                                                      >  1:8146370/1‑101 (MQ=32)
         aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc                                                      >  1:4011159/1‑101 (MQ=32)
         aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc                                                      <  1:7557719/101‑1 (MQ=32)
          cGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACg                                                     <  1:4351784/101‑1 (MQ=2)
          cGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACg                                                     <  1:9330761/101‑1 (MQ=2)
          cGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACg                                                     >  1:6833539/1‑101 (MQ=32)
          cGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACg                                                     <  1:1604695/101‑1 (MQ=32)
          cGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACg                                                     <  1:492114/101‑1 (MQ=32)
          cGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACGGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACg                                                     <  1:11979417/101‑1 (MQ=17)
            ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc                                                   <  1:419487/101‑1 (MQ=32)
            ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc                                                   >  1:12726102/1‑101 (MQ=32)
            ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc                                                   <  1:14583839/101‑1 (MQ=32)
             tGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcg                                                  <  1:2354616/101‑1 (MQ=33)
             tGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcg                                                  >  1:5812396/1‑101 (MQ=32)
             tGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcg                                                  >  1:9836238/1‑101 (MQ=32)
              gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc                                                 <  1:13698253/101‑1 (MQ=33)
               aaTACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGcgcg                                                >  1:5617732/1‑101 (MQ=32)
                 tACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGAt                                              <  1:2873872/101‑1 (MQ=33)
                   cGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATcc                                            >  1:9102884/1‑101 (MQ=32)
                    gACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt                                           <  1:6553429/101‑1 (MQ=33)
                    gACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt                                           >  1:3948897/1‑101 (MQ=32)
                     acacCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTg                                          <  1:3853430/101‑1 (MQ=32)
                     acacCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTg                                          >  1:6232983/1‑101 (MQ=32)
                      cacCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTgg                                         >  1:6931216/1‑101 (MQ=32)
                       acCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGt                                        >  1:817576/1‑101 (MQ=32)
                        ccccGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTa                                       >  1:9563240/1‑101 (MQ=32)
                         cccGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAg                                      <  1:14750726/101‑1 (MQ=33)
                          ccGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGt                                     >  1:86558/1‑101 (MQ=32)
                           cGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGtt                                    <  1:13816608/101‑1 (MQ=33)
                           cGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGtt                                    >  1:14884607/1‑101 (MQ=32)
                            gACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTtg                                   <  1:6719997/101‑1 (MQ=32)
                             aagcGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTcgc                                  >  1:13124605/4‑101 (MQ=2)
                             aCCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTtgc                                  >  1:7872361/1‑101 (MQ=32)
                              cccGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTtgct                                 >  1:10464696/1‑101 (MQ=32)
                              agcGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTcgcc                                 <  1:6676893/99‑2 (MQ=2)
                              agcGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTcgcc                                 <  1:3864909/99‑2 (MQ=2)
                                cGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTcgccgc                               >  1:8610903/1‑101 (MQ=2)
                                 gTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTcgccgct                              <  1:4372936/101‑2 (MQ=2)
                                   cACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGAc                            <  1:2209297/101‑1 (MQ=33)
                                   cACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGAc                            <  1:523957/101‑1 (MQ=33)
                                   cACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGAc                            <  1:1482704/101‑1 (MQ=33)
                                    aCTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTAcc                           >  1:6280574/1‑100 (MQ=2)
                                     cTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTg                          >  1:11808937/1‑101 (MQ=32)
                                       aCGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGAc                        <  1:2796243/101‑1 (MQ=33)
                                         gCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACgg                      >  1:872556/1‑101 (MQ=32)
                                         gCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACgg                      <  1:5306182/101‑1 (MQ=33)
                                          cacaCGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGc                     <  1:14974403/101‑1 (MQ=33)
                                             acGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCg                  <  1:10832286/101‑1 (MQ=32)
                                             acGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCg                  <  1:879611/101‑1 (MQ=32)
                                               gTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGAt                <  1:10773161/101‑1 (MQ=32)
                                               gTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGAt                <  1:13410462/101‑1 (MQ=32)
                                               cgtcgacgctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGa                 >  1:7103621/11‑101 (MQ=2)
                                                gtcgacgctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGAt                <  1:8431336/92‑1 (MQ=2)
                                                 tcgacgctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATg               <  1:11353713/93‑1 (MQ=2)
                                                 tcgacgctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATg               <  1:12163158/93‑1 (MQ=2)
                                                 aGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGc              <  1:6370311/101‑1 (MQ=32)
                                                 aGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGc              >  1:2391647/1‑101 (MQ=32)
                                                 aGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGc              >  1:14126358/1‑101 (MQ=32)
                                                  gACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGcc             <  1:8007610/101‑1 (MQ=32)
                                                  gACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGcc             <  1:4576142/101‑1 (MQ=32)
                                                   gacgctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGcc             >  1:11302910/7‑101 (MQ=2)
                                                    cTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGc           >  1:4675831/1‑101 (MQ=32)
                                                     tGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCa          <  1:573813/101‑1 (MQ=32)
                                                      gctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCa          <  1:7384246/98‑1 (MQ=2)
                                                      gctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCa          <  1:6620665/98‑1 (MQ=2)
                                                      gCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAg         >  1:7694205/1‑101 (MQ=32)
                                                      gCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAg         >  1:6478025/1‑101 (MQ=32)
                                                       ctccGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAg         >  1:1814097/3‑101 (MQ=2)
                                                       cccGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGa        <  1:11657359/101‑1 (MQ=32)
                                                       cccGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGa        >  1:3295817/1‑101 (MQ=32)
                                                         ccGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGAc       >  1:3506145/1‑101 (MQ=32)
                                                          cGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACt      <  1:12523390/101‑1 (MQ=32)
                                                           ggggCACGCCGACTATTTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTc     <  1:10080753/101‑1 (MQ=17)
                                                           ggggCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTc     >  1:76001/1‑101 (MQ=33)
                                                           ggggCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTc     <  1:10025256/101‑1 (MQ=32)
                                                           gggcCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGAccc     <  1:5103559/101‑3 (MQ=2)
                                                           gggcCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATTCTGGTTGTCGCCGCTACCGACGGTCCGATGCCGCAGAccc     >  1:10098299/1‑99 (MQ=2)
                                                            gggCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCg    <  1:4691324/101‑1 (MQ=32)
                                                            gggCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCg    >  1:10196757/1‑101 (MQ=33)
                                                              gCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGGCTGACGGCCCGATGCCGCAGACTCGTg  <  1:9562760/101‑1 (MQ=17)
                                                              gCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTg  <  1:11733752/101‑1 (MQ=32)
                                                              |                                                                                                    
ACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTG  >  NC_007779/4169277‑4169438

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.