Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,590,3340TG22.7% 34.9 / 11.7 22I316L (ATT→CTT) yjiYputative transporter
Reads supporting (aligned to +/- strand):  ref base T (8/9);  new base G (2/3);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACGCGCGTCGGTTTCATTAGCCAGCAGTTTTGGCGTCGTACCGGAAGAGATCAGCGCGTGGAAGCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTGGCAGCAGTGC  >  NC_000913/4590190‑4590472
                                                                                                                                                |                                                                                                                                          
gcactgctgccagtgtggctgatcctcgcaccgcgcgactatctggcaaccttcctgaaaatcggcgttatcgtcggcctggcgctgggtatcgtggtgctgaacccggaactgaaaatgcctgccatgacCCAGTACCGTCAATGTAc                                                                                                                                        <  8:21264‑M1/18‑1 (MQ=255)
   ctgctgccagtgtggctgatcctcgcaccgcgcgactatctggcaaccttcctgaaaatcggcgttatcgtcggcctggcgctgggtatcgtggtgctgaacccggaactgaaaatgcctgccatgacCCAGTACCGTCAATGTACTgg                                                                                                                                     >  7:21264‑M1/129‑149 (MQ=255)
                      cagcagTTTTGGCGTCGTACCGGAAGAGATCAGCGCGTGGAAGCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCa                                                                                                                  <  2:6560/149‑1 (MQ=255)
                                           ggAAGAGATCAGCGCGTGGAAGCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGa                                                                                             >  6:111243/1‑149 (MQ=255)
                                                     aGCGCGTGGAAGCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTc                                                                                                       <  4:85492/129‑1 (MQ=255)
                                                                gCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTcc                                                                                                             <  3:112908/112‑1 (MQ=255)
                                                                gCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTcc                                                                                                             >  4:112908/1‑112 (MQ=255)
                                                                             ccacAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGAtttt                                                           >  1:123822/1‑149 (MQ=255)
                                                                                 aGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGc                                                                  <  3:103835/138‑1 (MQ=255)
                                                                                 aGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGc                                                                  >  4:103835/1‑138 (MQ=255)
                                                                                    cGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAg                                                    <  5:111243/149‑1 (MQ=255)
                                                                                                       cGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTcgcgcg                                 <  2:123822/149‑1 (MQ=255)
                                                                                                         ccggaactgaaaatgcctgccatgacCCAGTACATTGACGGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTc                                      <  2:29829‑M1/116‑1 (MQ=255)
                                                                                                         ccggaactgaaaatgcctgccatgacCCAGTACATTGACGGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTc                                      >  1:29829‑M1/27‑142 (MQ=255)
                                                                                                                     atgcctgccatgacCCAGTACATTGACGGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAgc                                                                                     >  2:170770‑M1/15‑83 (MQ=255)
                                                                                                                     atgcctgccatgacCCAGTACATTGACGGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAgc                                                                                     <  1:170770‑M1/69‑1 (MQ=255)
                                                                                                                             catgacCCAGTACATTGACGGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTgg           <  2:125752‑M1/143‑1 (MQ=255)
                                                                                                                             aGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTgg           >  3:47252/1‑149 (MQ=255)
                                                                                                                                    CAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTGGCAGCAGt    <  5:127710‑M1/149‑1 (MQ=255)
                                                                                                                                      GTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTGGCAGCAGTGc  >  6:127710‑M1/1‑149 (MQ=255)
                                                                                                                                          cGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTGGCAGCAGt    <  2:90565/143‑1 (MQ=255)
                                                                                                                                          cGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTGGCAGCAGt    >  1:90565/1‑143 (MQ=255)
                                                                                                                                                |                                                                                                                                          
AACGCGCGTCGGTTTCATTAGCCAGCAGTTTTGGCGTCGTACCGGAAGAGATCAGCGCGTGGAAGCCAGATACCGCACCACAGGCGATGGTGATGAACAGGAACGGGAACAGAGCGCCTTTCCACAGCGGGCCAGTACCGTCAATGTACTGGGTCATGGCAGGCATTTTCAGTTCCGGGTTCAGCACCACGATACCCAGCGCCAGGCCGACGATAACGCCGATTTTCAGGAAGGTTGCCAGATAGTCGCGCGGTGCGAGGATCAGCCACACTGGCAGCAGTGC  >  NC_000913/4590190‑4590472

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