Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,920,796 T→C 17.2% A31A (GCA→GCG atpF ← F0 sector of membrane‑bound ATP synthase, subunit b

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,920,7960TC17.2% 270.3 / 46.6 122A31A (GCA→GCGatpFF0 sector of membrane‑bound ATP synthase, subunit b
Reads supporting (aligned to +/- strand):  ref base T (54/47);  new base C (8/13);  total (62/60)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 124 total aligned reads displayed.

CTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACA  >  NC_000913/3920697‑3920894
                                                                                                   |                                                                                                  
cTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCt                                                                                                    >  1:58230/1‑100 (MQ=255)
  ggtcggtcGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGc                                                                                                  >  1:2182098/1‑100 (MQ=255)
  ggtcggtcGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGc                                                                                                  >  1:1880865/1‑100 (MQ=255)
  ggtcggtcGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGc                                                                                                  >  2:2528016/1‑100 (MQ=255)
  ggtcggtcGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGc                                                                                                  >  1:582473/1‑100 (MQ=255)
   gtcggtcgCGCTGGGCTTTGCAAGGTCAAGGTACTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCGATTTCTTTTTGACGTTTTTCGATGGCTGc                                                                                                  <  2:1289007/99‑1 (MQ=255)
     cggtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCAt                                                                                               <  2:1458297/100‑1 (MQ=255)
     cggtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCAt                                                                                               <  2:1458059/100‑1 (MQ=255)
     cggtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCAt                                                                                               >  1:2549986/1‑100 (MQ=255)
     cggtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTACCAt                                                                                               <  2:1591419/100‑1 (MQ=255)
      ggtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCAtt                                                                                              >  2:535647/1‑100 (MQ=255)
       gtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTTTTTTTGACGTTTTTCGATGGCTGCCttt                                                                                              >  1:711827/1‑99 (MQ=255)
       gtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTa                                                                                             >  2:2476024/1‑100 (MQ=255)
            gcTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATggc                                                                                        <  2:1570903/100‑1 (MQ=255)
              tGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATggcgg                                                                                      <  2:1443289/100‑1 (MQ=255)
               ggCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATggcggc                                                                                     <  1:1701813/100‑1 (MQ=255)
               ggCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATggcggc                                                                                     <  1:2435422/100‑1 (MQ=255)
               ggCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATAGCTGGCATTAATgg                                                                                         >  1:1853162/1‑96 (MQ=255)
                gCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGcc                                                                                    <  1:2566042/100‑1 (MQ=255)
                 ccTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCa                                                                                   >  2:936089/1‑100 (MQ=255)
                 ccTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCa                                                                                   >  2:662281/1‑100 (MQ=255)
                  cTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCAt                                                                                  <  2:228010/100‑1 (MQ=255)
                   tttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCAtc                                                                                 <  2:1152000/100‑2 (MQ=255)
                   tttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCAtc                                                                                 <  1:878753/100‑2 (MQ=255)
                   tttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCAtc                                                                                 <  2:714266/100‑2 (MQ=255)
                   tttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCAtc                                                                                 <  2:146260/100‑2 (MQ=255)
                   tttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCAtc                                                                                 <  2:714198/100‑2 (MQ=255)
                   tttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCAtc                                                                                 <  2:1681448/100‑2 (MQ=255)
                    ttGCAAGGTCAAGGTCCTTCTGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCAtc                                                                                 <  1:1958918/99‑2 (MQ=255)
                    ttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATAc                                                                                <  2:374332/100‑1 (MQ=255)
                    ttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCAtcg                                                                                <  2:493305/100‑3 (MQ=255)
                    ttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTCTTGACGTTTTTCGATGGCCGCCAGTAATGGCAGCCAtc                                                                                 <  1:1603910/99‑2 (MQ=255)
                      gCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGt                                                                              >  2:160102/1‑100 (MQ=255)
                        aaggtcaaggtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTAc                                                                            >  2:1835287/1‑100 (MQ=255)
                         aggtcaaggtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGCTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACt                                                                           <  1:1385719/100‑1 (MQ=255)
                         aggtcaaggtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTCTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACt                                                                           >  2:1145576/1‑100 (MQ=255)
                           gtcaaggtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTc                                                                         >  1:201224/1‑100 (MQ=255)
                            tccgatctcTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCa                                                                        <  1:56763/92‑1 (MQ=255)
                             caaggtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCAt                                                                       <  2:349730/100‑1 (MQ=255)
                                   ccTTATGTGCTCTTTCTGCGTAAGCAAGGCCGTCCGCAATTTCTTTTTGACGTGTTTCGATGGCTGCCATTTATGTCGGCCATACGCACTTCATGcaga                                                                  <  1:925096/99‑1 (MQ=255)
                                    cTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCaagatcgga                                                                >  2:56764/1‑92 (MQ=255)
                                     ttATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcagaaca                                                               >  2:113505/1‑100 (MQ=255)
                                      tATGTTCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcagaacag                                                              <  2:137107/100‑1 (MQ=255)
                                       aTGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcagaacaga                                                             <  2:2211720/100‑1 (MQ=255)
                                       aTGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcagaacaga                                                             >  1:1189538/1‑100 (MQ=255)
                                        tgtgCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcagaacagaa                                                            >  2:1635732/1‑100 (MQ=255)
                                          tgCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacg                                                          >  2:2239969/1‑100 (MQ=255)
                                           gCTCGTTCTGCGGAAGCAAGGCCGTCAGCCATTGCTTTTTGACGGTGTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacga                                                         >  2:2572499/1‑100 (MQ=255)
                                             tCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacgaac                                                       <  2:2493711/100‑1 (MQ=255)
                                             tCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacgaac                                                       >  2:2148866/1‑100 (MQ=255)
                                             tCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacgaac                                                       <  2:2284121/100‑1 (MQ=255)
                                              cGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACa                                                      <  2:2232107/100‑1 (MQ=255)
                                               gTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAg                                                     <  2:822606/100‑1 (MQ=255)
                                                ttCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAgg                                                    <  2:2058027/100‑1 (MQ=255)
                                                 tCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGa                                                   >  2:682403/1‑100 (MQ=255)
                                                  cTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGAc                                                  <  1:1339926/100‑1 (MQ=255)
                                                  cTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGAc                                                  >  2:1881133/1‑100 (MQ=255)
                                                  cTGCGGAAGCAAGGCCGTCAGCAATTTCTTTATGACGTGTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAg                                                     <  1:642251/97‑1 (MQ=255)
                                                   tGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACa                                                 >  2:2128930/1‑100 (MQ=255)
                                                    gCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACaa                                                >  2:1868049/1‑100 (MQ=255)
                                                      ggAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAAc                                              >  1:2009551/1‑100 (MQ=255)
                                                         aGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAAcgcg                                           >  2:1373623/1‑100 (MQ=255)
                                                           cAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAt                                         >  1:46361/1‑100 (MQ=255)
                                                           cAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAt                                         >  1:1393601/1‑100 (MQ=255)
                                                             aGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtgg                                       >  2:1003590/1‑100 (MQ=255)
                                                               gCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcc                                     <  2:1223561/100‑1 (MQ=255)
                                                                ccGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcct                                    >  1:2458139/1‑100 (MQ=255)
                                                                 cGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcctg                                   >  2:151090/1‑100 (MQ=255)
                                                                  gTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcctgg                                  >  2:965376/1‑100 (MQ=255)
                                                                    cAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcctggcc                                <  2:1189542/100‑1 (MQ=255)
                                                                    cAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcctggcc                                <  1:1565081/100‑1 (MQ=255)
                                                                       cAATTTCTTTTTGACGTTTTTTTATTTTTTTCATTAATGGCGGCCATAAGTACTTGATGCAGAACAGAACGAACAGGACAAACGCGATGGACTGGCAGa                              >  1:498695/1‑99 (MQ=255)
                                                                       cAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAg                             >  2:1552744/1‑100 (MQ=255)
                                                                       cAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAg                             <  2:2176185/100‑1 (MQ=255)
                                                                         aTTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGa                           <  1:1827061/100‑1 (MQ=255)
                                                                             cTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgt                       >  2:1328016/1‑100 (MQ=255)
                                                                             cTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgt                       >  2:848621/1‑100 (MQ=255)
                                                                              tttttGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgtt                      >  2:1875743/1‑100 (MQ=255)
                                                                              tttttGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgtt                      >  1:355973/1‑100 (MQ=255)
                                                                              tttttGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgtt                      >  2:2463372/1‑100 (MQ=255)
                                                                               ttttGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgttg                     <  1:922205/100‑1 (MQ=255)
                                                                                tttGACATTTATCAATGGCTGCCATTAATGGCGGACATACCTACTTCATGCAGAACAGAACGAACAGGACAATCGCGTTGGCCTGGACGAGGATAGATGc                    <  1:637898/100‑1 (MQ=255)
                                                                                 ttGACGTTTTTCGACGGCTGCCATTAATGGCGGGCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgttg                     >  1:1650064/1‑98 (MQ=255)
                                                                                  tGACGTTTTTCGATGGCTGCCATTAATTGCGGCCATACGTACTTCATGCAGCACAGCACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCg                   >  1:933125/1‑99 (MQ=255)
                                                                                   gACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGtt                 <  2:2584567/100‑1 (MQ=255)
                                                                                   gACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGtt                 >  1:231251/1‑100 (MQ=255)
                                                                                    aCGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTa                <  1:1047783/100‑1 (MQ=255)
                                                                                    aCGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTa                >  1:1197655/1‑100 (MQ=255)
                                                                                        gtacgtaTGGCCGCCATTAATGGCAGGCATACGTACTTCATGCCGAACAGAACGACGAGGACCAACGAGATGGGCTGGCCGAGGATTGTTGCGTTAAg              >  1:10802/7‑98 (MQ=25)
                                                                                         tttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGAtt           <  2:699010/100‑1 (MQ=255)
                                                                                          ncgtaTGGCCGCCATTAATGGCAGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTc          >  2:1071356/5‑100 (MQ=39)
                                                                                            gtaTGGCCGCCATTAATGGCAGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTcac        >  1:2352655/3‑100 (MQ=255)
                                                                                            gtaTGGCCGCCATTAATGGCAGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTcac        >  1:1870100/3‑100 (MQ=255)
                                                                                            cGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTcac        <  1:581082/100‑1 (MQ=255)
                                                                                            cGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTcac        <  2:1708214/100‑1 (MQ=255)
                                                                                             taTGGCCGCCATTAATGGCAGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGACGAGGATTGTTGCGTTAAGATTcaca       >  1:7068/2‑100 (MQ=255)
                                                                                             taTGGCCGCCATTAATGGCAGCCAGACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTcaca       >  2:1157044/2‑100 (MQ=255)
                                                                                             gATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTcaca       >  1:1324478/1‑100 (MQ=255)
                                                                                               tGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGc     <  1:1144905/100‑1 (MQ=255)
                                                                                                ggCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGca    <  2:103229/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACA  >  NC_000913/3920697‑3920894

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: