Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,673,784 C→A 15.5% intergenic (+16/+32) yphA → / ← yphB DoxX family inner membrane protein/mutarotase superfamily protein, YphB family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,673,7840CA15.5% 332.1 / 41.4 116intergenic (+16/+32)yphA/yphBDoxX family inner membrane protein/mutarotase superfamily protein, YphB family
Reads supporting (aligned to +/- strand):  ref base C (47/51);  new base A (6/12);  total (53/63)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 121 total aligned reads displayed.

TCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGAT  >  NC_000913/2673686‑2673883
                                                                                                  |                                                                                                   
tCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcg                                                                                                    >  2:1792945/1‑100 (MQ=255)
   ggAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcggcg                                                                                                 >  2:1403979/1‑100 (MQ=255)
   ggAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcggcg                                                                                                 >  2:1554417/1‑100 (MQ=255)
    ggagaaTGTCAGTGTCGCTGGCGCGTTCTTGCTATTGGCAATGACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGt                                                                                                >  2:1751448/3‑100 (MQ=255)
    gaagaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGt                                                                                                >  1:17497/1‑100 (MQ=255)
     aagaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTg                                                                                               <  1:2247064/100‑1 (MQ=255)
       gaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGaa                                                                                             >  1:927892/1‑100 (MQ=255)
         aTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATTCGGCGTGAACg                                                                                           >  1:2470230/1‑100 (MQ=255)
         aTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACg                                                                                           <  2:1723776/100‑1 (MQ=255)
           gTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                         <  1:605655/100‑1 (MQ=255)
           gTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                         <  2:307355/100‑1 (MQ=255)
           gTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                         >  1:1283254/1‑100 (MQ=255)
           gTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                         <  2:513560/100‑1 (MQ=255)
               gTATCGCTGGCGCGTTCTTGCTATTGGCAATTCCCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTAt                                                                                     <  2:138076/100‑1 (MQ=255)
                  tCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                  <  1:1187778/100‑1 (MQ=255)
                  tCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                  <  1:1188069/100‑1 (MQ=255)
                  tCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                  <  2:491846/100‑1 (MQ=255)
                   cGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                 <  1:1527561/100‑1 (MQ=37)
                    gCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                <  1:2169623/100‑1 (MQ=255)
                    gCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                <  2:1589286/100‑1 (MQ=37)
                     cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                               >  1:801337/1‑100 (MQ=255)
                     cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                               <  2:269722/100‑1 (MQ=255)
                     cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                               <  2:2107150/100‑1 (MQ=255)
                     cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                               <  2:946900/100‑1 (MQ=255)
                     cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                               <  2:1192762/100‑1 (MQ=25)
                     cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                               <  2:1230965/100‑1 (MQ=25)
                      tGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                               <  2:1510205/99‑1 (MQ=255)
                      tGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                              <  2:2139717/100‑1 (MQ=25)
                      tGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                              <  2:1503070/100‑1 (MQ=25)
                       ggCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCAtt                                                                             >  1:1084174/1‑100 (MQ=255)
                       ggCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCAtt                                                                             <  2:84035/100‑1 (MQ=255)
                       ggCGCGTTCTTGCTATTGGCAATTACCGCGCCCTCGGCAATTTTTCTGTATCGTCGTCAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                              <  1:2285306/99‑1 (MQ=12)
                         cgcgTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAtctt                                                                           <  2:469291/100‑4 (MQ=21)
                          gcgTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaa                                                                          <  2:1960352/100‑1 (MQ=255)
                           cgTTCTTGCTATTGGCAATTCCCGAACCGGGGGCGATTTCTCTCGATCGGCGTTACGAAAGATTCCGGATAAGGCTTTCACGCCGCATCCGGCa                                                                               <  1:2006593/94‑1 (MQ=16)
                            gTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaata                                                                        >  1:1866357/1‑100 (MQ=255)
                               cTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAAt                                                                     >  2:1886372/1‑100 (MQ=255)
                                ttGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAtctttcctaac                                                                    <  1:306401/100‑11 (MQ=18)
                                   cTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACa                                                                 >  1:2410376/1‑100 (MQ=255)
                                   cTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACa                                                                 <  2:16816/100‑1 (MQ=255)
                                    tATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAg                                                                <  1:406896/100‑1 (MQ=255)
                                     aTTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGc                                                               >  2:1083381/1‑100 (MQ=255)
                                     aTTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGc                                                               <  2:806270/100‑1 (MQ=255)
                                       tGGCGATTACCGGGCAGGGGGTAATTTCTCTCGAGCGGCGTTAGGAAAGATGATGGATGCGGCGTGAACGCCTTATCCGGTATTTAATAAAGTACAg                                                                <  1:153136/97‑1 (MQ=37)
                                        ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                            >  2:1334588/1‑100 (MQ=255)
                                         gCAATTACCGCGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGCACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                            >  1:614298/1‑99 (MQ=255)
                                         gCAATTAACGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAct                                                           >  2:1362264/1‑100 (MQ=255)
                                           aaTTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCActct                                                         >  1:1562960/1‑100 (MQ=255)
                                            aTTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTa                                                        <  2:238367/100‑1 (MQ=255)
                                              tACCGGGCCGGGGGCAATTTCTCTCGATCTGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                      <  1:1554408/100‑1 (MQ=255)
                                              tACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                      <  1:2024491/100‑1 (MQ=255)
                                              tACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                      <  1:748740/100‑1 (MQ=255)
                                               aCCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATCTAATAAATTACAGCCACTCTAccc                                                     >  1:1811459/1‑100 (MQ=255)
                                                ccgggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCg                                                    <  2:112135/100‑1 (MQ=255)
                                                 cgggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCATTTAATAAATTACATCCACTCTAcccac                                                   <  2:1743253/100‑3 (MQ=255)
                                                    gccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAAc                                                <  1:2103956/100‑1 (MQ=255)
                                                    gccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAAc                                                <  2:328401/100‑1 (MQ=255)
                                                         gggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAt                                           >  2:1890553/1‑100 (MQ=255)
                                                          ggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtc                                          >  1:395230/1‑100 (MQ=255)
                                                            cAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtctc                                        >  1:2097985/1‑100 (MQ=255)
                                                            cAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtctc                                        >  2:1265579/1‑100 (MQ=255)
                                                              aTTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTg                                      <  1:1214839/100‑1 (MQ=255)
                                                                ttCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGa                                     <  1:2422140/99‑1 (MQ=255)
                                                                 tctctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAg                                   <  2:559584/100‑1 (MQ=255)
                                                                 tctctcGATCGGCGTTAGGAAAGATACCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAg                                   >  2:1280194/1‑100 (MQ=255)
                                                                  ctctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgt                                  >  2:833615/1‑100 (MQ=255)
                                                                  ctctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCAGCAACGACATCTCTGAAgt                                  >  2:109968/1‑100 (MQ=255)
                                                                   tctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgtt                                 >  2:880186/1‑100 (MQ=255)
                                                                   tctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgtt                                 >  1:970853/1‑100 (MQ=255)
                                                                   tctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgtt                                 >  2:620498/1‑100 (MQ=255)
                                                                    ctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttg                                >  1:1647814/1‑100 (MQ=255)
                                                                    ctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttg                                <  1:23060/100‑1 (MQ=255)
                                                                    ctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttg                                <  2:2312032/100‑1 (MQ=255)
                                                                      cGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttgtt                              >  1:2248210/1‑100 (MQ=255)
                                                                      cGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttgtt                              <  2:116401/100‑1 (MQ=255)
                                                                       gATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGttgttg                             >  2:1833517/1‑100 (MQ=255)
                                                                        aTCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGa                            >  1:2378930/1‑100 (MQ=255)
                                                                        aTCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGa                            <  1:912158/100‑1 (MQ=255)
                                                                         ttggCGGTAGGAAAGAGGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGAt                           <  2:752362/98‑1 (MQ=255)
                                                                             cGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTccc                       <  1:1300274/100‑1 (MQ=255)
                                                                                 aGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTACTAAATTACAGCCACGCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTgg                   >  2:624744/1‑100 (MQ=255)
                                                                                 aGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGCTTCCCCt                     >  2:2567491/1‑98 (MQ=255)
                                                                                  ggAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGt                  >  1:2385096/1‑100 (MQ=255)
                                                                                   gAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTg                 >  1:185015/1‑100 (MQ=255)
                                                                                   attaaaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGt                  >  1:1339714/6‑99 (MQ=21)
                                                                                     aaGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCt               <  2:2062333/100‑1 (MQ=255)
                                                                                       gATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAg             <  1:123183/100‑1 (MQ=255)
                                                                                       aaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAg             >  2:2075235/2‑100 (MQ=25)
                                                                                        aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc            <  1:1032890/100‑1 (MQ=255)
                                                                                        aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc            >  2:114052/1‑100 (MQ=25)
                                                                                        aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTCAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc            >  2:529475/1‑100 (MQ=25)
                                                                                         tGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgcg           <  2:581814/100‑1 (MQ=255)
                                                                                          gCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCTCTGGTGCTAgcgc          <  2:1436911/100‑1 (MQ=255)
                                                                                          gCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgcgc          >  2:1893513/1‑100 (MQ=255)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgcgc          >  2:24913/1‑100 (MQ=37)
                                                                                           ccGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCACCTGGTGCTAGCGCa         >  2:2552552/1‑100 (MQ=255)
                                                                                            cGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCaa        <  1:1403971/100‑1 (MQ=255)
                                                                                             ggATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAAt       >  1:1870149/1‑100 (MQ=255)
                                                                                                  cggcgTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGAt  >  2:1681465/1‑100 (MQ=255)
                                                                                                  cggcgTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGAt  >  2:2283716/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGAT  >  NC_000913/2673686‑2673883

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