Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,113,544 G→A 15.6% intergenic (+138/‑35) opgH → / → yceK OPG biosynthetic ACP‑dependent transmembrane UDP‑glucose beta‑1,2 glycosyltransferase; nutrient‑dependent cell size regulator, FtsZ assembly antagonist/outer membrane integrity lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,113,5440GA15.6% 341.7 / 43.0 140intergenic (+138/‑35)opgH/yceKOPG biosynthetic ACP‑dependent transmembrane UDP‑glucose beta‑1,2 glycosyltransferase; nutrient‑dependent cell size regulator, FtsZ assembly antagonist/outer membrane integrity lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (45/72);  new base A (7/15);  total (52/88)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 142 total aligned reads displayed.

CATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAG  >  NC_000913/1113445‑1113640
                                                                                                   |                                                                                                
cATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCg                                                                                                  >  2:96948/1‑100 (MQ=255)
 aTCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGt                                                                                                 >  2:2205419/1‑100 (MQ=255)
   cGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTa                                                                                               <  1:135161/100‑1 (MQ=255)
    gTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTAt                                                                                              <  1:821041/100‑1 (MQ=255)
       ccTGCAATTTATTGATTCTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGcc                                                                                           >  1:1707300/1‑100 (MQ=255)
         tGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCgg                                                                                         <  1:1054518/100‑1 (MQ=255)
         tGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCgg                                                                                         <  1:1870891/100‑1 (MQ=255)
          gCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGt                                                                                        <  1:123809/100‑1 (MQ=255)
           cAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTa                                                                                       >  2:334666/1‑100 (MQ=255)
             aTTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtt                                                                                     <  2:1975597/100‑1 (MQ=255)
              tttATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAttt                                                                                    <  1:1742663/100‑1 (MQ=25)
               ttATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttt                                                                                   <  1:2091055/100‑1 (MQ=37)
                tATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCGCGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttt                                                                                   <  1:1000556/99‑1 (MQ=37)
                 aTTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAtttttt                                                                                 >  1:1990421/1‑100 (MQ=255)
                 aTTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAtttttt                                                                                 <  1:1743310/100‑1 (MQ=255)
                 aTTGATTTTGCAAGACTTTGTGGGTCAGATAAGGCGTTCACGCCGCATCCGGCAAAAAAAACTAGCCAATAAAAATACCGGCATAACGCCGGTAtttt                                                                                   <  2:1298064/98‑1 (MQ=25)
                  ttGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttttat                                                                                <  2:1534909/100‑3 (MQ=37)
                  ttGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAttttt                                                                                  <  1:1537361/98‑1 (MQ=37)
                   tGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTa                                                                               <  1:401933/100‑1 (MQ=255)
                   tGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTa                                                                               <  1:2326731/100‑1 (MQ=255)
                   tGATTTTGCAAGACTTTGTGGGGCAGATACGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAttt                                                                                    <  2:779750/95‑1 (MQ=25)
                   tGATCTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTCTTTTTTTa                                                                               >  2:1115442/1‑100 (MQ=255)
                    gATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCATGCCGCATCCGGCAAAAACATCGAGCCAATAAAAATACCGGCATAACGCCGGTAtttttatt                                                                               <  2:586622/99‑4 (MQ=37)
                    gATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTAc                                                                              <  1:429435/100‑1 (MQ=255)
                    gATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTAc                                                                              <  2:1215172/100‑1 (MQ=255)
                    gATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTAc                                                                              >  1:2314726/1‑100 (MQ=255)
                    gATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttttattg                                                                              <  1:2421893/100‑5 (MQ=25)
                    gATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttttattg                                                                              <  1:1576332/100‑5 (MQ=25)
                     aTTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACg                                                                             >  2:454524/1‑100 (MQ=255)
                     aTTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttttattgg                                                                             <  2:419542/100‑6 (MQ=25)
                      ttttGCAAGACTTTGTGTGTCGGATAAGGCGTAAACGCCGCATCCGGCAAAAATAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTa                                                                               <  2:1714355/97‑1 (MQ=38)
                       tttGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCTTTATGCCGGTATTTTTTTACg                                                                             <  2:1159369/98‑1 (MQ=255)
                          gCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAaga                                                                        >  1:1940245/1‑100 (MQ=255)
                           cAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAagag                                                                       >  2:1157433/1‑100 (MQ=255)
                           cAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAagag                                                                       <  1:1761418/100‑1 (MQ=255)
                            aaGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAgg                                                                      <  1:1774684/100‑1 (MQ=255)
                              gACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTa                                                                    <  1:1161587/100‑1 (MQ=255)
                                cTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATc                                                                  >  2:1445897/1‑100 (MQ=255)
                                cTTTGGGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAACAAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATc                                                                  <  2:1591481/100‑1 (MQ=255)
                                 tttGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCa                                                                 <  1:91900/100‑1 (MQ=255)
                                  ttGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCaa                                                                <  1:1445908/100‑1 (MQ=255)
                                       ggTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCg                                                           >  2:1715614/1‑100 (MQ=255)
                                       ggTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCg                                                           <  2:1944559/100‑1 (MQ=255)
                                         tCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGtt                                                         <  1:1726318/100‑1 (MQ=255)
                                          cGGATAAGGCGTTCACGCCGCATCCGGAAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGtt                                                         >  2:1073517/1‑99 (MQ=255)
                                           ggATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTa                                                       <  1:887534/100‑1 (MQ=255)
                                               aaGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtg                                                   <  1:293838/100‑1 (MQ=255)
                                                aGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtgt                                                  >  1:1565431/1‑100 (MQ=255)
                                                  gCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtgg                                                <  1:1501639/100‑1 (MQ=255)
                                                  gCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtgg                                                <  2:344966/100‑1 (MQ=255)
                                                  gCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtgg                                                >  1:2457035/1‑100 (MQ=255)
                                                  gCGTTCACGCCGCATCCGGAAAAAACAACGAGCCAATAAAAATACCGGCATTATGCCTGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtgg                                                <  2:743148/100‑1 (MQ=255)
                                                   cGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtggt                                               >  2:1302996/1‑100 (MQ=255)
                                                    gTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtggtg                                              >  1:538512/1‑100 (MQ=255)
                                                    gTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtggtg                                              <  2:1834630/100‑1 (MQ=255)
                                                      tCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGGTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtggtg                                              >  2:485349/1‑98 (MQ=255)
                                                        ncgccgcATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGca                                          >  2:1804816/2‑100 (MQ=255)
                                                        aCGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAg                                            <  2:717531/98‑1 (MQ=255)
                                                         cgccgcATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGcat                                         >  1:1226852/1‑100 (MQ=255)
                                                         cgccgcATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGcat                                         <  1:2379195/100‑1 (MQ=255)
                                                         cgccgcATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTGAATTGTGGTGAGcat                                         >  1:991163/1‑100 (MQ=255)
                                                            cgcATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGcatcat                                      <  2:734988/100‑1 (MQ=255)
                                                             gcATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATg                                     <  1:577891/100‑1 (MQ=255)
                                                                 ccGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGa                                 <  1:2361614/100‑1 (MQ=255)
                                                                  cGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGAc                                >  1:1128569/1‑100 (MQ=255)
                                                                  cGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGAc                                <  2:1542832/100‑1 (MQ=255)
                                                                   ggCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGAcc                               <  1:384840/100‑1 (MQ=255)
                                                                   ggCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGAcc                               >  1:193252/1‑100 (MQ=255)
                                                                    gCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCt                              >  2:1047019/1‑100 (MQ=255)
                                                                     cAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCtt                             <  1:2358400/100‑1 (MQ=255)
                                                                     cAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCtt                             <  2:2370030/100‑1 (MQ=255)
                                                                       aaaaCAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGAATCAAATGCGTTTAATTGTTGTGAGCATCATGGTGACCTTGc                           >  2:1949537/1‑100 (MQ=255)
                                                                        aaaCAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCt                          <  1:587277/100‑1 (MQ=255)
                                                                            aacGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGt                      <  2:444413/100‑1 (MQ=255)
                                                                             acGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTg                     <  1:1325333/100‑1 (MQ=255)
                                                                             acGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTg                     <  2:143558/100‑1 (MQ=255)
                                                                             acGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTg                     <  2:1296986/100‑1 (MQ=255)
                                                                               gAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGc                   <  1:288631/100‑1 (MQ=255)
                                                                                 gCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCtg                 <  2:1896918/100‑1 (MQ=255)
                                                                                 gCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCtg                 <  1:1399845/100‑1 (MQ=255)
                                                                                  cgtaaaaaaaTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTACCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCtgt                >  1:1700074/6‑100 (MQ=37)
                                                                                   gtaaaaaaaTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCtgtg               >  1:1862323/5‑100 (MQ=37)
                                                                                    aaTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTgg              >  2:93564/1‑100 (MQ=255)
                                                                                    aaTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGAATCATGGTGACCTTGCTGAGTGGAtgtg               >  2:552125/1‑99 (MQ=255)
                                                                                     aaaaaaaTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgc             >  1:438604/3‑100 (MQ=38)
                                                                                     aTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgc             <  2:2409082/100‑1 (MQ=255)
                                                                                     aTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgc             <  2:337725/100‑1 (MQ=255)
                                                                                     aTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgc             <  1:1058107/100‑1 (MQ=255)
                                                                                     aTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgc             <  2:2032298/100‑1 (MQ=255)
                                                                                      tAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgca            >  1:1853822/1‑100 (MQ=255)
                                                                                      tAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgca            >  2:1223443/1‑100 (MQ=255)
                                                                                        aaaaTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgcagc          <  1:649944/100‑1 (MQ=255)
                                                                                        aaaaTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgcagc          >  1:973363/1‑100 (MQ=255)
                                                                                        aaaaTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgcagc          <  2:1595239/100‑1 (MQ=255)
                                                                                          aaTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCat        <  1:2126844/100‑1 (MQ=255)
                                                                                          aaTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGAGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgcagca         >  1:2137293/1‑99 (MQ=38)
                                                                                             aCCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCattat     <  2:794531/100‑1 (MQ=255)
                                                                                              ccGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCattatt    >  1:678272/1‑100 (MQ=255)
                                                                                              ccGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCattatt    <  2:155872/100‑1 (MQ=255)
                                                                                               cGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCAttatta   >  1:2272042/1‑100 (MQ=255)
                                                                                                   |                                                                                                
CATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAG  >  NC_000913/1113445‑1113640

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