Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,843,115 G→A 5.2% intergenic (+71/+33) PP_1642 → / ← PP_1643 hypothetical protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,843,1150GA 5.2% 866.5 / 22.3 266intergenic (+71/+33)PP_1642/PP_1643hypothetical protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (160/92);  new base A (7/7);  total (167/99)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 268 total aligned reads displayed.

GTTCAACCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCC  >  NC_002947/1842966‑1843263
                                                                                                                                                     |                                                                                                                                                    
gTTCAACCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGCGCAg                                                                                                                                                      >  2:1749425/1‑150 (MQ=255)
 ttCAACCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGt                                                                                                                                                     >  2:3467341/1‑150 (MQ=255)
   cAACCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTAc                                                                                                                                                   <  2:3116998/150‑1 (MQ=255)
    aaCCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTAcc                                                                                                                                                  >  2:6795892/1‑150 (MQ=255)
    aaCCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTAcc                                                                                                                                                  <  1:3453471/150‑1 (MQ=255)
     aCCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCt                                                                                                                                                 <  2:736067/150‑1 (MQ=255)
       cAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTAAAGGCCTGTGCAGTACCtgt                                                                                                                                               >  1:3631436/1‑150 (MQ=255)
         gTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTgg                                                                                                                                             <  1:2647780/150‑1 (MQ=255)
            gCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCCACCATTGAGCCTGGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCCCTTAAGGCCTGGGCCGTACCTGTGGGAg                                                                                                                                          >  2:5743265/1‑150 (MQ=255)
             cATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGc                                                                                                                                         <  1:4593136/150‑1 (MQ=18)
               tGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCgg                                                                                                                                       <  2:3703800/150‑1 (MQ=255)
                  agagagCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCg                                                                                                                                    >  2:4280219/1‑150 (MQ=255)
                      agCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGcc                                                                                                                                >  1:5807959/1‑150 (MQ=255)
                       gcgggcggCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGccc                                                                                                                               <  1:3694810/150‑1 (MQ=255)
                       gcgggcggCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGCGGGCATGccc                                                                                                                               <  2:5512125/150‑1 (MQ=17)
                               ccGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGaaaaa                                                                                                                       >  2:1897345/1‑150 (MQ=255)
                                    gtgCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTaa                                                                                                                  <  2:1550666/150‑1 (MQ=255)
                                     tgCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAAc                                                                                                                 >  2:654918/1‑150 (MQ=255)
                                            tGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTc                                                                                                          <  2:2739310/150‑1 (MQ=255)
                                              tCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTcgc                                                                                                        >  2:586149/1‑150 (MQ=255)
                                                ggTGATGCTGGCGCACCCCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTcgccc                                                                                                      <  2:6295506/150‑1 (MQ=255)
                                                 gTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTcgcccg                                                                                                     >  2:1010637/1‑150 (MQ=255)
                                                 gTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCGGGGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTcgcccg                                                                                                     >  1:103001/1‑150 (MQ=255)
                                                              acacCGTCACCATTTGAGCCTTGTGTTGCTCCCCTTCGCGGGCTTGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGGGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTc                                                                                        <  1:6353283/150‑1 (MQ=38)
                                                               cacCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCa                                                                                       >  1:2152774/1‑150 (MQ=255)
                                                               cacCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCGGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCa                                                                                       >  2:6128930/1‑150 (MQ=255)
                                                                  cGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGc                                                                                    >  1:3165143/1‑150 (MQ=255)
                                                                   gTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCt                                                                                   >  2:2140900/1‑150 (MQ=255)
                                                                    tCAACCTTTGCGCCTTGTGTTGCTCCCTTGCGCGGGCTGGCCCGCTCCCACAGTAAGCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGTGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTg                                                                                  <  1:3097764/150‑1 (MQ=38)
                                                                    ccaccccttGCGCCTTGTGGTGCTCCCCTGCGCGGGCAGGCCCGCTCCCACAGTAAACGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTg                                                                                  <  1:5529289/143‑1 (MQ=38)
                                                                     cAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTgg                                                                                 <  1:2256597/150‑1 (MQ=255)
                                                                      aaCCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGa                                                                                >  1:3509115/1‑150 (MQ=255)
                                                                       aCCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGaa                                                                               >  2:4567157/1‑150 (MQ=255)
                                                                         cATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACt                                                                             >  2:2395820/1‑150 (MQ=255)
                                                                         cATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACt                                                                             <  1:5106676/150‑1 (MQ=255)
                                                                          aTTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACtt                                                                            >  2:4570005/1‑150 (MQ=255)
                                                                           ttGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTc                                                                           <  1:2971376/150‑1 (MQ=255)
                                                                             gAGCCTTGTGTTGCTCCCAGTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCa                                                                         <  2:4726672/150‑1 (MQ=255)
                                                                              aGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAg                                                                        <  1:6007813/150‑1 (MQ=255)
                                                                                ccTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAgcg                                                                      >  1:1514010/1‑150 (MQ=255)
                                                                                  ttGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAATATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCa                                                                    >  1:4143616/1‑150 (MQ=255)
                                                                                     tgtTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGt                                                                 >  1:3870245/1‑150 (MQ=255)
                                                                                      gtTGCTCCCATTCGCGGGCATGCCCGCTCCCCCAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTa                                                                <  2:3330655/150‑1 (MQ=255)
                                                                                      gtTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTa                                                                >  2:6504442/1‑150 (MQ=255)
                                                                                         gCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAg                                                             >  1:200388/1‑150 (MQ=255)
                                                                                           tCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTa                                                           >  1:6209431/1‑150 (MQ=255)
                                                                                           tCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTa                                                           >  1:395558/1‑150 (MQ=255)
                                                                                           gcccAGTCGCGGGCAGGCCCGCTCCCACAGTACTCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGCGCGGGCGTGCCCGCGAAAAATCGCACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTa                                                           <  1:410027/149‑1 (MQ=255)
                                                                                            cccATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAg                                                          >  1:1103657/1‑150 (MQ=255)
                                                                                             ccATGCGCGGGCTTGCCCGCTCCCCCAGTACTCGCACCGCACTTACGGCCTGGGCAGGACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGc                                                         <  2:5187018/150‑1 (MQ=255)
                                                                                              cATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGcc                                                        >  2:380230/1‑150 (MQ=255)
                                                                                               aTTCGCGGGCATGCCCGCTCCCCCAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGccc                                                       <  1:5305150/150‑1 (MQ=255)
                                                                                               aTTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGccc                                                       >  1:3320280/1‑150 (MQ=255)
                                                                                               aTTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGccc                                                       <  1:1557928/150‑1 (MQ=255)
                                                                                               aTTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACGTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCGAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACGTCCAGCGCACGTACAGTAGccc                                                       >  2:1074/1‑150 (MQ=255)
                                                                                                 tCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCgg                                                     >  1:3349351/1‑150 (MQ=255)
                                                                                                  cgcgGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGc                                                    <  1:4570005/150‑1 (MQ=255)
                                                                                                    cgGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCaa                                                  >  2:5120098/1‑150 (MQ=255)
                                                                                                    cgGGCACGCCCGCTCCCACAGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGCGGGCATGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCaa                                                  >  1:6739310/1‑150 (MQ=17)
                                                                                                      ggcgcgCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCACTg                                                >  1:822025/6‑150 (MQ=255)
                                                                                                      ggCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGGGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATg                                                <  2:2276816/150‑1 (MQ=255)
                                                                                                       gCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGa                                               <  1:3014829/150‑1 (MQ=255)
                                                                                                         acgCCCGCTCCCACAGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCACTGAAc                                             >  1:6831371/3‑150 (MQ=17)
                                                                                                           gCCCGCTCCCACAGTAATCGCACCGCGCTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAg                                           <  2:4121420/150‑1 (MQ=255)
                                                                                                           gCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAg                                           <  2:2238470/150‑1 (MQ=255)
                                                                                                           gCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAg                                           >  1:1992931/1‑150 (MQ=255)
                                                                                                           gCCCGCTCCCACAGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGCGGGCATGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAg                                           >  2:1134971/1‑150 (MQ=17)
                                                                                                             ccGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGcc                                         >  2:2589441/1‑150 (MQ=255)
                                                                                                             ccGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGcc                                         >  1:1916275/1‑150 (MQ=255)
                                                                                                             ccGCTCCCACAGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGcc                                         >  1:5017074/1‑150 (MQ=18)
                                                                                                              cGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGccc                                        >  1:255416/1‑150 (MQ=255)
                                                                                                               gCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGcccc                                       >  2:5628258/1‑150 (MQ=255)
                                                                                                                 tCCCACAGTAATCGCACCGCACTTAAGGCCTGAGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCAAGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAg                                     >  2:6209740/1‑150 (MQ=255)
                                                                                                                  cccACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAgg                                    <  1:5046292/150‑1 (MQ=255)
                                                                                                                   ccACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATc                                                                                                                     >  2:5018663/1‑68 (MQ=255)
                                                                                                                   ccACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATc                                                                                                                     <  1:5018780/68‑1 (MQ=255)
                                                                                                                    cacaGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCt                                  >  1:1632792/1‑150 (MQ=255)
                                                                                                                    cacaGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCt                                  >  1:1289128/1‑150 (MQ=255)
                                                                                                                     acaGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTc                                 >  1:844017/1‑150 (MQ=255)
                                                                                                                     acaGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTc                                 >  2:2066539/1‑150 (MQ=255)
                                                                                                                       aGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCAc                               <  1:1877756/150‑1 (MQ=255)
                                                                                                                       aGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCAc                               <  2:1367917/150‑1 (MQ=255)
                                                                                                                         tAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCa                             >  2:846855/1‑150 (MQ=255)
                                                                                                                          aaTCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAg                            <  1:3161127/150‑1 (MQ=255)
                                                                                                                            tCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCa                          >  1:5492211/1‑150 (MQ=255)
                                                                                                                              gcaccgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACt                        >  2:2881535/1‑150 (MQ=255)
                                                                                                                               caccgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACtt                       >  1:3827084/1‑150 (MQ=255)
                                                                                                                                accgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCACGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGAGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTc                      >  2:70794/1‑150 (MQ=255)
                                                                                                                                 ccgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCa                     <  1:5175940/150‑1 (MQ=255)
                                                                                                                                  cgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCaa                    >  2:1101100/1‑150 (MQ=255)
                                                                                                                                       ttAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAagccag               >  2:4446229/1‑150 (MQ=255)
                                                                                                                                       ttAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAagccag               >  1:563344/1‑150 (MQ=255)
                                                                                                                                         aaGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAagccagcc             >  2:2575559/1‑150 (MQ=255)
                                                                                                                                            ggcTGTGCAGTACCTGGGGGTGGGGGGGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGaa          <  1:42679/148‑1 (MQ=255)
                                                                                                                                            aaaagtccAGTACCTGAGGGAGCGGGCGTGCACGCGAAAAATCTAACTGCCTTAGCCCGCCATGCGCCGCTCATGCAGGAACTACCAGCGCACGAACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACATCAAGCCAGCCGaa          >  2:5807335/8‑150 (MQ=255)
                                                                                                                                                  gCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATAccgc    >  1:472193/1‑150 (MQ=255)
                                                                                                                                                  gCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATAccgc    <  2:3422193/150‑1 (MQ=255)
                                                                                                                                                  gCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATAccgc    >  2:4978522/1‑150 (MQ=255)
                                                                                                                                                   cAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATAccgcc   >  1:5608765/1‑150 (MQ=255)
                                                                                                                                                    aGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATAccgccc  >  2:6685525/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTCAACCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCCGAATACCGCCC  >  NC_002947/1842966‑1843263

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