Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,377,810 C→G 9.5% intergenic (+268/+24) PP_2970 → / ← PP_2971 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,377,8100CG 9.5% 585.7 / 15.8 201intergenic (+268/+24)PP_2970/PP_2971hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base C (108/71);  new base G (3/19);  total (111/90)
Fisher's exact test for biased strand distribution p-value = 4.51e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 209 total aligned reads displayed.

GTTGTTTATGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACTTGCCGCCCAT  >  NC_002947/3377661‑3377959
                                                                                                                                                      |                                                                                                                                                     
gttgttTATGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAgc                                                                                                                                                       >  1:3452047/1‑150 (MQ=255)
gttgttTATGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAgc                                                                                                                                                       <  1:317596/150‑1 (MQ=255)
    tttATGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAgcggcg                                                                                                                                                   >  1:3991882/1‑150 (MQ=255)
     ttATGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGa                                                                                                                                                  >  2:5521070/1‑150 (MQ=255)
        tGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTc                                                                                                                                               <  2:3682131/150‑1 (MQ=255)
        tGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTc                                                                                                                                               <  2:1164283/150‑1 (MQ=255)
        tGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTc                                                                                                                                               <  2:3992975/150‑1 (MQ=255)
        tGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTc                                                                                                                                               >  2:3981680/1‑150 (MQ=255)
         gggCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCg                                                                                                                                              <  1:5538508/150‑1 (MQ=255)
         gggCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCg                                                                                                                                              >  1:1759196/1‑150 (MQ=255)
         gggCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCg                                                                                                                                              <  2:359971/150‑1 (MQ=255)
           gCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGc                                                                                                                                            <  2:1194054/150‑1 (MQ=255)
           gCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGc                                                                                                                                            >  2:2978827/1‑150 (MQ=255)
           gCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCGCTGCAGAGGGCGATTCgga                                                                                                                                            <  2:5041830/150‑2 (MQ=255)
            cACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGcc                                                                                                                                           <  2:3100610/150‑1 (MQ=255)
             aCCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGccc                                                                                                                                          <  2:2640001/150‑1 (MQ=255)
               cTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCAc                                                                                                                                        >  1:4711767/1‑150 (MQ=255)
                tCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACt                                                                                                                                       >  1:1973487/1‑150 (MQ=255)
                tCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACt                                                                                                                                       >  2:340370/1‑150 (MQ=255)
                tCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACt                                                                                                                                       >  1:3002950/1‑150 (MQ=255)
                  gAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTcc                                                                                                                                     >  2:3120768/1‑150 (MQ=255)
                      aaCGCGAGGTTTTCAGTGAATCGGACCGCCAAATGATGGCCATTATCGATTTCGGTTCTATTTCTGCGGCTAATCACCGGTTTCTGGGCGGTGTCGCCCCTTTTTTCTTCCGAATCGCCCTCTGCTGGGGCGATTCGGCCCACTCCGGGc                                                                                                                                 <  2:3573438/150‑1 (MQ=255)
                          cgAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCa                                                                                                                             >  1:1879081/1‑150 (MQ=255)
                             gtttttCTGTGAATCCGACCGCAAAATGAGGGCCATGATCTAATTCGGTTCGATTTGTGCGGCGAATCACGGCTTTCTGGGGGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCTATTCGGCCCACTCCGGGCGTCAGaa                                                                                                                          <  2:4279086/148‑1 (MQ=255)
                                tttCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGc                                                                                                                       >  2:4404820/1‑150 (MQ=255)
                                   ctgTGAACCCGACCGCCATTTGTTGGCCCGGATCTACTTGGGTACGATTTCGGCGCCGAACCCCCGCTTTGGGGGGGGGGGCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCTATTCGGCCCACTCCGGGCGTCAGAAGGCCtt                                                                                                                    <  2:3306406/148‑1 (MQ=255)
                                   cAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCtt                                                                                                                    >  1:3317720/1‑150 (MQ=255)
                                    aGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTg                                                                                                                   >  1:4393521/1‑150 (MQ=255)
                                    aGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTg                                                                                                                   >  1:1281890/1‑150 (MQ=255)
                                     gTGATTCCGACGGCCAATGGTGGGCCATGAGCGAATGGGGTGCGTTTTCGGGGGCGACTCCCCGGTTGCGGGGGGGGGGCGCCCGGTTTTGCTTCCGAATCGCCGTCTGCCGGGGGGCTGGGGCCCACTCCGGGCGTCAGAAGGCCTTGa                                                                                                                  <  2:2834532/150‑1 (MQ=255)
                                     gTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGa                                                                                                                  >  1:4609064/1‑150 (MQ=255)
                                     gTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGa                                                                                                                  >  1:2651430/1‑150 (MQ=255)
                                      tGAATCCGACCGCCAAATGATGGCCATGATCGAATTGGGTTCGTTTTCTGCGGCGAATCACCGCTTTCGGGGCGGTGTCGCCCCTTTTTTCTTCCGAATCGCCCTCTGCAGGGGCGTTTCGGCCCACTCCGGGCGTCAGAAGGCCTTGAt                                                                                                                 <  1:5043393/150‑1 (MQ=255)
                                       gAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGAtt                                                                                                                <  2:4644493/150‑1 (MQ=255)
                                        aaTCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTc                                                                                                               >  2:174756/1‑150 (MQ=255)
                                         aTCCCACCGCAAATTGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCTTTTTTCTTCCTAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTcc                                                                                                              <  2:810984/150‑1 (MQ=255)
                                             gACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGcc                                                                                                          >  1:5486165/1‑150 (MQ=255)
                                             gACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGcc                                                                                                          >  1:4411788/1‑150 (MQ=255)
                                                 gCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCt                                                                                                      >  2:181237/1‑150 (MQ=255)
                                                  ccAAATGTGGGCCATTTTCGAATTCGGTTGGATTTCTGCGGGGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTGGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTg                                                                                                     <  2:2244847/150‑1 (MQ=255)
                                                  ccAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTg                                                                                                     >  2:5390083/1‑150 (MQ=255)
                                                     tttttgggCCTTTATCGATTTCGGTTCTTTTTCTGCGGCGAAACACCTCTTGCTGGGGGGTGGCGCCCCTTTTTTCTTCCCAACCGCCCTCTGCCGCGGCTATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGaa                                                                                                  <  2:2385722/144‑1 (MQ=255)
                                                          tgGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTaca                                                                                             >  2:4079558/1‑150 (MQ=255)
                                                             ccATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCa                                                                                          >  2:5333912/1‑150 (MQ=255)
                                                                           ggTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTg                                                                            >  2:3768171/1‑150 (MQ=255)
                                                                            gTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGGGTTGTCGCCCCTTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGt                                                                           <  2:4445632/150‑1 (MQ=255)
                                                                             ttCTATTTCTGGGGCGAATCACCGCTTTGTGGGGGGTGTCGCCCGGTTTTTCTTCCTAATCGCCCTCTGCTGCGGCTTTTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTa                                                                          <  1:3261347/150‑1 (MQ=255)
                                                                             ttCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTa                                                                          >  2:1469737/1‑150 (MQ=255)
                                                                             ttCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTa                                                                          >  2:1775699/1‑150 (MQ=255)
                                                                             ttCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTa                                                                          >  2:5183730/1‑150 (MQ=255)
                                                                             ttCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTa                                                                          >  2:2050658/1‑150 (MQ=255)
                                                                             tgcGATTTCTGCGGCGAACCCCCGATTTCGGGGCGGTGGCGCCCCGTTTTTCTCCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTa                                                                          <  2:3849857/148‑1 (MQ=255)
                                                                              tCGTTTTCTGCGGCTAATCACCGCTTTCTGGGCGGGGTCGCCCCTTTTTCCGTCCGAATCGCCCTCTGCAGCGGCTATTCGCCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAg                                                                         <  2:2482017/150‑1 (MQ=255)
                                                                                gttttCTGCGGCTAATCACCGATTTCTGGGGGGTGTCGCCCCTTTTTACTTCCGAATCGCCCTCTGCCGCGGCGTTCCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCa                                                                       <  2:2573451/148‑1 (MQ=255)
                                                                                 aTTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAg                                                                      >  1:4912498/1‑150 (MQ=255)
                                                                                   ttCTGCGGCGAATCACCGCTTTCGGGGCGGGGTCGCCCCTTTTTTCTTCCGAATCGCCCTCGGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGgc                                                                    <  2:731335/150‑1 (MQ=255)
                                                                                       gcggcgAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGgccgcc                                                                >  1:804581/1‑150 (MQ=255)
                                                                                         ggctAATCCCGGCTTTGGGGGGGGTGGCGCCCCTTTTTTCTTCCTACGCGCCGTGCGCCGGGGCGCTGCGCCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCcat                                                              <  1:1165363/150‑1 (MQ=255)
                                                                                           ctaaTCCCCGTTTTGTGGGGGGTGTCGCCCCTTTTTCCTACCGAAACGCCCTCTGCCGGGGCTATGCGCCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCcatca                                                            <  1:4377968/148‑1 (MQ=255)
                                                                                            gAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCcatcat                                                           >  2:1960999/1‑150 (MQ=255)
                                                                                             aaTCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCTATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCatcatc                                                          <  2:4066291/150‑1 (MQ=255)
                                                                                               tCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgt                                                        >  2:2836314/1‑150 (MQ=255)
                                                                                                cACGGCTTTGTGGGGGGTGTCGCCCCTTTTTTCTTCCGAATCGCCCTCTGCAGCGGGGATTGGCCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgtc                                                       <  1:5469364/150‑1 (MQ=255)
                                                                                                 aCCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgtca                                                      >  2:1869590/1‑150 (MQ=255)
                                                                                                  ccGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgtcat                                                     >  1:3746981/1‑150 (MQ=255)
                                                                                                   cGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCGGCAGCGGCTATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgtcatc                                                    <  1:5118389/150‑1 (MQ=255)
                                                                                                    gCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgtcatcg                                                   <  2:3720996/150‑1 (MQ=255)
                                                                                                     cTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGc                                                  >  1:3703781/1‑150 (MQ=255)
                                                                                                        tcgggggggTGTCGCCCCTTTTTTCTTCCGAATCGCCCTGTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCaaa                                               <  2:4711767/147‑1 (MQ=255)
                                                                                                          tGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGt                                             >  1:2462403/1‑150 (MQ=255)
                                                                                                            gggggTGTCGCCCTGTTTTTCTTCCGAAGCGCCCTCTGCCGGGGCTTTGGGCCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTc                                           <  1:3531570/147‑1 (MQ=255)
                                                                                                             gCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTcg                                          >  2:4589467/1‑150 (MQ=255)
                                                                                                              gggTGTCGCCCCTTTTTTCTTCCGAAGCGCCCTGTGCGGGGGGTTTCCGCCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTcgc                                         <  1:1157595/149‑1 (MQ=255)
                                                                                                              cGGTGTCGCCCCTTTTTTCTTCCGACTCGCCCTCTGCCGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTcgc                                         <  2:3002950/150‑1 (MQ=255)
                                                                                                                gtgtCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTcgcgc                                       <  2:3223302/150‑1 (MQ=255)
                                                                                                                  gtCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCgcgcgc                                     <  2:3452047/150‑1 (MQ=255)
                                                                                                                    cGCCCCTTTTTTCTCCCGATTCGCCCTCTGCTGGGGCGATTCGCCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGccc                                   <  1:5419802/150‑1 (MQ=255)
                                                                                                                      ccccTTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCt                                 <  2:1920453/150‑1 (MQ=255)
                                                                                                                      ccccGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCCTCATCGTCATCGCACAGTTCGCGCGccccg                                 >  2:5437632/1‑149 (MQ=255)
                                                                                                                      ccccGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCt                                 >  1:1026164/1‑150 (MQ=255)
                                                                                                                       cccGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTg                                >  2:554899/1‑150 (MQ=255)
                                                                                                                       cccGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTg                                >  1:216021/1‑150 (MQ=255)
                                                                                                                       cccGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTAGCGCGCCCCTg                                >  2:5232379/1‑150 (MQ=255)
                                                                                                                         cGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTc                              <  2:1274100/150‑1 (MQ=255)
                                                                                                                          gTTTTTCTTCCGAATCGCCCTCTGCCGGGGCTTTTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTcc                             <  1:1482609/150‑1 (MQ=255)
                                                                                                                          gTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTcc                             >  2:2600621/1‑150 (MQ=255)
                                                                                                                             tttCTTCCGAATCGCCCTGTTCAGCGGCGTTTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGAt                          <  2:4652475/150‑1 (MQ=255)
                                                                                                                             tttCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGAt                          <  2:1606176/150‑1 (MQ=255)
                                                                                                                                 ttcCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTa                      >  1:1539364/1‑150 (MQ=255)
                                                                                                                                 ggccGAATCGCCGCTGCAGAGGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTa                      >  1:165464/3‑150 (MQ=255)
                                                                                                                                   ccGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACg                    >  1:2865311/1‑150 (MQ=255)
                                                                                                                                    ctaaTCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACgg                   <  1:856027/148‑1 (MQ=255)
                                                                                                                                    cGAATCGCCGCTGCAGAGGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAgg                                                                     <  2:4324263/100‑1 (MQ=38)
                                                                                                                                    cGAATCGCCGCTGCAGAGGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAgg                                                                     >  1:4324263/1‑100 (MQ=38)
                                                                                                                                      aaTCGCCCTCTGCCGCGGCGATTCGCCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGAt                 <  2:505663/150‑1 (MQ=255)
                                                                                                                                      aaTCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGAt                 >  2:4215041/1‑150 (MQ=255)
                                                                                                                                            ccTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACtt           >  1:419554/1‑150 (MQ=255)
                                                                                                                                               ctGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTATAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACTTgcc        >  1:5370728/1‑150 (MQ=255)
                                                                                                                                                   aGCGGCTATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACTTGCCGccc    <  1:340370/150‑1 (MQ=255)
                                                                                                                                                    gcggcgATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACTTGCCGCCCa   >  2:3515660/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GTTGTTTATGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACTTGCCGCCCAT  >  NC_002947/3377661‑3377959

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: