Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,449,256 T→G 18.0% A19A (GCT→GCG PP_5556 → bacteriophage lysis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,449,2560TG18.0% 315.9 / 28.0 180A19A (GCT→GCGPP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base T (73/67);  new base G (10/26);  total (86/94)
Fisher's exact test for biased strand distribution p-value = 1.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

ATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATC  >  NC_002947/3449111‑3449404
                                                                                                                                                 |                                                                                                                                                    
aTCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcg                                                                                                                                                  <  1:2122567/150‑1 (MQ=255)
 tCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcgg                                                                                                                                                 <  2:3811090/150‑1 (MQ=255)
  ccTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggc                                                                                                                                                <  1:3037355/150‑1 (MQ=255)
   cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGgcggcg                                                                                                                                               >  2:2197073/1‑150 (MQ=255)
    tGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGgcgc                                                                                                                                              <  2:1329078/150‑1 (MQ=255)
    tGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGgcgc                                                                                                                                              >  2:2229055/1‑150 (MQ=255)
       cGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCTGGCGGCGCCCg                                                                                                                                           >  1:782384/1‑150 (MQ=255)
        gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCgg                                                                                                                                          >  1:3000256/1‑150 (MQ=255)
        gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCgg                                                                                                                                          >  2:3886040/1‑150 (MQ=255)
         aTCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCggg                                                                                                                                         <  2:3692060/150‑1 (MQ=255)
           cACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCa                                                                                                                                       <  2:2989526/150‑1 (MQ=255)
            aCCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGGGccg                                                                                                                                      >  1:3197924/1‑148 (MQ=255)
             ccAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGc                                                                                                                                     >  2:1056347/1‑150 (MQ=255)
              cAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGcc                                                                                                                                    >  1:127467/1‑150 (MQ=255)
               aaaCGTATCAATGGTGGTAGCAACGGTCTTGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGAGCTGCTGCTGGCCGCCAGCCCAGGAGGCGCCCCGGCAGCCt                                                                                                                                   >  2:1739052/1‑150 (MQ=255)
                 aCGTATCTTTGGTGGTAGCATGGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGTGGTGCTGCTGTGCGCGGTTGGGCCACCCGCCTGGTCGCGGTGGCGCTGCTGCTGGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTgg                                                                                                                                 <  1:3242088/150‑1 (MQ=255)
                 aCGTATCATTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGCGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCGTGGTCGCGGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTgg                                                                                                                                 <  1:1045279/150‑1 (MQ=255)
                     ctCAAGGGTGGTAGCAACGGTCTGGAGGATCGCACGGCCCTCTACAAGCGCGCGCTCGAGGTGCTGCAGGGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGGGCTGCTGCTGGCCGCCTGCGCTGGCGGCGCCCGGGCAGCCTGGgtgt                                                                                                                             <  2:3678875/149‑1 (MQ=255)
                     atctttggtggTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCTGTGCGCGCTTGGGCCACCCGCCTGGTCGCGGTGGGGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGgtgt                                                                                                                             <  1:2146218/145‑1 (MQ=255)
                       cattggtggTAGCAAGGGTCTGGAGGATCGCAAGGCCCTCTACACGCGAGCGCTCGTGGTGCTGCCGTGAGCGCTTGGGCCCCCCGCCTGGTCGCGGTGGCGCTGGTGCGGGCCGCCCGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTgg                                                                                                                           <  1:2406279/147‑1 (MQ=255)
                       cAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTgg                                                                                                                           <  2:734844/150‑1 (MQ=255)
                        aaTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGGGGCGCCCGGGCCGCCGGGGGGTGGc                                                                                                                          >  1:26565/1‑150 (MQ=255)
                          ggggggtAGCACCGGTCGGGAGGAGCGCACGGCCCTCGACACGCGCGCGCTCGCGGTGCGGCGGGGAGCGCTTGGGCCACCCGCCTGGTCGCCGGGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAg                                                                                                                        <  2:2146001/146‑1 (MQ=255)
                            gtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGc                                                                                                                      >  2:2274432/1‑150 (MQ=255)
                              ggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCa                                                                                                                    >  2:1993253/1‑150 (MQ=255)
                                gaGCAACGGTCTGGAGGATCGCACGGCCCTCTCCCCGCGAGCGCGCGCGGTGCTGCAGTGAGCGCTGGGGCCACCCGCCTGGTCGCCGTGGCGCTGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAAc                                                                                                                  <  1:279899/149‑1 (MQ=255)
                                  gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGc                                                                                                                >  2:4091429/1‑150 (MQ=255)
                                   cAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAAcccc                                                                                                               >  1:2786829/1‑150 (MQ=255)
                                   cAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCGCCCGGGCCGCCTCGGGGTGGCAGGCCAAcccc                                                                                                               >  1:1237629/1‑150 (MQ=255)
                                    aaCGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGGAGGCCAACGCCt                                                                                                              >  1:343034/1‑150 (MQ=255)
                                    aaCGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCt                                                                                                              >  2:4142775/1‑150 (MQ=255)
                                    aaCGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGGGTGGCAGGCCAACGCCt                                                                                                              >  2:3778191/1‑150 (MQ=255)
                                      cGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCGGGCGGCGCCCGGGCCGCCGGGGGGGGGCAGGCCCACCCCgat                                                                                                            >  1:3844431/1‑150 (MQ=255)
                                      cGGTCTGGAGGAGCGCACGGCCCTCTACAAGCGAGCGCTCGCGGTGCTGCCGTGCGCGCTGGGGCCACCCGCCTGGACGCCGTGGCGGTGCTGCTGGCCGCCAGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCtat                                                                                                            <  2:1190330/150‑1 (MQ=255)
                                        gTCTGGTGGATCGCACGGCCCTCCCCACGCGTGCGCTCGTGGTGCTGCAGGGCGCGCTTGGCCCACCCGCGTGGTCGCCGTGGGGGTGGGGGGGGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATaa                                                                                                          <  1:3429825/150‑1 (MQ=255)
                                         tcgggaggaTCGCAAGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAg                                                                                                         <  2:4309456/147‑1 (MQ=255)
                                             gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgcgc                                                                                                     >  2:728168/1‑150 (MQ=255)
                                              aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCa                                                                                                    >  2:4343984/1‑150 (MQ=255)
                                                      aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGCCAACGCCAAAAAGGCGCAGCTCGCTg                                                                                            >  1:1580294/1‑150 (MQ=255)
                                                           ccTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCGGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGCCAACGCCAAAAAGGCGCAGCGCGCGGATCAg                                                                                       >  1:1884817/1‑150 (MQ=255)
                                                            ctctACACGCGAGCGCGCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAgg                                                                                      <  2:1922951/150‑1 (MQ=255)
                                                             tctACAAGCGAGCGCTCGAGTTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGc                                                                                     <  2:815612/150‑1 (MQ=255)
                                                             tctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGc                                                                                     >  2:3447888/1‑150 (MQ=255)
                                                               tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCCCCCCGCCAGCCAGGGTGTGGCAGCCCAACCCCGATAAGGCGCAGCTCGCTGAGCAGGccc                                                                                   >  1:3235563/1‑149 (MQ=255)
                                                                    gcgagcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGAc                                                                              >  2:1284367/1‑150 (MQ=255)
                                                                      gagcgGTCGAGGTGCTGCTGGGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGGGCTGCTGCGGGCCGCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACtt                                                                            <  2:2378113/150‑1 (MQ=255)
                                                                      gagcgCTCGTGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACtt                                                                            <  1:4510761/150‑1 (MQ=255)
                                                                          gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCGGGCCGCCTGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGAAGTCATCGGc                                                                        <  2:3002591/150‑1 (MQ=255)
                                                                            tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCaa                                                                      <  2:2188429/150‑1 (MQ=255)
                                                                            tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCaa                                                                      >  2:2880996/1‑150 (MQ=255)
                                                                            tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCaa                                                                      >  2:49359/1‑150 (MQ=255)
                                                                            tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCaa                                                                      >  2:3156463/1‑150 (MQ=255)
                                                                              gAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGc                                                                    >  1:2606209/1‑150 (MQ=255)
                                                                              gAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCGGATCAGGCCGATGACGTCGGCAAGc                                                                    >  1:2606073/1‑150 (MQ=255)
                                                                               tggTGCGGCTGTGTGCGCTTGGGCCACCCGCCTGGTCGCGGTGGGGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCa                                                                   <  1:4311873/149‑1 (MQ=255)
                                                                                     tgcTGTGCGCGCTTGGGCCCCCCGCCTGGTCGCCGTGGCGCGGCTGCTGGCCGCCTGGGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGc                                                             <  1:441073/150‑1 (MQ=255)
                                                                                     tgcCGTGCGCGGTTGGGCCCCCCGGGGGGGCGCCGTGGCGGGGGGGGGGGCCGCCCGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGc                                                             <  1:1970934/150‑1 (MQ=255)
                                                                                     tgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGc                                                             >  2:4193841/1‑150 (MQ=255)
                                                                                       cAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCgg                                                           >  2:1609950/1‑150 (MQ=255)
                                                                                        tgTGCGCGGTGGGGCCCCCCGCCTGGTCGCGGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGga                                                          <  1:1000297/149‑1 (MQ=255)
                                                                                        aGTGAGCGCTTGGGCCACCCGCCTGGTCGCGGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGga                                                          <  1:647280/150‑1 (MQ=255)
                                                                                        aGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCggc                                                          >  1:2729751/1‑149 (MQ=255)
                                                                                          ggcgcgcTGGGGCCACCGGGGTGGGGGGGGTGGCGCTGCTGCGGGCCGCGCGGGCGGGGGTCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGTCTTCGGCAAGCAGTTGGCGgaga                                                        <  1:1363293/147‑1 (MQ=255)
                                                                                           gAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTGGGCGGAGaa                                                       <  2:1204594/150‑1 (MQ=255)
                                                                                            aGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAg                                                      <  2:545904/150‑1 (MQ=255)
                                                                                             ggggttGGGCCACCCGCGGGGTCGGGGTGGGGCTGGTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAgg                                                     <  1:4229479/146‑1 (MQ=255)
                                                                                              cgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGa                                                    <  2:3289513/150‑1 (MQ=255)
                                                                                              cgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGa                                                    >  2:3054732/1‑150 (MQ=255)
                                                                                                cTTGGGCCACCCGCGTGGTCGCGGGGGGGCGGCGGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATc                                                  <  1:3024689/150‑1 (MQ=255)
                                                                                                cTTGGGCCACCCGCGTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATc                                                  <  1:643613/150‑1 (MQ=255)
                                                                                                cTTGGGCCACCCGCCTGGTCGCCGTGGCGCGGCTGGGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATc                                                  <  1:28394/150‑1 (MQ=255)
                                                                                                  tGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgc                                                <  1:3348374/150‑1 (MQ=255)
                                                                                                          cggcggggTCGCGGGGGCGCTGCTGCGGGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACgg                                        <  1:4265775/144‑1 (MQ=255)
                                                                                                          ccGCGTGGTCGGCGTGGGGGTGCTGCTGGGGGGCAGGGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCGTTTAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCATGCAGTTGGCGGAGAAGGATCGCGCGCACgg                                        <  1:3355918/150‑1 (MQ=255)
                                                                                                          ccGCCTGGTCGCCGTGGCGGTGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCTGGCCGTTGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACgg                                        <  2:364967/150‑1 (MQ=255)
                                                                                                          ccGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCCAGCAGGTGGCGGGGAAGGATCGCGCGCACgg                                        >  1:1714091/1‑150 (MQ=255)
                                                                                                            gcgtGGTCGCGGTGGCGGTGCTGGTGGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTc                                      <  1:3978575/147‑1 (MQ=255)
                                                                                                              gtGGTCGGGGTGGCGCTGCTGCTGGCCGCCTGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTg                                    <  1:3597160/149‑1 (MQ=255)
                                                                                                                ggTCGGGGTGGGGCTGCTGCTGGCCGCCCGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGa                                  <  1:3627103/150‑1 (MQ=255)
                                                                                                                    gCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGCTGACTTCGGCAAGCACTTGGCGGAGAAGGATAGCGCGCACGGTCTGGAGCgg                              >  1:3481674/1‑150 (MQ=255)
                                                                                                                    gCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAATGCGCAGCTAGCAGATCAGGCCGATGACTACGGCAAGCAGTTGGCGGAGAATGATCGCGCGCACGGTCTTGAGCgg                              >  1:2196470/1‑150 (MQ=255)
                                                                                                                      cGTGGCGCTGCTGCTGGCCGCCAGCGCTTGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAAAGCCTATAATGCGCAGCTCGCTGATCAGGCCGATTAATTAGGCAATCAGTTTGCGGAGAAGTATAGCGCGCACGGTCTGGAGCGGga                            >  1:2934528/1‑150 (MQ=255)
                                                                                                                      cGTGGCGCTGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGga                            <  2:1172472/150‑1 (MQ=255)
                                                                                                                         ggCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa                         <  1:1404951/150‑1 (MQ=255)
                                                                                                                         ggCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa                         <  2:1551232/150‑1 (MQ=255)
                                                                                                                         ggCGCTGCGGCTGGCCGCCACCGCTGCCGGCGCCCGGGCAGCCTGCGTGTGGCTGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGCAGAAGGATCGCGCGAACGGTCTGGAGCGGGAGaa                         <  2:1551241/150‑1 (MQ=255)
                                                                                                                          gcgcGGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAg                        <  2:574756/150‑1 (MQ=255)
                                                                                                                           cgcgggtgcggGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGACGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAgg                       <  1:3149729/144‑1 (MQ=255)
                                                                                                                            gctgctgctGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGc                      >  2:3219127/1‑150 (MQ=255)
                                                                                                                             ctgctgctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCt                     <  1:376615/150‑1 (MQ=255)
                                                                                                                                ctgctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgca                  >  2:2331707/1‑150 (MQ=255)
                                                                                                                                ctgcggGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgca                  <  1:63249/150‑1 (MQ=255)
                                                                                                                                cggctgGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgca                  <  1:3262217/148‑1 (MQ=255)
                                                                                                                                 ggcgggCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcag                 <  1:457656/146‑1 (MQ=255)
                                                                                                                                  gctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCTGGAGAAGGCTgcagc                >  1:92344/1‑150 (MQ=255)
                                                                                                                                    tgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcag              >  2:3745532/1‑150 (MQ=255)
                                                                                                                                    tgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcag              >  1:4124730/1‑150 (MQ=255)
                                                                                                                                        cgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCagcagcagc          <  1:2519930/150‑1 (MQ=255)
                                                                                                                                          ccAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAgcgc        <  1:305912/150‑1 (MQ=255)
                                                                                                                                                cTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATc  <  2:2136060/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATC  >  NC_002947/3449111‑3449404

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: