Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,169,297 G→C 13.0% D71H (GAC→CAC)  tufA → protein chain elongation factor EF‑Tu

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077794,169,2970G→C13.0% 38.7 ‑4.0 ‑0.5108D71H (GAC→CAC) tufAprotein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand):  new base (2/12):  ref base (47/47):  total (49/59)
Fisher's exact test strand distribution p-value = 1.90e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.27e-01

Only 100 of 115 total aligned reads displayed.

TAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCT  >  NC_007779/4169197‑4169397
                                                                                                    |                                                                                                    
tAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACg                                                                                                      >  1:2599932/1‑101 (MQ=32)
 aaaaCCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGa                                                                                                     <  1:3216004/101‑1 (MQ=32)
  aaaCCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGac                                                                                                    >  1:1839761/1‑101 (MQ=33)
  aaaCCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGac                                                                                                    <  1:3160919/101‑1 (MQ=32)
   aaCCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGaca                                                                                                   >  1:9371398/1‑101 (MQ=32)
   aaCCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGaca                                                                                                   >  1:5427766/1‑101 (MQ=32)
   aaCCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGaca                                                                                                   <  1:12028523/101‑1 (MQ=32)
   aaCCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGaca                                                                                                   <  1:3421420/101‑1 (MQ=32)
     ccTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAcc                                                                                                 <  1:829161/101‑1 (MQ=32)
      cTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAccc                                                                                                >  1:3256524/1‑101 (MQ=32)
      cTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAccc                                                                                                >  1:13547029/1‑101 (MQ=32)
       tACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAcccc                                                                                               <  1:6566669/101‑1 (MQ=32)
       tACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAcccc                                                                                               <  1:13746992/101‑1 (MQ=32)
        aCGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCg                                                                                              >  1:1889069/1‑101 (MQ=32)
        aCGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCg                                                                                              >  1:7478188/1‑101 (MQ=32)
         cggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGa                                                                                             <  1:11039692/101‑1 (MQ=32)
         cggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGa                                                                                             >  1:8834744/1‑101 (MQ=32)
          ggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAc                                                                                            >  1:1243737/1‑101 (MQ=32)
          ggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAc                                                                                            >  1:8640225/1‑101 (MQ=32)
          ggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAc                                                                                            <  1:5276358/101‑1 (MQ=32)
          ggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAc                                                                                            >  1:487070/1‑101 (MQ=33)
           gcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAcc                                                                                           <  1:2737201/101‑1 (MQ=32)
           gcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAcc                                                                                           >  1:1403550/1‑101 (MQ=32)
             ggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCg                                                                                         >  1:4912007/1‑101 (MQ=32)
             ggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCg                                                                                         <  1:5843119/101‑1 (MQ=32)
              gTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGt                                                                                        >  1:3945278/1‑101 (MQ=32)
               tgctgcTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTc                                                                                       <  1:13181247/101‑1 (MQ=32)
               tgctgcTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTc                                                                                       <  1:4915180/101‑1 (MQ=32)
                gctgctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAGAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCa                                                                                      <  1:11455975/101‑1 (MQ=17)
                 ctgctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTAAc                                                                                     <  1:4101343/101‑1 (MQ=17)
                 ctgctCGTGAATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCAc                                                                                     >  1:10130267/1‑101 (MQ=17)
                   gctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTa                                                                                   >  1:3697553/1‑101 (MQ=32)
                   gctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTa                                                                                   <  1:4904307/101‑1 (MQ=32)
                   gctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTa                                                                                   >  1:1225245/1‑101 (MQ=32)
                   gctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGCTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTa                                                                                   <  1:11858125/101‑1 (MQ=17)
                    ctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTAc                                                                                  <  1:150981/101‑1 (MQ=32)
                    ctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTAc                                                                                  >  1:8861442/1‑101 (MQ=32)
                    ctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTAc                                                                                  <  1:10242226/101‑1 (MQ=32)
                    ctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTAc                                                                                  <  1:10234179/101‑1 (MQ=32)
                     tCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACg                                                                                 <  1:1361812/101‑1 (MQ=32)
                      cGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGc                                                                                >  1:11291474/1‑101 (MQ=32)
                       gTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGca                                                                               >  1:14010239/1‑101 (MQ=31)
                       gTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGca                                                                               >  1:5277374/1‑101 (MQ=31)
                       gTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGca                                                                               >  1:15194421/1‑101 (MQ=31)
                        tGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGcac                                                                              >  1:8683457/1‑101 (MQ=11)
                        tGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGcac                                                                              <  1:5346277/101‑1 (MQ=11)
                                                                                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)
                                                                                 cgtctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTc                     <  1:3284404/94‑1 (MQ=2)
                                                                                  gtctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCa                    <  1:5636874/95‑1 (MQ=2)
                                                                                  acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                    >  1:10121127/1‑101 (MQ=32)
                                                                                  acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                    <  1:5905181/101‑1 (MQ=32)
                                                                                  acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                    >  1:5219838/1‑101 (MQ=32)
                                                                                   tctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAg                   <  1:6424250/96‑1 (MQ=2)
                                                                                   tctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAg                   <  1:3633632/96‑1 (MQ=2)
                                                                                   tctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAg                   <  1:7199432/96‑1 (MQ=2)
                                                                                    ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                  <  1:14000799/101‑1 (MQ=32)
                                                                                    ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                  >  1:6627991/1‑101 (MQ=32)
                                                                                    ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                  <  1:11780201/101‑1 (MQ=32)
                                                                                    ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                  >  1:10896206/1‑101 (MQ=32)
                                                                                    ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                  <  1:9747902/101‑1 (MQ=32)
                                                                                     tctcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGAt                 <  1:6112270/101‑1 (MQ=32)
                                                                                      ctcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATg                <  1:11286160/101‑1 (MQ=32)
                                                                                      cgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATg                <  1:6368232/99‑1 (MQ=2)
                                                                                       tcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATgg               <  1:4860957/101‑1 (MQ=32)
                                                                                       gcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATgg               >  1:2413898/2‑101 (MQ=2)
                                                                                        cACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGa              <  1:14935360/101‑1 (MQ=2)
                                                                                        cACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGa              >  1:213703/1‑101 (MQ=32)
                                                                                        cACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGa              <  1:4649622/101‑1 (MQ=32)
                                                                                         aCGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGAc             <  1:6647398/101‑1 (MQ=2)
                                                                                         aCGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGAc             <  1:4176891/101‑1 (MQ=2)
                                                                                         aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc             >  1:8146370/1‑101 (MQ=32)
                                                                                         aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc             <  1:7557719/101‑1 (MQ=32)
                                                                                         aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc             >  1:4011159/1‑101 (MQ=32)
                                                                                          cGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACg            <  1:9330761/101‑1 (MQ=2)
                                                                                          cGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACg            >  1:6833539/1‑101 (MQ=32)
                                                                                          cGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACg            <  1:1604695/101‑1 (MQ=32)
                                                                                            ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc          <  1:14583839/101‑1 (MQ=32)
                                                                                            ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc          >  1:9242016/1‑101 (MQ=32)
                                                                                            ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc          <  1:419487/101‑1 (MQ=32)
                                                                                            ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc          >  1:12726102/1‑101 (MQ=32)
                                                                                             tGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcg         >  1:9836238/1‑101 (MQ=32)
                                                                                             tGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcg         <  1:2354616/101‑1 (MQ=33)
                                                                                             tGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcg         >  1:5812396/1‑101 (MQ=32)
                                                                                              gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc        <  1:5305805/101‑1 (MQ=33)
                                                                                              gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc        <  1:13698253/101‑1 (MQ=33)
                                                                                              gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc        <  1:14875435/101‑1 (MQ=33)
                                                                                                aTACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGa      <  1:7132742/101‑1 (MQ=32)
                                                                                                aTACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGa      <  1:9742270/101‑1 (MQ=33)
                                                                                                 tACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGAt     <  1:2873872/101‑1 (MQ=33)
                                                                                                  aCGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATc    >  1:5043964/1‑101 (MQ=32)
                                                                                                  aCGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATc    >  1:8941787/1‑101 (MQ=32)
                                                                                                   cGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATcc   >  1:7634334/1‑101 (MQ=32)
                                                                                                    gACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt  >  1:7505216/1‑101 (MQ=32)
                                                                                                    gACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt  <  1:11323146/101‑1 (MQ=33)
                                                                                                    gACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt  <  1:4242597/101‑1 (MQ=33)
                                                                                                    gACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt  >  1:3948897/1‑101 (MQ=32)
                                                                                                    gACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt  <  1:6553429/101‑1 (MQ=33)
                                                                                                    |                                                                                                    
TAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCT  >  NC_007779/4169197‑4169397

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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.