Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,157,900 +A coding (32/1434 nt) ptsG → fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,157,8961.A90.9% 26.9 / ‑3.1 11coding (28/1434 nt)ptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (6/4);  total (6/5)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

AAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCC  >  NC_000913/1157768‑1158030
                                                                                                                                 |                                                                                                                                       
aaataagcgaattgtAGATGCCCTATACACGGCGATGCTCTCCCCCCGTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGTATGCATTTGCTAACCTGCAAAAAGGTCGGTaa                                                                                                                           <  2:241716/129‑1 (MQ=255)
                 aTGCCCTGTACACGGCCAGGCTCTACCCCCTTGCCACGCGTGAGAACGTAAAAACAGCACCCATACTCACCAGCACTCTCAATTATGTATAAGAATGCATTTCCTAACCTGCAAAAAGGTCCGTAAATCGCTAAtgctgcagc                                                                                                          >  2:254246/1‑140 (MQ=255)
                       tttACACGGCGAGGCTCTCCCCCCTTTCCACTCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCg                                                                                                    <  1:78181/141‑1 (MQ=255)
                                 gAGGGTCTCCCCCCGTGCCACGCGTGAGAATGTAACAAAAGCACCCATACTCAAGAGCACTCTCAAGTATGTTTAAGAATGCATTTGCTGACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTAt                                                                                          <  2:57874/143‑1 (MQ=255)
                                                          tataaCGTCAAGAAAGCACCCGTTCTCCGGAGCACTCTCAATTATGCTTAAGAATGCATTTGCTAACCTGCAAATAGGTCGGTAAAGCGCTGATGCTGCCGGTATCCGTATTGTCTATCCCAGGTATTCTGCTGGGCGTCGGt                                                                 <  2:198222/140‑1 (MQ=255)
                                                               cgtaacagtaGCCCCCATCCTTCGGAGCACTCCCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCcgc                                                            <  2:149243/134‑1 (MQ=255)
                                                                                                cAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTCTCCGTACTGCCTATCGCATGTATTCTGCTGGGAAGCGAGAACGCGAGTTTCAGCTGGCTGCCCGCCGTTGTATCGCa                           >  1:295206/1‑143 (MQ=255)
                                                                                                          tAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCTGTATCCGTACTGCCTATCGCAGGTATTCTTCTGTGCGTCGGTTCCACGAATTTCAGCTAGCTGACCGCCGTTGTATCGCATGTTATGTCa                 >  1:190141/1‑143 (MQ=255)
                                                                                                                  cATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGACGGTATCCGTACTGCCTATCGCAGGTATTCTGCTCGGCGTCGGTTCCGCGAATTTCAGCTCGCTGCCCGACGTGGTATCGCAGGTTCTGGCGGAAGCAgg         >  2:8826/1‑143 (MQ=255)
                                                                                                                         tAACCTGCAAAAAGGTCGGTACATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCCTTCTATCTCATGTTCTGGCAGAACCAGGCCGTTcc  >  2:77488/1‑143 (MQ=255)
                                                                                                                         tAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATAGCAAAAGCAGGCggaagc  >  1:292791/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                       
AAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCC  >  NC_000913/1157768‑1158030

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: