Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,452,278 A→C 11.8% 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,452,2780AC11.8% 424.0 / 6.1 219S18R (AGC→CGC) PP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base A (96/87);  new base C (16/17);  total (112/107)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

AAGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTC  >  NC_002947/3452129‑3452427
                                                                                                                                                     |                                                                                                                                                     
aaGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCa                                                                                                                                                       <  2:528815/150‑1 (MQ=255)
  gggCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAgc                                                                                                                                                     >  2:5672396/1‑150 (MQ=255)
   ggCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCcgcg                                                                                                                                                    >  1:2380894/1‑150 (MQ=255)
     cGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCt                                                                                                                                                  <  2:6258292/150‑1 (MQ=255)
      gACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTg                                                                                                                                                 <  2:5764301/150‑1 (MQ=255)
        cctCCTGGCGATCACCCCACGGCCCCAGGGAGGTAGCAACGGTCTGGAGGATCGCACGGCCCACGACACGCGAGCGCCCGAGGTGCTGCAGAGAGCGCTTGCGCCACCCGCCCGGTCGCCGTGGCGCAGCGGCTGGCCGCCAGCGCGggc                                                                                                                                               <  2:2620088/148‑1 (MQ=255)
          tCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCgcgggcgg                                                                                                                                             >  1:3052423/1‑150 (MQ=255)
          tCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcgg                                                                                                                                             <  2:3705258/150‑1 (MQ=255)
            cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACACGCGACCCCTCCAGGTGCTGCAGTTAGCGCTTGGGCCACCCGCCTGGTCGACGAGGCGCTGCAGCTGGCCGCCAGCGCTGgcggcg                                                                                                                                           >  2:3852814/1‑150 (MQ=255)
            cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGgcggcg                                                                                                                                           >  2:430905/1‑150 (MQ=255)
            cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGgcggcg                                                                                                                                           >  2:5207460/1‑150 (MQ=255)
            cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGgcggcg                                                                                                                                           >  2:2995183/1‑150 (MQ=255)
             tGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggcca                                                                                                                                          <  2:4423860/150‑3 (MQ=255)
             tGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGgcgc                                                                                                                                          >  1:509003/1‑150 (MQ=255)
                  aTCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGTATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCggg                                                                                                                                     <  2:475651/150‑1 (MQ=255)
                     aCCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAg                                                                                                                                  >  2:4696390/1‑150 (MQ=255)
                      ccAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGc                                                                                                                                 >  1:5064156/1‑150 (MQ=255)
                       cAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGCGGTGCTGCAGTGCGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGcc                                                                                                                                <  1:2235804/150‑1 (MQ=255)
                         aaCGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCCGgccgccag                                                                                                                              >  1:3507560/1‑148 (MQ=255)
                          aCGTATCAATGGTGGTAGCAACGTTCTGGAGTATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTgg                                                                                                                             <  2:5975924/150‑1 (MQ=255)
                          aCGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGCGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTgg                                                                                                                             <  2:6832204/150‑1 (MQ=255)
                          aCGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTgg                                                                                                                             >  1:4875343/1‑150 (MQ=255)
                            gTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGgt                                                                                                                           >  2:5121140/1‑150 (MQ=255)
                             tATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGgtg                                                                                                                          >  1:3361181/1‑150 (MQ=255)
                              accAATGGTGTTAGCAACGGTCTGGAGTATCGCACGGCCCTCTACACGCGAGCGCTGGAGGGGCTGCTGGGAGCGCTTGGGCCACCCGCGTGGTCGCCGTGGCGCTGCGGGTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGgtgt                                                                                                                         <  2:5738092/148‑1 (MQ=255)
                                cAATGGGGGTCGCACCGGTGGGGAGGATCGCAAGGCCCTCTCCACGCGCGCGCTCTAGGGGCTGCAGTGAGCGCTGGGGCCACCCGCCGGGCCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTgg                                                                                                                       <  2:5355480/150‑1 (MQ=255)
                                  agggtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACCCGCGCGCGCTCGCGGTGCTGCAGTGCGCGCTTGGGCCACCCGCCTGGTCGCGGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCa                                                                                                                     <  1:2211236/148‑1 (MQ=255)
                                  aTGGTGGTAGCAACGGTCTGGAGGATCGCACGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCa                                                                                                                     <  1:6820571/150‑1 (MQ=255)
                                  aTGGTGGTAGCAACGGTCGGGAGGATCGCAAGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCGGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCa                                                                                                                     <  2:4429388/150‑1 (MQ=255)
                                   tggtggTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAg                                                                                                                    >  2:6100561/1‑150 (MQ=255)
                                   tggtgcTAGCAACGGTCTGGAGTATCCCACGGCCCTCAACCAGCGCGCGCTCGAGGTGCTGCTGTGCGCGCTTGGCCCACCCGCCTGGTCGCCGGGGGGCGGCTGGTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAg                                                                                                                    <  1:3002517/150‑1 (MQ=255)
                                      tggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGcc                                                                                                                 >  2:3931658/1‑150 (MQ=255)
                                      gggtAGCAACGGACTGGAGGATCGCAAGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCATCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGcc                                                                                                                 >  2:1562626/2‑150 (MQ=255)
                                              aCGGTCGGGCGGAGCGCACGGCCCCCTACCCGCGCGCGCTCGAGGTGCTGCAGGGCGCGCTGGGGCCCCCCGCCTGGTCGCGGTGGCGGGGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCta                                                                                                         <  1:3275828/150‑1 (MQ=255)
                                                ggTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCCGGCCAACGCCgata                                                                                                       >  1:4326636/1‑150 (MQ=255)
                                                ggTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggc                                                                                                                                               >  1:1159915/1‑110 (MQ=255)
                                                 gTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATaa                                                                                                      <  2:1007407/150‑1 (MQ=255)
                                                  tCTGGAGTATCGCCCGGCCCTCTCCCCGCGCGCGCGCGGGGTGCTGGTGTGAGCGCTTGGGCCACGCGCCTGGGGGCGGTGGCGCTGCGGGTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAg                                                                                                     <  1:24389/150‑1 (MQ=255)
                                                   cTGGAGGATCGCAAGGCCCTCTACACGCGCGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAgg                                                                                                    <  2:4239151/150‑1 (MQ=255)
                                                     ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgcg                                                                                                  >  2:5459428/1‑150 (MQ=255)
                                                       aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGgcgcc                                                                                                >  1:3782953/1‑149 (MQ=255)
                                                       aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCa                                                                                                >  2:6671251/1‑150 (MQ=255)
                                                         gaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGc                                                                                              >  2:4468320/1‑150 (MQ=255)
                                                              cAAGGCCCTCTACACGCGAGCGCTGGGGGTGCTGCAGTGAGCGCGTGGGCCACCCGCGTGGGCGCGGGGGCGCTGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCt                                                                                         <  1:4318537/150‑1 (MQ=255)
                                                              cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCt                                                                                         >  2:2154632/1‑150 (MQ=255)
                                                              cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCCGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCt                                                                                         >  1:849667/1‑150 (MQ=255)
                                                               acggCCCTCTACAAGCGCGCGCTCGCGGTGCTGCAGTGCGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTg                                                                                        <  2:3052455/148‑1 (MQ=255)
                                                               aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCTATAAGGCGCAGCTCGCTg                                                                                        >  2:2329850/1‑150 (MQ=255)
                                                               aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTCTAAGGCGCAGCTCGCTg                                                                                        >  2:3487993/1‑150 (MQ=255)
                                                                aGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCGGCTCGCTGa                                                                                       >  1:4598563/1‑150 (MQ=255)
                                                                  gCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCGCGCTGATc                                                                                     >  2:736148/1‑150 (MQ=255)
                                                                    ccTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCTATAAGGCGCAGCGCGCTGATCAg                                                                                   >  1:4990423/1‑150 (MQ=255)
                                                                    ccTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGGGGGGCACGCCCACCCCGCTACGGCGCAGCGCGCTGATCAg                                                                                   >  1:5317783/1‑150 (MQ=255)
                                                                    ccTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCCGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCGCGCTGATCAg                                                                                   >  1:3104749/1‑150 (MQ=255)
                                                                    ccTCAACACGGGAGCGCTCGCGGTGCTGCCGGGAGCGCTTGGGCCACCCGCGTGGTCGCGGTGGGGCGGCTGCTGGCCGCCTGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAg                                                                                   <  2:878479/150‑1 (MQ=255)
                                                                     ctctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAgg                                                                                  >  1:3871662/1‑150 (MQ=255)
                                                                        tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCg                                                                               >  2:2800806/1‑150 (MQ=255)
                                                                        tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACCCCTATAAGGCGCAGCTCGCTGATCAGGCCg                                                                               >  1:5405244/1‑150 (MQ=255)
                                                                        tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCACCCTGGGTGTGGGAGGCCCACGCCTATAATGCGCAGCTCGCGGATAAGGccc                                                                               >  1:2873709/1‑149 (MQ=255)
                                                                         aCATGCGTGGGCTCGCGGTGCTGCAGTGAGCGCTTGGGCCCCGCGCGTGGTGGCTGTGGCGCTTCTGCTGGCGGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGa                                                                              <  1:5715247/150‑1 (MQ=255)
                                                                         aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGGGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCGGATCAGGCCGa                                                                              >  1:819105/1‑150 (MQ=255)
                                                                         aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCGGATCAGGCCGa                                                                              >  1:5804458/1‑150 (MQ=255)
                                                                          cacgcgtgcgCTCGAGGGGCTGCAGTGAGCGCTGGGGCCCCCCGCCTGGCCGCGGGGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGAt                                                                             <  1:2614115/147‑1 (MQ=255)
                                                                           aaGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCCGCCGGGGTGTGGCAGGCCCACCCCGATAAGGCGCCGCTCGCGGCTCAGGCCGCTg                                                                            >  1:895697/1‑150 (MQ=255)
                                                                           aaGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATg                                                                            >  2:1160059/1‑150 (MQ=255)
                                                                                tgcggTCGTGGTGCTGCAGTGAGCGCTTGGGCCACCCGCGTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGGGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTc                                                                       <  1:6195013/149‑1 (MQ=255)
                                                                                 gcgcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCg                                                                      >  2:6495325/1‑150 (MQ=255)
                                                                                   gcgggcggTGCTGCTGTGAGCGCTTGGGCCCCCCGGGTGGTCGCGGTGGCGCTGGTGGTGGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGc                                                                    <  1:2652647/144‑1 (MQ=255)
                                                                                    cTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCCATAAGGCGCAGCTCGCTGAGCAGGCCCATGACTTCGGCa                                                                   >  1:2563612/1‑150 (MQ=255)
                                                                                     tCGAGGTGCTGCTGTGAGCGCTTGGGCCACCCGCGTGGTCGCGGTGGGGCTGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCaa                                                                  <  1:1405275/150‑1 (MQ=255)
                                                                                      ggcggTGCTGCTGTGCGCGGGGGGGCCACCGGCCTGGTCGCGGGGGCGGTGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAg                                                                 <  1:2934558/147‑1 (MQ=255)
                                                                                           tgctgcAGGGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCGGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGtt                                                            <  2:3171444/150‑1 (MQ=255)
                                                                                            gctgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTAGAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTg                                                           >  2:4011327/1‑150 (MQ=255)
                                                                                            gctgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGGAGGGCCACGCCTAAAAGGCGCAGCTCGCTGATCAGGCCGATGACTTAGGCAAGCAGTTg                                                           >  1:1955059/1‑150 (MQ=255)
                                                                                             gtgcTGTGAGCGCTGGGGCCACCCGCGTGGGCGGGGGGGCGGGGGTGGTGGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTgg                                                          <  1:5249706/149‑1 (MQ=255)
                                                                                              tgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGc                                                         <  2:6225204/150‑1 (MQ=255)
                                                                                                 agggAGCGCTTGGGCCACCCGCATGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGTCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGga                                                      <  2:4609094/147‑1 (MQ=255)
                                                                                                 aGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGga                                                      <  2:3640610/150‑1 (MQ=255)
                                                                                                  gTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGgag                                                     >  2:1533077/1‑150 (MQ=255)
                                                                                                          ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcg                                             >  2:5535591/1‑150 (MQ=255)
                                                                                                            gggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcg                                           >  2:6209761/1‑150 (MQ=255)
                                                                                                                cACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCa                                       >  2:6409463/1‑150 (MQ=255)
                                                                                                                  cccGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACg                                     <  2:2441411/150‑1 (MQ=255)
                                                                                                                         ggTCGCCGTGGCGCTGCTGCTGGCCGCCTGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGa                              <  2:5064156/150‑1 (MQ=255)
                                                                                                                         ggTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGa                              <  1:6689044/150‑1 (MQ=255)
                                                                                                                             gCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGAAATCGGCAAGCACTTGGCGGAGAAGGATCGCGCCCACGGTCTGGAGCgg                          >  1:3597410/1‑150 (MQ=255)
                                                                                                                              ccGTGGCGCTGCTGGTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCggg                         <  1:1667900/150‑1 (MQ=255)
                                                                                                                                 tGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGgaga                      <  2:5762997/150‑1 (MQ=255)
                                                                                                                                    cgcTGCTGCGGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAgg                   <  1:100568/150‑1 (MQ=255)
                                                                                                                                     gctgctgctGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGc                  >  2:5220227/1‑150 (MQ=255)
                                                                                                                                      ctgctgcggGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCt                 <  2:3611140/150‑1 (MQ=255)
                                                                                                                                       ggctggtgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTg                <  2:2472806/149‑1 (MQ=255)
                                                                                                                                         ccgggcgccgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAAGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGATGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgca              >  1:4698850/7‑150 (MQ=255)
                                                                                                                                          tgctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcag             >  2:4300206/1‑150 (MQ=255)
                                                                                                                                             gcgccgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCCAGCAGTTGGCGGAGaa                                                   >  1:6157278/3‑109 (MQ=255)
                                                                                                                                              cgccgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcagc         >  2:6282124/2‑150 (MQ=255)
                                                                                                                                                   ccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAgcgc    <  1:5591337/150‑1 (MQ=255)
                                                                                                                                                   ccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAgcgc    >  1:5903793/1‑150 (MQ=255)
                                                                                                                                                    cAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCt   >  2:201415/1‑150 (MQ=255)
                                                                                                                                                    aaGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCACGCAGTTGGCGGAGAAGGATCGCGCGCCCGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCt   >  2:5428457/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTC  >  NC_002947/3452129‑3452427

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: