Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,406,2370GA27.6% 44.7 / 11.4 29F288F (TTC→TTT) puuPputrescine permease
Reads supporting (aligned to +/- strand):  ref base G (9/12);  new base A (8/0);  total (17/12)
Fisher's exact test for biased strand distribution p-value = 8.92e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCACTGGCGTCTTGTAGCGCGAATGCAGGCGGGCGAAGACGCTGGCCGGGATCACGTTGTCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGATGAAGTACGACACGGTGATGAACACCACACCG  >  NC_002947/1406092‑1406369
                                                                                                                                                 |                                                                                                                                    
gCACTGGCGTCTTGTAGCGCGAATGCAGGCGGGCGAAGACGCTGGCCGGGATCACGTTGTCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGaa                                                                                                                                    <  8:10612/148‑1 (MQ=255)
    tGGCGTCTTGTAGCGCGAATGCAGGCGGGCGAAGACGCTGGCCGGGATCACGTTGTCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCAAGGCGATGGTGTTGATAACGGTGCACGCGATAAAAAt                                                                                                                                 >  4:105172/1‑147 (MQ=255)
       cGTCTTGTAGCGCGAATGCAGGCGGGCGAAGACGCTGGCCGGGAGCACGTTGGCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGAc                                                                                                                             <  5:115788/148‑1 (MQ=255)
                  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)
                                                                                                                                                 |                                                                                                                                    
GCACTGGCGTCTTGTAGCGCGAATGCAGGCGGGCGAAGACGCTGGCCGGGATCACGTTGTCGCGGCCCATTACATACAGCAGGCGTGACACGCTGGTTTGCGAGGCCAGGCCCGAGGCGATGGTGTTGATCACGGTGCAGGCGATGAAGATCGACTGGAAGAGCTTGCCGCCAACGTACAGGGCAATTTCCGGCAGGGCTGCTTCCTGGTCGTGGAAGCGGGCCATGGTCGGGAAGTAGGCCTGGATGAAGTACGACACGGTGATGAACACCACACCG  >  NC_002947/1406092‑1406369

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: ‑