Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,449,251 A→C 13.4% 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,2510AC13.4% 237.3 / 4.2 126S18R (AGC→CGC) PP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base A (61/44);  new base C (12/8);  total (73/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

AAGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGC  >  NC_002947/3449102‑3449398
                                                                                                                                                     |                                                                                                                                                   
aaGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCa                                                                                                                                                     <  2:3450629/150‑1 (MQ=255)
aaGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCa                                                                                                                                                     <  2:5179159/150‑1 (MQ=255)
aaGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCa                                                                                                                                                     >  2:2289364/1‑150 (MQ=255)
  gggCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGccgcccgc                                                                                                                                                   >  1:3331029/1‑150 (MQ=255)
       aCATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGAGCTGCTGCTGGCCGCCAGCGCTgg                                                                                                                                              <  2:3874697/150‑1 (MQ=255)
        cATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGggc                                                                                                                                             >  1:1561635/1‑150 (MQ=255)
        cATCCGGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggc                                                                                                                                             <  1:1842667/150‑1 (MQ=255)
          tCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcgg                                                                                                                                           <  2:1700913/150‑1 (MQ=255)
            cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGgcggcg                                                                                                                                         >  1:1503689/1‑150 (MQ=255)
            cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGgcggcg                                                                                                                                         >  1:1861368/1‑150 (MQ=255)
             tGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGgcgc                                                                                                                                        <  2:5262205/150‑1 (MQ=255)
                cGATCACCAAACGTATCAATGGGGGTAGCAAGGGTGTGGTGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGTGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCg                                                                                                                                     <  1:2679669/150‑1 (MQ=255)
                     aCCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAg                                                                                                                                >  1:3939281/1‑150 (MQ=255)
                     aCCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAg                                                                                                                                >  2:1611432/1‑150 (MQ=255)
                     aCCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCATGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAg                                                                                                                                >  2:3969816/1‑150 (MQ=255)
                      ccAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGc                                                                                                                               >  1:4546587/1‑150 (MQ=255)
                       cAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGcc                                                                                                                              >  2:1302640/1‑150 (MQ=255)
                        aaaCGTATCTCTGGGGGTAGCACGGGTCGGGAGGAGCGCACGGCCCTCTACACGCGAGCGCTCGTGGTGCTGCAGTGAGCGCTGGGGCCACCCGCGTGGGCGCCGGGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCt                                                                                                                             <  1:486391/150‑1 (MQ=255)
                        aaaCGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCt                                                                                                                             >  2:4699006/1‑150 (MQ=255)
                          aCGTATCATTGGTGTTAGAAACGGTCTGGAGGATCGCATGGCCCTCTACAAGCTAGCGCTCGTGGTGCTGCTGTGAGGGGTTGGGCCACCGGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTgg                                                                                                                           <  2:1407337/150‑1 (MQ=255)
                              aTCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCACCAGCGCTGGCGGCGCCCGGGCAGCCTGGgtgt                                                                                                                       >  1:1865306/1‑150 (MQ=255)
                                 aaTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGc                                                                                                                    >  1:2223460/1‑150 (MQ=255)
                                  ggtgccatAGCACCGGTCTGGAGGAGCGCAAGGCCCTCAACAAGCGCGCGCTCGCGGTGCTGCTGTGCGGGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCa                                                                                                                   <  1:4278949/143‑1 (MQ=255)
                                  aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCa                                                                                                                   <  2:3729692/150‑1 (MQ=255)
                                   tggtgtTAGCAACGTTCTGTTGGATCGAATGGCCCTCTACAAGCGTGCGCTCTAGGTGCTGCATTTAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGGTGGCCGCATGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAg                                                                                                                  <  2:198485/150‑1 (MQ=255)
                                   tggtggTCGCAAGGGTCTGGATGATCGCACGGCCCTCTACAAGCGAGCGCTGGTGTTGCCGCTTTGAGCGCTTGGGCCACCCGCCTGGGCGCGGTGGCGCTGCTGCTGGCCGCCAGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAg                                                                                                                  <  2:691329/150‑1 (MQ=255)
                                     gtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGc                                                                                                                >  1:584468/1‑150 (MQ=255)
                                     gtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGc                                                                                                                >  2:2061510/1‑150 (MQ=255)
                                     gtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCGACCCGCGAGCGCTCGAGGTGCTGCTGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGCGGCGCTGCTGCGGCCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGc                                                                                                                <  2:3193369/150‑1 (MQ=255)
                                      tggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCGGGcc                                                                                                               >  1:1480733/1‑150 (MQ=255)
                                      tggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGcc                                                                                                               >  1:2016000/1‑150 (MQ=255)
                                         tAGCACCGGGCTGGCGGATCGCAAGGCCCTCTACACGCGCGCGCTCGCGGTGCGGCTGGGAGCGCTGGGGCCACCCGCGTGGGCGCCGTGGCGCTGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAAc                                                                                                            <  2:485565/150‑1 (MQ=255)
                                         tAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAAc                                                                                                            >  2:5351824/1‑150 (MQ=255)
                                          aGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACg                                                                                                           >  2:4736181/1‑150 (MQ=255)
                                           gCAACGGTCGGGAGGATCGCACGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCCGGGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGc                                                                                                          <  2:5379368/150‑1 (MQ=255)
                                               cGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCtat                                                                                                      <  2:3844872/150‑1 (MQ=255)
                                                   gtGGAGGATCGCAAGGCCCTCTACACGCGAGCGCTCGTGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAgg                                                                                                  <  2:1882791/149‑1 (MQ=255)
                                                      gtggaTCGCACGGCCCTCTACACGCGAGCGCTCGCGGTGCTGCTGTGTGCGCTTGGGCCACCCGCGTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgcgc                                                                                               <  1:1334076/148‑1 (MQ=255)
                                                      gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGCCAACGCCAATAAGgcgc                                                                                               >  1:4922425/1‑150 (MQ=255)
                                                      gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCGATAAGgcgc                                                                                               >  1:2661377/1‑150 (MQ=255)
                                                      gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgcgc                                                                                               >  2:4968816/1‑150 (MQ=255)
                                                       aggaTCGCAAGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCa                                                                                              <  2:2635222/150‑1 (MQ=255)
                                                           tCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggccag                                                                                                                                       <  1:3969365/105‑4 (MQ=255)
                                                           tCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggccag                                                                                                                                       >  2:3969365/1‑102 (MQ=255)
                                                              cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCt                                                                                       >  2:455093/1‑150 (MQ=255)
                                                                aGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGAGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGACCAACGCCAATAAGGAGCAACTCGCTGa                                                                                     >  1:5160461/1‑150 (MQ=255)
                                                                 ggCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGAt                                                                                    >  2:2674610/1‑150 (MQ=255)
                                                                     ctctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCAGCCAGGGTGTGGCAGGCCCACGCCGATAAGGCGCCCCTCGCTGATCAgg                                                                                >  1:2952912/1‑150 (MQ=255)
                                                                      tctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGAAAAGGCGCCGCACGCAGATCAGGc                                                                               >  1:4273000/1‑150 (MQ=255)
                                                                        tACACGCGCGCGCGCGTGGTGCTGCTGGGCGCGCTTGGGCCACCCGCCTGGTCGCGGGGGGGCGGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCg                                                                             <  1:4215268/150‑1 (MQ=255)
                                                                         aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCGGGGGGTGGCAGGCCCACGCCTATAAGGCGCAGCTCGCGGATCAGGCCGa                                                                            >  1:3055595/1‑150 (MQ=255)
                                                                         aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCGGGTCGCGGTGGGGCGGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGa                                                                            <  2:4230737/150‑1 (MQ=255)
                                                                          cacgcgagcgCGCGCGGGGCGGCGGTGAGCGGTTGGGCCACCCGCCGGGGCGCGGGGGCGCTGCGGCTGGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGAt                                                                           <  1:3835739/147‑1 (MQ=255)
                                                                          cAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAATGCGCAGCTCGCTGATCAGGCCGAt                                                                           >  1:4314951/1‑150 (MQ=255)
                                                                             gcgagcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggc                                                                                                                                          <  1:212892/84‑1 (MQ=255)
                                                                             gcgagcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCTAAAAGGCGCAGCTCGCTGATCAGGCCGAGGAc                                                                        >  1:1441583/1‑150 (MQ=255)
                                                                              cgagcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACt                                                                       >  2:2677815/1‑150 (MQ=255)
                                                                                cgcgcTCGCGGTGCTGCTGTGCGCGCTGGGGCCACCCGCGTGGGCGCGGGGGCGGTGCTGCGGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTc                                                                     <  1:1013977/149‑1 (MQ=255)
                                                                                   gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCCGGCCAACCCCGATAAGGCGCAGCGCGCGGATCAGGCCGAGGACTTCGGc                                                                  >  1:1119366/1‑150 (MQ=255)
                                                                                     tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGGGAAGGCGCAGCTCGCTGATCCGGCCGATGACGTCGGCaa                                                                >  2:4140304/1‑150 (MQ=255)
                                                                                             ctgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTgg                                                        <  1:1319288/150‑1 (MQ=255)
                                                                                              tgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGc                                                       >  2:242002/1‑150 (MQ=255)
                                                                                              tgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGc                                                       >  1:4583588/1‑150 (MQ=255)
                                                                                               gccgTGCGCGCTGGGGCCACCCGCCGGGTCGCCGTGGCGGGGCTGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCg                                                      <  2:2313147/147‑1 (MQ=255)
                                                                                               gcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGAGGCAGGCCAACGCCGATAAGGCGCAGCTCCCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCg                                                      >  2:3658491/1‑150 (MQ=255)
                                                                                               gcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGGTGGCg                                                      >  1:3166164/1‑150 (MQ=255)
                                                                                                  gTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCCACCCCGATAATGCGCAGCGCGCGGATAAGGCCGCGGACGTCGGCAAGCAGTTGGCGgag                                                   >  1:485705/1‑150 (MQ=255)
                                                                                                   tGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGgaga                                                  >  2:778469/1‑150 (MQ=255)
                                                                                                    gAGCGCTTGGGCCACCCGCCGGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGaa                                                 <  2:4201530/150‑1 (MQ=255)
                                                                                                    gAGCGCTGGGGCCCCCCGCCTGGTCGCCGTGGCGCTGCGGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGaa                                                 <  1:2663204/150‑1 (MQ=255)
                                                                                                     tgcgcTGGGGCCACCCGCCTGGTCGCCGTGGCGCGGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAg                                                <  2:4583588/149‑1 (MQ=255)
                                                                                                     aGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAg                                                <  2:369284/150‑1 (MQ=255)
                                                                                                     aGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAg                                                >  1:3141643/1‑150 (MQ=255)
                                                                                                     aGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCGGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAg                                                >  1:3141694/1‑150 (MQ=255)
                                                                                                          ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcg                                           >  2:3943601/1‑150 (MQ=255)
                                                                                                           tGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgc                                          <  2:123384/150‑1 (MQ=255)
                                                                                                            gggCCACCCGCCTGGTCGCCGTGGAGCTGCTGCTGGGCGCCAGCGCTGGGGGAGCCCGGGCAGCCTGGGTGTGGCAGGCCAAAGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGTTagcg                                         >  2:1380276/1‑150 (MQ=255)
                                                                                                             ggCCACCCGCCGGGTCGCCGGGGCGCGGGTGCGGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                                        <  1:3344706/150‑1 (MQ=255)
                                                                                                             gcccACCCCCCTGCTCGCCGCGTCGCTGCTGCTGCCGGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAACGCTCAGCTCGCTGATCAGGCCGTTGACTTCGGCAAGCAGTTGGCGGATAAGGATcgcgc                                        <  2:1168768/148‑1 (MQ=255)
                                                                                                              gCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACGCCTAGAAGGCGCCGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgcg                                       >  2:2861624/1‑150 (MQ=255)
                                                                                                              gCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgcg                                       >  2:4522998/1‑150 (MQ=255)
                                                                                                                    cGCCTGGGCGCCGGGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGt                                 <  2:1607025/150‑1 (MQ=255)
                                                                                                                     gCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTAGGCGGAGAAGGATCGCGCGCACGGTc                                >  1:4427233/1‑150 (MQ=255)
                                                                                                                       cggggcGCCGTGGCGCGGCGGCGGGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTg                              <  1:4699949/145‑1 (MQ=255)
                                                                                                                       cTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGAGGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCCCGGTCTg                              >  1:4402583/1‑150 (MQ=255)
                                                                                                                         ggTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGa                            <  1:3879765/150‑1 (MQ=255)
                                                                                                                                  ggCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa                   <  1:3042052/150‑1 (MQ=255)
                                                                                                                                  ggCGCTGCTGCTGCCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa                   <  1:3563464/150‑1 (MQ=255)
                                                                                                                                        gcggcgggCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgc             <  2:3331029/144‑1 (MQ=255)
                                                                                                                                               gccgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcagca      >  1:1146928/1‑150 (MQ=255)
                                                                                                                                                ccgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGcagcagcag     >  2:3019579/1‑150 (MQ=255)
                                                                                                                                                 cgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGc                                                                               <  1:379118/75‑1 (MQ=255)
                                                                                                                                                 cgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGc                                                                               >  2:379118/1‑75 (MQ=255)
                                                                                                                                                 cgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCg                                                      <  1:3508531/100‑1 (MQ=255)
                                                                                                                                                 cgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCg                                                      >  2:3508531/1‑100 (MQ=255)
                                                                                                                                                  gccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAgcg   >  1:677936/1‑150 (MQ=255)
                                                                                                                                                  gccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAgcg   >  2:5393155/1‑150 (MQ=255)
                                                                                                                                                   ccAGCGCTGGCGGAGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGATGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAgcgc  >  2:532917/1‑150 (MQ=255)
                                                                                                                                                     aGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCg                                                      >  2:1170710/1‑96 (MQ=255)
                                                                                                                                                     aGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCg                                                      <  1:1170710/96‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AAGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGC  >  NC_002947/3449102‑3449398

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: