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.A88.9% 20.0 / ‑3.0 9coding (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/2);  total (6/3)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

ATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCC  >  NC_000913/1157785‑1158030
                                                                                                                |                                                                                                                                       
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)
                                                                               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)
                                                                                                                |                                                                                                                                       
ATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCC  >  NC_000913/1157785‑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: