Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,717,840 G→C 4.5% A339A (GCG→GCC)  atpA → F1 sector of membrane‑bound ATP synthase subunit alpha

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,717,8400G→C 4.5% 23.5 ‑0.3 ‑0.1243A339A (GCG→GCC) atpAF1 sector of membrane‑bound ATP synthase subunit alpha
Reads supporting (aligned to +/- strand):  new base (6/5):  ref base (115/117):  total (121/122)
Fisher's exact test strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.89e-01

Only 100 of 280 total aligned reads displayed.

CCGAATACGTTGAAGCCTTCACCAAAGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAACCCGGGTATTT  >  NC_007779/3717740‑3717940
                                                                                                    |                                                                                                    
ccGAATACGTTGAAGCCTTCACCAAAGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGcg                                                                                                      >  1:5398030/1‑101 (MQ=255)
 cGAATACGTTGAAGCCTTCACCAAAGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGcgt                                                                                                     >  1:13407223/1‑101 (MQ=255)
  gAATACGTTGAAGCCTTCACCAAAGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGcgtt                                                                                                    >  1:1269200/1‑101 (MQ=255)
   aaTACGTTGAAGCCTTCACCAAAGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGcgttc                                                                                                   >  1:10367021/1‑101 (MQ=255)
     tACGTTGAAGCCTTCACCAAAGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGcgttcgt                                                                                                 >  1:8245007/1‑101 (MQ=255)
     tACGTTGAAGCCTTCACCAAAGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGcgttcgt                                                                                                 >  1:11358673/1‑101 (MQ=255)
        gTTGAAGCCTTCACCAAAGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTcc                                                                                              <  1:11653015/101‑1 (MQ=255)
          tGAAGCCTTCACCAAAGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGa                                                                                            >  1:10252422/1‑101 (MQ=255)
            aaGCCTTCACCAAAGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGGTGATTATCGAAACTCAGGCGGGTGATGTTTCTGAGTTCGTTCCCAcc                                                                                          >  1:13907865/1‑101 (MQ=255)
              gCCTTCACCAAAGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCaa                                                                                        <  1:2058604/101‑1 (MQ=255)
                   cACCAAAGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTaa                                                                                   >  1:363977/1‑101 (MQ=255)
                       aaaGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAAtctc                                                                               <  1:5005431/101‑1 (MQ=255)
                           gtgaagtgaaAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCAtt                                                                           <  1:1853384/101‑1 (MQ=255)
                             gaagtgaaAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTAc                                                                         >  1:2452135/1‑101 (MQ=255)
                              aagtgaaAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTAcc                                                                        >  1:9822558/1‑101 (MQ=255)
                               agtgaaAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCg                                                                       >  1:3285945/1‑101 (MQ=255)
                                gtgaaAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGa                                                                      >  1:12518751/1‑101 (MQ=255)
                                gtgaaAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGa                                                                      >  1:8947187/1‑101 (MQ=255)
                                  gaaAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATg                                                                    <  1:7556541/101‑1 (MQ=255)
                                   aaaGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATgg                                                                   <  1:11581413/101‑1 (MQ=255)
                                    aaGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGt                                                                  <  1:14006636/101‑1 (MQ=255)
                                    aaGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGt                                                                  >  1:6819656/1‑101 (MQ=255)
                                     aGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTc                                                                 <  1:7586650/101‑1 (MQ=255)
                                     aGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTc                                                                 <  1:2531197/101‑1 (MQ=255)
                                      gggAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCa                                                                <  1:3001217/101‑1 (MQ=255)
                                       ggAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAg                                                               <  1:13750446/101‑1 (MQ=255)
                                         aaaaCCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGAt                                                             <  1:10121344/101‑1 (MQ=255)
                                           aaCCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCt                                                           >  1:3710034/1‑101 (MQ=255)
                                             caggctcgctGACCGCACTGCCGATCATCGAGACCAAGGCAAACGACGTGTCGGCCTTCATTCCGACCAACGTCATTTCCATTACCGACGGCCAGGTCTTc                                                         <  1:8901233/93‑1 (MQ=255)
                                                gTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTg                                                      <  1:466154/101‑1 (MQ=255)
                                                 ttCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTgg                                                     <  1:5137875/101‑1 (MQ=255)
                                                 ttCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTgg                                                     <  1:8097466/101‑1 (MQ=255)
                                                 ctcgctGACCGCACTGCCGATCATCGAGACCAAGGCAAACGACGTGTCGGCCTTCATTCCGACCAACGTCATTTCCATTACCGACGGCCAGGTCTTCCtcg                                                     >  1:8608024/5‑99 (MQ=255)
                                                  tctcTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGa                                                    <  1:12948574/101‑1 (MQ=255)
                                                  tcgctGACCGCACTGCCGATCATCGAGACCAAGGCAAACGACGTGTCGGCCTTCATTCCGACCAACGTCATTTCCATTACCGACGGCCAGGTCTTCCTCGa                                                    <  1:1057631/98‑1 (MQ=255)
                                                   ctctGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGaa                                                   >  1:7112431/1‑101 (MQ=255)
                                                   ctctGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGaa                                                   >  1:2675365/1‑101 (MQ=255)
                                                   cgctGACCGCACTGCCGATCATCGAGACCAAGGCAAACGACGTGTCGGCCTTCATTCCGGCCAACGTCATTTCCATTACCGACGGCCAGGTCTTCCTCGag                                                   <  1:6142809/99‑2 (MQ=255)
                                                   cgctGACCGCACTGCCGATCATCGAGACCAAGGCAAACGACGTGTCGGCCTTCATTCCGACCAACGTCATTTCCATTACCGACGGCCAGGTCTTCCTCGag                                                   >  1:6294615/3‑100 (MQ=255)
                                                    gctGACCGCACTGCCGATCATCGAGACCAAGGCAAACGACGTGTCGGCCTTCATTCCGACCAACGTCATTTCCATTACCGACGGCCAGGTCTTCCTCGagt                                                  >  1:3308315/2‑99 (MQ=255)
                                                      tGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAaacc                                                <  1:1662436/101‑1 (MQ=255)
                                                      tGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAaacc                                                >  1:6628506/1‑101 (MQ=255)
                                                      tGACCGCACTGCCGATCATCGAGACCAAGGCAAACGACGTGTCGGCCTTCATTCCGACCAACGTCATTTCCATTACCGACGGCCAGGTCTTCCTCGagtcc                                                >  1:5883006/1‑97 (MQ=255)
                                                       gACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAaacca                                               >  1:7880067/1‑101 (MQ=255)
                                                       gACCGCACTGCCGATCATCGAGACCAAGGCAAACGACGTGTCGGCCTTCATTCCGACCAACGTCATTTCCATTACCGACGGCCAGGTCTTCCTCGagtccg                                               >  1:2028419/1‑96 (MQ=255)
                                                         ccGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAaaccaac                                             >  1:960394/1‑101 (MQ=255)
                                                          cGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAaaccaacc                                            >  1:4713638/1‑101 (MQ=255)
                                                            cACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGCTGGTCAGATCTTCCTGGAAACCAACCTg                                          >  1:13084577/1‑101 (MQ=255)
                                                            cACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTg                                          >  1:12849819/1‑101 (MQ=255)
                                                              cTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGtt                                        >  1:5302721/1‑101 (MQ=255)
                                                              cTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGtt                                        <  1:4149305/101‑1 (MQ=255)
                                                               tGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTc                                       >  1:2617970/1‑101 (MQ=255)
                                                               tGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTc                                       >  1:4785086/1‑101 (MQ=255)
                                                                gCCGATTATCGAAACTCCGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCa                                      >  1:7840969/1‑101 (MQ=255)
                                                                 ccGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCaa                                     <  1:7273259/101‑1 (MQ=255)
                                                                   gATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCTTGGAAACCAACCTGTTCAACg                                   <  1:11917555/101‑1 (MQ=255)
                                                                   gATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACg                                   <  1:741986/101‑1 (MQ=255)
                                                                   gATTATCGAAACTCAGGCGGGCGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACg                                   >  1:8755953/1‑101 (MQ=255)
                                                                      tATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCg                                >  1:9381577/1‑101 (MQ=255)
                                                                      tATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCg                                <  1:5204610/101‑1 (MQ=255)
                                                                         cGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTa                             >  1:10976187/1‑101 (MQ=255)
                                                                          gAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTAt                            >  1:1735871/1‑101 (MQ=255)
                                                                          gAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTAt                            >  1:7730576/1‑101 (MQ=255)
                                                                          gAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGGAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTAt                            >  1:3115905/1‑101 (MQ=255)
                                                                            aaCTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTc                          <  1:10631572/101‑1 (MQ=255)
                                                                            aaCTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTc                          <  1:2706756/101‑1 (MQ=255)
                                                                            aaCTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTc                          >  1:6055026/1‑101 (MQ=255)
                                                                            aaCTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTc                          <  1:12445707/101‑1 (MQ=255)
                                                                              cTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGt                        <  1:2272719/101‑1 (MQ=255)
                                                                              cTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGt                        <  1:7228397/101‑1 (MQ=255)
                                                                                 aGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCt                     >  1:8300373/1‑101 (MQ=255)
                                                                                  ggCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTg                    <  1:1585730/101‑1 (MQ=255)
                                                                                   gCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGc                   <  1:10441120/101‑1 (MQ=255)
                                                                                   gCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGc                   <  1:4517320/101‑1 (MQ=255)
                                                                                   gCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGc                   >  1:1080318/1‑101 (MQ=255)
                                                                                     gggTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCgg                 <  1:11230851/101‑1 (MQ=255)
                                                                                     gggTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCgg                 >  1:10508953/1‑101 (MQ=255)
                                                                                      ggTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGt                >  1:925817/1‑101 (MQ=255)
                                                                                       gTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGtt               <  1:9166367/101‑1 (MQ=255)
                                                                                       gTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGtt               <  1:5330687/101‑1 (MQ=255)
                                                                                        tGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTa              <  1:1114825/101‑1 (MQ=255)
                                                                                         gACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTaa             <  1:13524111/101‑1 (MQ=255)
                                                                                         gACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTaa             <  1:1421744/101‑1 (MQ=255)
                                                                                         gACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTaa             >  1:13155556/1‑101 (MQ=255)
                                                                                         gACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTaa             <  1:12790707/101‑1 (MQ=255)
                                                                                         gACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTaa             >  1:898109/1‑101 (MQ=255)
                                                                                          aCGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAAc            >  1:6969670/1‑101 (MQ=255)
                                                                                          aCGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAAc            >  1:373491/1‑101 (MQ=255)
                                                                                          aCGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAAc            >  1:3050497/1‑101 (MQ=255)
                                                                                            gTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAAccc          <  1:1571869/101‑1 (MQ=255)
                                                                                            gTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAAccc          <  1:1505880/101‑1 (MQ=255)
                                                                                             tttCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAACCCg         >  1:9224236/1‑101 (MQ=255)
                                                                                             tttCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAACCCg         <  1:1100229/101‑1 (MQ=255)
                                                                                              ttCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAACCCgg        <  1:1741046/101‑1 (MQ=255)
                                                                                              ttCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAACCCgg        <  1:12238280/101‑1 (MQ=255)
                                                                                               tCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAACCCggg       <  1:8147770/101‑1 (MQ=255)
                                                                                               tCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAACCCggg       >  1:2875275/1‑101 (MQ=255)
                                                                                                 tGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAACCCGGGTa     >  1:10744890/1‑101 (MQ=255)
                                                                                                    gttcgttcCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAACCCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAACCCGGGTAttt  <  1:10056244/101‑1 (MQ=255)
                                                                                                    gttcgttcCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAACCCGGGTAttt  >  1:5626586/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CCGAATACGTTGAAGCCTTCACCAAAGGTGAAGTGAAAGGGAAAACCGGTTCTCTGACCGCACTGCCGATTATCGAAACTCAGGCGGGTGACGTTTCTGCGTTCGTTCCGACCAACGTAATCTCCATTACCGATGGTCAGATCTTCCTGGAAACCAACCTGTTCAACGCCGGTATTCGTCCTGCGGTTAACCCGGGTATTT  >  NC_007779/3717740‑3717940

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