Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947923,6490TG58.8% ‑1.1 / 12.3 18Q516P (CAG→CCG) PP_0804protein secretion ABC transporter permease/ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base G (10/0);  total (12/6)
Fisher's exact test for biased strand distribution p-value = 5.66e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCGATGCCGTCCAGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCGGTACTCGCCGTGAATGGCCGGCAGGTGCACGCGGGCTTCGCCTTCGGGGTGATCGACCGG  >  NC_002947/923530‑923782
                                                                                                                       |                                                                                                                                     
tGGCGATGCCGTCCAGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCGCGGCGGGGGATACAGCCGGcc                                                                                                                       >  2:312979/1‑136 (MQ=255)
tGGCGATGCCGTCCAGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCGCCGCGGGTTATACCGCCGGcg                                                                                                                       >  1:395988/1‑135 (MQ=255)
   cGATGCCGTCCAGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATCCAGCCGGCCAAt                                                                                                                    <  2:287112/136‑1 (MQ=255)
    gATGCCGTCCAGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATACCGCCGGCCAATg                                                                                                                   >  2:331501/1‑136 (MQ=255)
              aGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCGGGGTATACAGCCGGCCCATGTTGAGCACTg                                                                                                         >  2:339673/1‑136 (MQ=255)
                       tcGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCGCCGCGGGGTATGCAGCCCGCCAAAGTTGAGCACAGGGGGGTTGt                                                                                                >  1:121486/1‑136 (MQ=255)
                       tcGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGGGGGGTATGCAGCCCGCCAAACTTCAGCACAGGCGGGGAGt                                                                                                >  1:349574/1‑136 (MQ=255)
                          cccGGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCtt                                                                                             <  2:210135/136‑1 (MQ=255)
                                          tGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCCGCGGGGTATAAACCCGGCCCATGTTAAGCACAGGCGGGTAGTCTTAGCTGTATCGGAAGtt                                                                             >  2:381143/1‑136 (MQ=255)
                                            gCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCCGCAGGATATAAACCCCCCCAAAGTTAAACACAGGCCGGGAGTCTTAGCAGTATCGGAAAtttt                                                                           >  1:59869/1‑134 (MQ=255)
                                            gCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCCCCGGGGTATACAGCCCGCCCAATATGAGCACAGGCGGGGAGTAATACGTGTTGAGGAAGTTgg                                                                           >  2:374674/1‑136 (MQ=255)
                                             ccccACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCCGCGGGATATACACCCCGCCAAAGTTAAGCACAGGCGGGGAGTATTCGCGGTAGCGGGAGtttgg                                                                          >  2:102270/1‑133 (MQ=255)
                                                         ccTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCg                                                              <  1:310721/136‑1 (MQ=255)
                                                         ccTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCg                                                              <  1:312979/136‑1 (MQ=255)
                                                              agcaAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCGCGGCGGGGATTACAGCCCGCCCATGTTCAGCACTGGCGGGGATTAATAGCTGTAGCGGGAGTTTGACAGGCGCCGGCGGGAca                                                         >  2:286422/1‑135 (MQ=255)
                                                                                                            gcTCGCCCCGGGGGGTATACAGCCGGCCCATGTTGAGCACTGGGGGGGTGTGTTGGGTGTAGCGGGAGTTGGGCTGGCGCAGGCGGTTCTCGCCGTGAATGGCCGGCAGGGGGACGCGGGCGTCGCCGTCGGGGTg           >  2:372789/1‑136 (MQ=255)
                                                                                                                    ggCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCGGTACTCGCCGTGAATGGCCGGCAGGTGCACGCGGGCTTCGCCTTCGGGGTGATCGACCg   <  2:138321/136‑1 (MQ=255)
                                                                                                                    ggCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCGGTACTCGCCGTGAATGGCCGGCAGGTGCACGCGGGCTTCGCCTTCGGGGTGATCGACCg   <  2:135600/136‑1 (MQ=255)
                                                                                                                     ggtGGATATCCAGCCGGCCAATGTTCAGCACTGGCTGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCGGTACTCGCCGTGAATGGCCGGCAGGTGCACGCGGGCTTCGCCTTCGGGGTGATCGACCgg  <  1:58374/134‑1 (MQ=255)
                                                                                                                       |                                                                                                                                     
TGGCGATGCCGTCCAGGCTGATCTCGCCCTGCACCAGATCCATGGCCCCACCCAGCGCCTGCAGCAAGGTTGATTTGCCGGCGCCGTTGCGGCCCAGCACGGCAATGCGCTCGCCCGGCTGGATATCCAGCCGGCCAATGTTCAGCACTGGCGGGTAGTCTTCGCTGTAGCGGAAGTTGGCCTGGCGCAGGCGGTACTCGCCGTGAATGGCCGGCAGGTGCACGCGGGCTTCGCCTTCGGGGTGATCGACCGG  >  NC_002947/923530‑923782

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: