Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,406,2490GA25.0% 45.2 / 10.6 28F284F (TTC→TTT) puuPputrescine permease
Reads supporting (aligned to +/- strand):  ref base G (10/11);  new base A (7/0);  total (17/11)
Fisher's exact test for biased strand distribution p-value = 2.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGAATGCAGGCGGGCGAAGACGCTGGCCGGGATCACGTTGTCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGATGAAGTACGACACGGTGATGAACACCACACCGCCAATCAGTGCGGTGAG  >  NC_002947/1406110‑1406386
                                                                                                                                           |                                                                                                                                         
gcgAATGCAGGCGGGCGAAGACGCTGGCCGGGATCACGTTGTCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCtt                                                                                                                                   <  2:32813/148‑1 (MQ=255)
                gAAGACGCTGGCCGGGATCACGTTGTCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAgg                                                                                                                   <  8:145506/148‑1 (MQ=255)
                  aGACGCTGGCCGGGATCACGTTGTCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGa                                                                                                                                          >  7:111794/1‑123 (MQ=255)
                   gACGCTGGCCGGGATCACGTTGTCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGaa                                                                                                                                         <  8:111794/123‑1 (MQ=255)
                        tGGCCGGGATCACGTTGTCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTc                                                                                                           >  3:113102/1‑148 (MQ=255)
                             gggATCACGTTGTCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCa                                                                                                      >  1:113344/1‑148 (MQ=255)
                                  cACGTTGTCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGGGAGGCCAGGCCGGAGGCGATGGTGTTGATGTCGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGGCGCCAACGTACAGGGCAATTTCCGGCAGGgct                                                                                                 <  1:179019/148‑1 (MQ=255)
                                  cACGTTGACCCGGCCCATTACATAAAGCAGGCGTCACACGCTGGTTTGCTAGGCCAGGCCCGAGGCGATGGTGATGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACACGGCAATTTCCGGCAGGgct                                                                                                 <  4:99277/148‑1 (MQ=255)
                                          cgcgGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTg                                                                                         <  3:51815/148‑1 (MQ=255)
                                           gcgGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTAGATCACGGTGCACGCGATAAAAATAAAATTAAAAAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTgg                                                                                        >  6:33014/1‑148 (MQ=255)
                                                 cATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTAGATCACGGTGCTCGCCATAAAAATAAAAAAAAAAAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGt                                                                                    >  4:64953/1‑146 (MQ=255)
                                                   ttACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGaa                                                                                >  5:162998/1‑148 (MQ=255)
                                                     acatacaGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATAGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGc                                                                              >  2:60166/1‑148 (MQ=255)
                                                     acatacaGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTAGATCACTGTGCACGCCATAAAAATAAAAAAAAAAAGATTGCCGCCAAAGTACAGGGCAATTTCCGGCAGGTCTGCTTCCTGGTCGTGGAAGc                                                                              >  6:35216/1‑148 (MQ=255)
                                                     acatacaGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTAGATCACGGTGCACGCCATAAAAATAAAAAAAAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGc                                                                              >  6:35215/1‑148 (MQ=255)
                                                       atacaGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCgg                                                                            >  7:33337/1‑148 (MQ=255)
                                                                  gTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTAGATCACGGTGCACGCCATAAAAATAAAAAAAAAGAGCTAGCCGCCAACGTCCAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCg                                                                  >  6:94526/1‑147 (MQ=255)
                                                                      cacGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGaa                                                              <  1:37668/147‑1 (MQ=255)
                                                                                   gAGGCCAGGCCCGAGGCGATGGTGTAGATCACGGTGCACGCCATAAAAATAAAAAAAAAAAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGAt                                                  >  4:11588/1‑146 (MQ=255)
                                                                                    aggccaggccCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGATGaa                                               <  6:77179/148‑1 (MQ=255)
                                                                                     ggccaggccCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGATGAAg                                              <  2:143099/148‑1 (MQ=255)
                                                                                         aggccCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGATGAAGTACg                                          >  5:100782/1‑148 (MQ=255)
                                                                                         aggccCGAGGCGATGGTGTAGATCACGGTGCACGCCATAAAAATAAAAAAAAAAATCTTGACGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTGGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGATGAAGTACg                                          >  2:46018/1‑148 (MQ=255)
                                                                                                  gCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGATGAAGTACGACACGGtga                                 >  5:145974/1‑148 (MQ=255)
                                                                                                                cACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGATGAAGTACGACACGGTGATGAacaccacacc                    <  2:113344/147‑1 (MQ=255)
                                                                                                                 aCGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGATGAAGTACGACACGGTGATGAACACCACACCg                   >  1:166306/1‑147 (MQ=255)
                                                                                                                                 aGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGATGAAGTACGACACGGTGATGAACACCACACCGCCAATCAGTGCGGTGa   <  6:162998/147‑1 (MQ=255)
                                                                                                                                       aCTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGATGAAGTACGACACGGTGATGAACACCACACCGCCAATCAGTGCGGTGAg  >  1:9931/1‑142 (MQ=255)
                                                                                                                                           |                                                                                                                                         
GCGAATGCAGGCGGGCGAAGACGCTGGCCGGGATCACGTTGTCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGATGAAGTACGACACGGTGATGAACACCACACCGCCAATCAGTGCGGTGAG  >  NC_002947/1406110‑1406386

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑