Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,788,194 G→A 5.5% P250P (CCG→CCA ppsR → PEP synthase kinase and PEP synthase pyrophosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,788,1940GA 5.5% 381.8 / 7.3 109P250P (CCG→CCAppsRPEP synthase kinase and PEP synthase pyrophosphorylase
Reads supporting (aligned to +/- strand):  ref base G (52/51);  new base A (2/4);  total (54/55)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

CAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACC  >  NC_000913/1788045‑1788342
                                                                                                                                                     |                                                                                                                                                     
cAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCg                                                                                                                                                       <  2:2175498/150‑1 (MQ=255)
 atcaaaaaTTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCATCTCATCTTCAGTGCAGGATGGAAGTCGAGGAAGTGGAAGCCTTGTAACGTAAAAATCAGATCACga                                                                                                                                                      >  1:943615/6‑149 (MQ=255)
 aGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGt                                                                                                                                                      >  1:943582/1‑150 (MQ=255)
 aGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGt                                                                                                                                                      >  1:743335/1‑150 (MQ=255)
   cATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGa                                                                                                                                                   >  2:318823/1‑150 (MQ=255)
     tAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCATGGAt                                                                                                                                                  >  1:939737/1‑150 (MQ=255)
     tAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCATGGAt                                                                                                                                                  >  1:2128134/1‑150 (MQ=255)
        aTTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTaa                                                                                                                                               <  1:445888/150‑1 (MQ=255)
        aTTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTaa                                                                                                                                               <  1:2100739/150‑1 (MQ=255)
        aTTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTaa                                                                                                                                               <  1:1603543/150‑1 (MQ=255)
        aTTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTaa                                                                                                                                               <  1:144378/150‑1 (MQ=255)
         ttGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAAc                                                                                                                                              >  2:1096259/1‑150 (MQ=255)
              cGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTAc                                                                                                                                         <  1:1079254/150‑1 (MQ=255)
               ggCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTAcc                                                                                                                                        >  1:1257870/1‑150 (MQ=255)
               ggCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTAcc                                                                                                                                        >  1:284457/1‑150 (MQ=255)
               ggCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTAcc                                                                                                                                        >  1:230616/1‑150 (MQ=255)
                gCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCa                                                                                                                                       >  1:1264648/1‑150 (MQ=255)
                    gACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAatta                                                                                                                                   >  1:1798737/1‑150 (MQ=255)
                    gACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAatta                                                                                                                                   >  1:1798740/1‑150 (MQ=255)
                         tCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCgg                                                                                                                              <  2:717400/150‑1 (MQ=255)
                           gACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTa                                                                                                                            >  2:1312631/1‑150 (MQ=255)
                           gACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTa                                                                                                                            >  1:2227861/1‑150 (MQ=255)
                            aCCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTag                                                                                                                           >  2:1890856/1‑150 (MQ=255)
                              ccGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTagaa                                                                                                                         <  2:1061360/150‑1 (MQ=255)
                              ccGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTagaa                                                                                                                         <  2:527598/150‑1 (MQ=255)
                                 gAACGTCTGGCGGCGATTCGCGAGGAACGGCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAagag                                                                                                                      >  2:2207377/1‑150 (MQ=255)
                                    cGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGAtt                                                                                                                   <  2:1308167/150‑1 (MQ=255)
                                    cGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGAtt                                                                                                                   <  1:1686548/150‑1 (MQ=255)
                                     gTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTg                                                                                                                  >  1:1143757/1‑150 (MQ=255)
                                     gTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTg                                                                                                                  >  1:1652206/1‑150 (MQ=255)
                                        tGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGcca                                                                                                               >  1:799317/1‑150 (MQ=255)
                                            ggcgATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCaa                                                                                                           <  2:1143342/150‑1 (MQ=255)
                                             gcgATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAg                                                                                                          >  2:863850/1‑150 (MQ=255)
                                             gcgATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAg                                                                                                          >  2:857005/1‑150 (MQ=255)
                                               gATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGAt                                                                                                        <  1:532374/150‑1 (MQ=255)
                                                   cgcgAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAACTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTc                                                                                                    <  1:1948826/150‑1 (MQ=255)
                                                    gcgAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTACAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGGAGAAGAGATTGCCACCAAGATCCTCg                                                                                                   >  2:611142/1‑150 (MQ=255)
                                                    gcgAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCg                                                                                                   >  1:1285803/1‑150 (MQ=255)
                                                      gAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGag                                                                                                 >  2:1464357/1‑149 (MQ=255)
                                                      gAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCGCGat                                                                                                 >  1:654061/1‑150 (MQ=255)
                                                       aGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGGAGAAGAGATTGCCACCAAGGTCCTCgtga                                                                                                >  1:896368/1‑147 (MQ=255)
                                                         gAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATc                                                                                              >  2:286295/1‑150 (MQ=255)
                                                               cGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGc                                                                                        <  2:2225545/150‑1 (MQ=255)
                                                               cGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGc                                                                                        <  2:1380805/150‑1 (MQ=255)
                                                                  gagaACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCtt                                                                                     >  1:1219486/1‑150 (MQ=255)
                                                                      aCAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTc                                                                                 >  1:1559055/1‑150 (MQ=255)
                                                                        aGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGc                                                                               <  2:1911863/150‑1 (MQ=255)
                                                                         gTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGcc                                                                              <  1:588177/150‑1 (MQ=255)
                                                                         gTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGcc                                                                              <  1:249516/150‑1 (MQ=255)
                                                                            gCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGaa                                                                           <  1:1274805/150‑1 (MQ=255)
                                                                                 gCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTAc                                                                      <  1:702705/150‑1 (MQ=255)
                                                                                    tCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTag                                                                   >  2:2234857/1‑150 (MQ=255)
                                                                                     cGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTaga                                                                  <  1:1312631/150‑1 (MQ=255)
                                                                                     cGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTaga                                                                  <  2:683566/150‑1 (MQ=255)
                                                                                     cGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTaga                                                                  <  2:297707/150‑1 (MQ=255)
                                                                                        ttCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGaa                                                               >  2:1352944/1‑150 (MQ=255)
                                                                                         tCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTagttctc                                                              <  2:1723780/150‑6 (MQ=255)
                                                                                         tCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTagttctc                                                              <  1:1528480/150‑6 (MQ=255)
                                                                                         tCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAAc                                                              <  2:1915621/150‑1 (MQ=255)
                                                                                            tCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAg                                                           <  1:126605/150‑1 (MQ=255)
                                                                                            tCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAg                                                           <  2:768771/150‑1 (MQ=255)
                                                                                              aGTGCAGGATGGAAGTCGCGGAAGTGGTAGCCTTGTACCGTAAAAATCAGATCCCATGGATTTACTGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGATTGTACTAGAGAACTAGTg                                                         <  2:922760/150‑1 (MQ=255)
                                                                                              aGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTg                                                         >  2:1793754/1‑150 (MQ=255)
                                                                                              aGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCATGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTg                                                         <  2:2111279/150‑1 (MQ=255)
                                                                                                tGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGACGAGATTGCCACCAAGATGCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGACCTAGTGCa                                                       >  1:2056201/1‑150 (MQ=255)
                                                                                                ggCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCa                                                       <  1:1265515/149‑1 (MQ=255)
                                                                                                  cAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCAtt                                                     <  1:1582729/150‑1 (MQ=255)
                                                                                                  cAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCAtt                                                     >  1:1635021/1‑150 (MQ=255)
                                                                                                     gATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGc                                                  <  1:1482388/150‑1 (MQ=255)
                                                                                                      aTGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCt                                                 <  1:1340826/150‑1 (MQ=255)
                                                                                                      aTGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCt                                                 >  2:1647647/1‑150 (MQ=255)
                                                                                                         gAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTCGTGCATTAGCTTAt                                             >  2:1818094/1‑150 (MQ=255)
                                                                                                          aaGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTAtt                                            <  2:1486208/150‑1 (MQ=255)
                                                                                                           aGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTAttt                                           <  1:318823/150‑1 (MQ=255)
                                                                                                             tCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATGTCGTGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTAttttt                                         >  1:174509/1‑150 (MQ=255)
                                                                                                             tCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTAttttt                                         >  1:715883/1‑150 (MQ=255)
                                                                                                                   aaGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTa                                   <  2:799317/150‑1 (MQ=255)
                                                                                                                     gTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGTTCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATc                                 >  2:1485301/1‑150 (MQ=255)
                                                                                                                     gTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACCAGTGCATTAGCTTATTTTTTTGTTATc                                 >  2:1485222/1‑150 (MQ=255)
                                                                                                                         aaGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGc                             >  2:65614/1‑150 (MQ=255)
                                                                                                                             cTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGCTATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTTGTTATCATGCTaa                          >  2:746927/1‑150 (MQ=255)
                                                                                                                             cTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAac                         >  2:2255825/1‑150 (MQ=255)
                                                                                                                             cTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAac                         >  1:197276/1‑150 (MQ=255)
                                                                                                                             cTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAac                         <  1:1064399/150‑1 (MQ=255)
                                                                                                                             cTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATAATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTTTTTATAATGCTaa                          >  2:1800533/1‑150 (MQ=255)
                                                                                                                                gTACCGTAAAAATCAGATCCCATGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAaccac                      <  2:2128134/150‑1 (MQ=255)
                                                                                                                                  aCCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCAccc                    <  2:1798740/150‑1 (MQ=255)
                                                                                                                                  aCCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCAccc                    <  2:1798737/150‑1 (MQ=255)
                                                                                                                                           aTCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGgtg           >  2:863499/1‑150 (MQ=255)
                                                                                                                                           aTCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGgtg           <  2:1219486/150‑1 (MQ=255)
                                                                                                                                           aTCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGAGCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTCGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGgtg           >  2:863473/1‑150 (MQ=255)
                                                                                                                                            tCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGTTCGTCGTTATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCTTGCTAACCACCCGGCGAGgtgt          >  1:2107736/1‑150 (MQ=255)
                                                                                                                                            tCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGgtgt          >  2:529659/1‑150 (MQ=255)
                                                                                                                                              agagCCCGTGGATTAACAGTACCACTTATTCGGTAGAAGAGATTGCCACCAAGCTCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGa        <  1:611142/150‑1 (MQ=255)
                                                                                                                                              aGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGCTATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGa        >  2:523294/1‑150 (MQ=255)
                                                                                                                                               gATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGac       <  2:230616/150‑1 (MQ=255)
                                                                                                                                               gATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGac       <  2:1276477/150‑1 (MQ=255)
                                                                                                                                                  cccGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGacaca    <  1:117909/150‑1 (MQ=255)
                                                                                                                                                   ccGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGacacac   <  1:1647647/150‑1 (MQ=255)
                                                                                                                                                    cGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACAcc  <  1:929836/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACC  >  NC_000913/1788045‑1788342

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