Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,449,251 A→C 14.5% S18R (AGC→CGC)  PP_5556 → bacteriophage lysis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,449,2510AC14.5% 360.5 / 13.2 197S18R (AGC→CGC) PP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base A (92/67);  new base C (15/18);  total (107/90)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

AAGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTC  >  NC_002947/3449102‑3449400
                                                                                                                                                     |                                                                                                                                                     
aaGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCa                                                                                                                                                       <  1:781500/150‑1 (MQ=255)
  gggCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAgc                                                                                                                                                     <  1:7610163/150‑1 (MQ=255)
    gCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAgcgc                                                                                                                                                   >  1:6609327/1‑150 (MQ=255)
        cATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggc                                                                                                                                               <  1:1320717/150‑1 (MQ=255)
          tCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcgg                                                                                                                                             <  1:4203664/150‑1 (MQ=255)
           ccTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggc                                                                                                                                            <  1:684967/150‑1 (MQ=255)
           ccTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggc                                                                                                                                            <  1:6318076/150‑1 (MQ=255)
            cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGgcggcg                                                                                                                                           >  2:8466181/1‑150 (MQ=255)
             tGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTTGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGgcgc                                                                                                                                          >  2:4956847/1‑150 (MQ=255)
             tGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGgcgc                                                                                                                                          <  2:7065996/150‑1 (MQ=255)
             tGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGgcgc                                                                                                                                          >  2:4791758/1‑150 (MQ=255)
             tGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGgcgc                                                                                                                                          >  2:6615362/1‑150 (MQ=255)
             gggCGAGCACCAAACGTATCAAGGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCGGCTGCTGGCCGCCAGCGCTGGCGgcgc                                                                                                                                          <  2:4616868/149‑1 (MQ=255)
                 gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGccccg                                                                                                                                      >  1:5660140/1‑148 (MQ=255)
                    cACCAAAGGTTTCAATGGTGGTAGAAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCa                                                                                                                                   <  2:4437636/150‑1 (MQ=255)
                    cACCAAACGTATCACTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCa                                                                                                                                   <  2:4775066/150‑1 (MQ=255)
                     aCCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAg                                                                                                                                  >  2:2507442/1‑150 (MQ=255)
                      ccAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGc                                                                                                                                 >  1:3170580/1‑150 (MQ=255)
                       cAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGTTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCCCCCGGgccgcc                                                                                                                                >  2:5691153/1‑150 (MQ=255)
                         accGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCACTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTg                                                                                                                              <  2:1509761/148‑1 (MQ=255)
                            gTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGgt                                                                                                                           <  2:8163771/150‑1 (MQ=255)
                                 aaTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGc                                                                                                                      >  2:3969147/1‑150 (MQ=255)
                                     gtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGc                                                                                                                  >  2:7491694/1‑150 (MQ=255)
                                     gtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGc                                                                                                                  >  2:3226860/1‑150 (MQ=255)
                                         tAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAAc                                                                                                              >  2:5503720/1‑150 (MQ=255)
                                           gCAACGGTCTGGAGGAGCGCCAGGCCCTCTCCCCGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTGGGGCCCCCCGCCGGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGc                                                                                                            <  1:7955797/150‑1 (MQ=255)
                                            cAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGGGTGGCAGGCCAACGcc                                                                                                           >  2:5633654/1‑150 (MQ=255)
                                             aaCGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAAccccg                                                                                                          >  1:4988856/1‑149 (MQ=255)
                                             aaCGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCaagatcg                                                                                                          >  1:3102622/1‑145 (MQ=255)
                                             aaCGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGCCGGCCACCCCCt                                                                                                          >  1:3986408/1‑150 (MQ=255)
                                                  tCTGGTGGATCGCAAGGCCCTCTACCCGCGAGCGCTCGTGGTGCGGCCGGGCGCGCTTGGGCCACCCGCCTGGTCGCGGTGGCGCGGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAg                                                                                                     <  1:2825951/150‑1 (MQ=255)
                                                  tCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAg                                                                                                     >  2:3056920/1‑150 (MQ=255)
                                                  tCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAg                                                                                                     >  2:7315123/1‑150 (MQ=255)
                                                     ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgcg                                                                                                  >  1:2534127/1‑150 (MQ=255)
                                                        ggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCCGGCCAACGCCTCGAAGGCGCAg                                                                                               >  1:4112368/1‑150 (MQ=255)
                                                         gaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGc                                                                                              <  2:6609327/150‑1 (MQ=255)
                                                         gaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGc                                                                                              >  1:8544678/1‑150 (MQ=255)
                                                         gaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCCCCTGGGTGGGGCAGGCCCACGCCAAAAAGGCGCAGc                                                                                              >  1:3939676/1‑150 (MQ=255)
                                                          ttCGCAAGGCCATCTACAAGCGAGCGCTGGTGGTGCTGCTGTGCGCGGTTGGGCCACCCGCCTGGTCGCGGTGGCGGTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCt                                                                                             <  1:7335317/149‑1 (MQ=255)
                                                             gCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGc                                                                                          >  2:6863996/1‑150 (MQ=255)
                                                              cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCCGGCAGCCTGGGGGTGGCAGGCCAACCCCTAAAACGCGCAGCTCGCt                                                                                         >  1:1695175/1‑150 (MQ=255)
                                                               aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGAAAAGGCGCAGCTCGcgg                                                                                        >  1:8224139/1‑148 (MQ=255)
                                                                 ggCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCGGTGGCGCTGCTGCTGGCCGCAAGCGCTGGCGGCTCCCGGGCAGCCTGTGTGTGGCAGTCCAACGCCTATAAGGCGCAGCTCGCTGAt                                                                                      <  2:9100207/150‑1 (MQ=255)
                                                                  gCCCTCTACAAGGGCGGGATCGAGGTGCTGCAGTGAGCGCTTGGGCAACCGGCCTGGGCGCGGTGGCGCTGCTGGTGGCCGCCAGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATc                                                                                     <  1:7444824/150‑1 (MQ=255)
                                                                      tctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACCCCGATAAGGCGCAGCTCGCGGATCAGGc                                                                                 >  1:8983836/1‑150 (MQ=255)
                                                                       ctACAAGCCAGCGCTCGAGGTGCTCAAGTCAGCGCTTGGGCCACCCGCCTGGTCCCCGTGGCGCTGCTGCTGGCCCCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGcc                                                                                <  2:1308357/150‑1 (MQ=255)
                                                                        tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGGGGCAGGCCAACCCCTATAAGGCGCAGCGCGCGGATCAGGCCg                                                                               >  1:4468427/1‑150 (MQ=255)
                                                                         aCAAGGGAGCGCTCGTGGTGCTGCTGTGCGCGCTTGGGCCACCCGCCTGGGCGCCGTGGCGGTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGa                                                                              <  1:4225663/150‑1 (MQ=255)
                                                                                cgcgcgcgcGGTGCGGCTGGGAGCGCGTGGGCCCCCCGCCTGGGCGCCGGGGCGGTGCGGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTc                                                                       <  1:4938171/149‑1 (MQ=255)
                                                                                agcgCTCGCGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCCGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTc                                                                       <  1:9166771/150‑1 (MQ=255)
                                                                                 gcgcTCGTGGTGCTGCCGTGAGCGCTGGGGCCACCCGCCTGGTCGCGGTGGCGCTGCTGGTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCg                                                                      <  1:6728582/150‑1 (MQ=255)
                                                                                    cTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCCGCCTGGGGGTGGCAGGCCCCCGCCTATAAGGCGCCGCTCGCTGATCAGGCCGAGGTCTTCGGCa                                                                   >  1:2015426/1‑150 (MQ=255)
                                                                                     tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGTCGTCGGCaa                                                                  >  2:8705347/1‑150 (MQ=255)
                                                                                       gcggTGCTTCAGTGAGCGGTTGGGCCCCCGGCCTGGGGGCGGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGc                                                                <  1:7528768/148‑1 (MQ=255)
                                                                                       gAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTACTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGc                                                                >  2:1709901/1‑150 (MQ=255)
                                                                                       gAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTACTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGc                                                                >  2:1710092/1‑150 (MQ=255)
                                                                                         ggTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCcagcag                                                              >  1:6556949/1‑150 (MQ=255)
                                                                                             ctgcCGTGCGCGGTGGGGCCACCGGCGTGGTCGCCGTGGGGGTGGTGCTGGCCGCCCGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTgg                                                          <  1:1084043/150‑1 (MQ=255)
                                                                                              tgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCTACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGc                                                         >  2:2974352/1‑150 (MQ=255)
                                                                                              tgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCa                                                                                    >  1:1338204/1‑123 (MQ=255)
                                                                                              tgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCa                                                                                    <  2:1338204/123‑1 (MQ=255)
                                                                                              tgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCACGCAGTTGGc                                                         >  2:5719406/1‑150 (MQ=255)
                                                                                               gcAGTGAGCGCTTGGGCCACCCGCCTGGTGGCCGTGGCGCTGCTGCTGGCCGCCTGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCg                                                        <  2:3127531/150‑1 (MQ=255)
                                                                                               gcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCCCGCCTATAAGGCGCAGCGCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCg                                                        >  1:3647034/1‑150 (MQ=255)
                                                                                                 aGTGCGCGCTTGGGCCACCCGCGTGGGCGCCGGGGCGCTGGTGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGga                                                      <  1:872677/150‑1 (MQ=255)
                                                                                                   tGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGgaga                                                    <  2:5268213/150‑1 (MQ=255)
                                                                                                   ggcgcgcGTGGCCCCCCCGCCGGGGCGCGGTGGCGCTGCGGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGgaga                                                    <  1:1928452/147‑1 (MQ=255)
                                                                                                    gcgcgcTTGGGCCACCCGCCGGGGCGCCGTGGCGCGGCGGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGaa                                                   <  1:513891/148‑1 (MQ=255)
                                                                                                    gcgcgcTGGGGCCCCCCGCCTGGGCGCCGTGGCGGTGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGaa                                                   <  1:4663234/148‑1 (MQ=255)
                                                                                                            gggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCgg                                                       <  2:6464767/138‑1 (MQ=255)
                                                                                                            gggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCgg                                                       >  1:6463648/1‑138 (MQ=255)
                                                                                                             ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                                          <  2:1987857/150‑1 (MQ=255)
                                                                                                              gCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCCGGCCGCCTGGGGGTGGGAGGCCAACCCCGATAAGGCGCAGCGCGCGGGTCAGGCCGAGGACATCGGCAAGCAGTTGGCGGGGCAGGATcgcgcg                                         >  1:2463643/1‑150 (MQ=255)
                                                                                                               ccACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCGCTAAGGCGCAGCTCGCTGATCAGGCCGATGACTACGGCAAGCAGTTGGCGGAGAAGGATCgcgcgc                                        >  1:3299785/1‑150 (MQ=255)
                                                                                                                cccccGCGGGGGCGCGGGGGCGCTGCGGCGGGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCACGCGCa                                       <  1:3517087/148‑1 (MQ=255)
                                                                                                                 aCCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGGGTGGCCGGCCACCGCCTATAAGGCGCCGCTCGCGGATCAGGCCCATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCAc                                      >  1:1929232/1‑150 (MQ=255)
                                                                                                                    cGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCGGATCAGGCCGAGGACTTCGGCCAGCAGTTGGCGGAGAAGGATCGCGCGCACGGt                                   >  1:4714084/1‑150 (MQ=255)
                                                                                                                       ccgggcGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGCTGACTTCGGCAAGCAGTGGGCGGAGAAGGATCGCGCGCACGGTCTg                                >  1:7784420/6‑150 (MQ=255)
                                                                                                                       cTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTg                                <  2:3442258/150‑1 (MQ=255)
                                                                                                                         ggTCGCCGTGGCGCGGCTGCTGGCCGCCTGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATACGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTggg                              <  2:845727/150‑2 (MQ=255)
                                                                                                                            ggCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCg                           <  1:7631086/149‑1 (MQ=255)
                                                                                                                            cGCCGGGGCGCTGCTGCTGCCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCg                           <  2:4118553/150‑1 (MQ=255)
                                                                                                                             gCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCCGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg                          >  2:6222777/1‑150 (MQ=255)
                                                                                                                                gTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGgag                       <  2:8367317/150‑1 (MQ=255)
                                                                                                                                gTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCACCCAGGGTGTGGCAGGCCAACCCCTAAAATGCGCAGCTCGCGGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGgag                       >  1:7121434/1‑150 (MQ=255)
                                                                                                                                  ggCGCTGCGGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa                     <  1:774462/150‑1 (MQ=255)
                                                                                                                                     gctgctgctGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGc                  >  2:3328153/1‑150 (MQ=255)
                                                                                                                                       tgctgctgGCCGCCTGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTg                <  2:1223432/150‑1 (MQ=255)
                                                                                                                                       gctgcccggCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGTAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTg                >  2:6439158/8‑150 (MQ=255)
                                                                                                                                         ccgggcgccgccAGCGCAGGCG‑‑GCCCGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAATGCGCAGCACGCAGATCAGGCCGATGACTTAGGCAAGCACTTTGCGGATAACGATAGCGCGCACGGTCTGGAGCGGGAGAAAGCTgcagca           >  2:7941426/7‑150 (MQ=255)
                                                                                                                                          tgctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcag             >  1:6044331/1‑150 (MQ=255)
                                                                                                                                          ggcgggCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcag             <  1:221728/146‑1 (MQ=255)
                                                                                                                                            ctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagca           >  2:2643589/1‑150 (MQ=255)
                                                                                                                                             tgGCCGCCAGCGGTGGCGGCGCCCGTGCAGCCTGGGTGTGACAGGCCAACGCCTATTAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCTGTTGGCTGAGAATGATCGCGCGCACGGTCTGGTGCGGGAGAATGCTgcagcag          <  2:5360790/150‑1 (MQ=255)
                                                                                                                                              ggCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcagc         >  2:6941332/1‑150 (MQ=255)
                                                                                                                                               gccgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcagca        >  2:5030300/1‑150 (MQ=255)
                                                                                                                                               gccgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcagca        >  2:4278927/1‑150 (MQ=255)
                                                                                                                                                ccgcccgcGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGcagcagcag       <  1:18215/150‑1 (MQ=255)
                                                                                                                                                  gccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAgcg     >  2:115549/1‑150 (MQ=255)
                                                                                                                                                     aGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTc  <  2:3647034/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTC  >  NC_002947/3449102‑3449400

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