Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,020,779 A→G 6.9% intergenic (‑45/+189) PP_5679 ← / ← PP_4421 hypothetical protein/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,020,7790AG 6.9% 680.8 / 26.0 217intergenic (‑45/+189)PP_5679/PP_4421hypothetical protein/aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (96/106);  new base G (7/8);  total (103/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

Only 100 of 220 total aligned reads displayed.

ACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAA  >  NC_002947/5020630‑5020928
                                                                                                                                                     |                                                                                                                                                      
   cACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCg                                                                                                                                                     >  2:144983/1‑150 (MQ=255)
    aCCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCgg                                                                                                                                                    >  2:3218606/1‑150 (MQ=255)
           gCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCgg                                                                                                                                             >  2:4776966/1‑150 (MQ=255)
           gCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCgg                                                                                                                                             >  1:4960647/1‑150 (MQ=255)
            ccGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGgc                                                                                                                                            >  1:4519183/1‑150 (MQ=255)
              gCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGgcgc                                                                                                                                          >  2:1653718/1‑150 (MQ=255)
                ttGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCGCGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCgg                                                                                                                                        <  2:3350020/150‑1 (MQ=255)
                ttGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCgg                                                                                                                                        <  2:2144996/150‑1 (MQ=255)
                 tGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCGCCGGAGACCGGTGCggg                                                                                                                                       <  2:4895824/150‑1 (MQ=255)
                    agagTGCCTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCGCCGGAGACCGGTGCGGGCtt                                                                                                                                    <  1:3890686/150‑1 (MQ=255)
                       gTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCACACCGGTCTCCGGCGCGGGCTTTTa                                                                                                                                 <  1:2253442/150‑1 (MQ=255)
                           ctgctgGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCGCCGGAGACCGGTGCGGGCttattcatc                                                                                                                             <  1:3457443/150‑8 (MQ=39)
                               tgGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCgg                                                                                                                         >  2:4070526/1‑150 (MQ=255)
                                 ggCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGc                                                                                                                       >  1:5483664/1‑150 (MQ=255)
                                 ggCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGc                                                                                                                       <  1:1206042/150‑1 (MQ=255)
                                 ggCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGc                                                                                                                       >  2:4543864/1‑150 (MQ=255)
                                  gCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGcc                                                                                                                      >  2:2109538/1‑150 (MQ=255)
                                    aCGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCtt                                                                                                                    <  2:5399199/150‑1 (MQ=255)
                                      gAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTgg                                                                                                                  >  2:2892254/1‑150 (MQ=255)
                                      gAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTgg                                                                                                                  >  2:1592596/1‑150 (MQ=255)
                                         gCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACa                                                                                                               <  2:4495704/150‑1 (MQ=255)
                                         gCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACa                                                                                                               <  2:3795694/150‑1 (MQ=255)
                                             ctcGGCAGCCCTGGCCGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCGTAGCGCTCTCGTTTGAGTCCACGATTAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGTCAGTAt                                                                                                           <  2:563750/150‑1 (MQ=255)
                                             ctcGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCGTCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTAt                                                                                                           <  2:1419728/150‑1 (MQ=255)
                                              tcGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATc                                                                                                          <  2:289836/150‑1 (MQ=255)
                                                  cAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGAc                                                                                                      <  1:261031/150‑1 (MQ=255)
                                                   aGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGAcc                                                                                                     <  1:1972352/150‑1 (MQ=255)
                                                     cccTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCg                                                                                                   >  2:3660251/1‑150 (MQ=255)
                                                      ccTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGa                                                                                                  <  1:2146253/150‑1 (MQ=255)
                                                      ccTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCACGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGa                                                                                                  >  2:2913568/1‑150 (MQ=255)
                                                        tGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGaaa                                                                                                >  2:3829771/1‑150 (MQ=255)
                                                        tGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGaaa                                                                                                <  2:2254933/150‑1 (MQ=255)
                                                        tGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGaaa                                                                                                <  2:3698623/150‑1 (MQ=255)
                                                        tGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGaaa                                                                                                <  2:4256518/150‑1 (MQ=255)
                                                            cGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTAc                                                                                            <  1:5415130/150‑1 (MQ=255)
                                                                  ggTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCtt                                                                                      <  1:1788342/150‑1 (MQ=255)
                                                                   gTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTtg                                                                                     >  1:2636842/1‑150 (MQ=255)
                                                                   gTAGCCACTGCTGCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTCTCTGACCCGACAGTACGTCCTtg                                                                                     >  2:1023442/1‑150 (MQ=255)
                                                                    tAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTtgt                                                                                    <  1:3689025/150‑1 (MQ=255)
                                                                    tAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGACCTtgt                                                                                    >  1:2839746/1‑150 (MQ=255)
                                                                       ccACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGaa                                                                                 <  1:4627511/150‑1 (MQ=255)
                                                                       ccACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGaa                                                                                 <  1:5072491/150‑1 (MQ=255)
                                                                        cACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTGTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAAc                                                                                >  1:4904491/1‑150 (MQ=255)
                                                                        cACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAAc                                                                                <  1:398348/150‑1 (MQ=255)
                                                                            gCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACtgat                                                                            <  1:1415320/150‑1 (MQ=255)
                                                                            gCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACtgat                                                                            <  1:4176334/150‑1 (MQ=255)
                                                                             cttcttTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACtgatg                                                                           >  1:1240730/1‑150 (MQ=255)
                                                                              ttcttTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACtgatga                                                                          >  2:581599/1‑150 (MQ=255)
                                                                              ttcttTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACtgatga                                                                          >  2:1177562/1‑150 (MQ=255)
                                                                                  ttGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGc                                                                      >  1:1034699/1‑150 (MQ=255)
                                                                                    gATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCttt                                                                                                                                   <  2:266944/87‑1 (MQ=255)
                                                                                    gATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCttt                                                                                                                                   >  1:266944/1‑87 (MQ=255)
                                                                                    gATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGt                                                                    <  1:1005223/150‑1 (MQ=255)
                                                                                     aTCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTc                                                                   <  1:3646518/150‑1 (MQ=255)
                                                                                     aTCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTc                                                                   <  2:4411976/150‑1 (MQ=255)
                                                                                      tCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCt                                                                  >  2:2614035/1‑150 (MQ=255)
                                                                                       cGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCtt                                                                 <  2:881114/150‑1 (MQ=255)
                                                                                       cGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCtt                                                                 <  1:856165/150‑1 (MQ=255)
                                                                                       cGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCtt                                                                 <  1:3555582/150‑1 (MQ=255)
                                                                                        ggAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTg                                                                <  1:3432360/150‑1 (MQ=255)
                                                                                         gAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTgg                                                               >  2:3901866/1‑150 (MQ=255)
                                                                                          aGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGc                                                              <  1:4810304/150‑1 (MQ=255)
                                                                                           gCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCt                                                             >  1:3347589/1‑150 (MQ=255)
                                                                                             cAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCttt                                                           <  2:4321298/150‑1 (MQ=255)
                                                                                               gCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGa                                                         <  1:3187343/150‑1 (MQ=255)
                                                                                                 tCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGAtt                                                       >  1:688171/1‑150 (MQ=255)
                                                                                                     tCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGa                                                   <  1:3327140/150‑1 (MQ=255)
                                                                                                       aTCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATa                                                 <  2:4680016/150‑1 (MQ=255)
                                                                                                            aGCGCTCTCGATAGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCAATTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACt                                            <  2:5132263/150‑1 (MQ=255)
                                                                                                             gcgcTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTg                                                                >  2:2806055/1‑129 (MQ=255)
                                                                                                             gcgcTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTg                                                                <  1:2806055/129‑1 (MQ=255)
                                                                                                                 tctcGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTGTGATTGTGATAGCACTACaaa                                      >  2:5399457/1‑150 (MQ=255)
                                                                                                                  ctcGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaa                                     <  1:3218606/150‑1 (MQ=255)
                                                                                                                   tcGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaa                                    >  2:660512/1‑150 (MQ=255)
                                                                                                                   tcGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaa                                    >  1:3142432/1‑150 (MQ=255)
                                                                                                                    cGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaa                                   <  1:4289587/150‑1 (MQ=255)
                                                                                                                     gTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaaa                                  <  2:1037895/150‑1 (MQ=255)
                                                                                                                         gAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaaaaaaa                              <  1:4605822/150‑1 (MQ=255)
                                                                                                                          aGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAg                             >  1:5357357/1‑150 (MQ=255)
                                                                                                                             ccACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACATTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGccc                          <  2:4960899/150‑1 (MQ=255)
                                                                                                                             ccACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaaaaaaaagcc                          >  1:1451726/1‑146 (MQ=255)
                                                                                                                              cACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAATAAAATAGCCCg                         >  2:937559/1‑150 (MQ=255)
                                                                                                                                 gatgagatgaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATTGCACTACAAAAAAAAAAGCCCGCCTa                     <  1:3524539/150‑1 (MQ=255)
                                                                                                                                 gatgagatgaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCt                      >  2:3504466/1‑150 (MQ=255)
                                                                                                                                   tgagatgaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTaaa                   <  1:2109538/150‑1 (MQ=255)
                                                                                                                                       atgaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGc                >  1:140739/1‑150 (MQ=255)
                                                                                                                                         gaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCgg              >  1:5304849/1‑150 (MQ=255)
                                                                                                                                          taaatgCCCGCGCCGGAGACCGGTGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAATCCCGCCTAAAGGCggg             >  1:1563952/6‑150 (MQ=255)
                                                                                                                                          aaTAAGCCCGCACCGGTCTCCGGCGCGGTCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCTAAAGTACGTCCTTGTGAACTGATGACAGATTATTGTCTTTTATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCggg             <  1:1104696/150‑1 (MQ=255)
                                                                                                                                          aaTAAGCCCGCACCGGTATCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCggg             <  2:412116/150‑1 (MQ=255)
                                                                                                                                             aaGCCCGCGCCGGAGACCGGTGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACa                                        >  2:4651814/1‑120 (MQ=255)
                                                                                                                                              aGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACTGTATCTGCCCCGAAAGTACGTCCTGGTGAAGTGATGACAGCGTCTTGGCTTGGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTTTAGGCGGGCttt         <  2:3418513/150‑2 (MQ=255)
                                                                                                                                               gCCCGCGCCGGAGACCGGTGCGGGCTTTTACTTTCCGGGCCTTGGACCGTAGCTGACCCGAAAGGACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCGAACGGCGGGCGTa         >  2:3560139/1‑150 (MQ=255)
                                                                                                                                               gCCCGCGCCGGAGACCGGTGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGa                                                   >  2:1628986/1‑108 (MQ=39)
                                                                                                                                               gCCCGCGCCGGAGACCGGTGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGa                                                   <  1:1628986/108‑1 (MQ=39)
                                                                                                                                               gCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAACGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTCCAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTaa        >  2:1929835/1‑150 (MQ=255)
                                                                                                                                                 ccGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTTTAGGCGGGCtttttt      <  1:4978698/150‑5 (MQ=255)
                                                                                                                                                 ccGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTaaaa      >  2:1319243/1‑150 (MQ=255)
                                                                                                                                                 ccGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTaaaa      <  1:3829771/150‑1 (MQ=255)
                                                                                                                                                     aCCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGaa  <  1:4550771/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
ACTCACCCGTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAA  >  NC_002947/5020630‑5020928

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: