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

AAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGC‑AAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCC  >  NC_000913/1157768‑1158030
                                                                                                                                 |                                                                                                                                       
aaataagcgaattgtAGATGCCCTATACACGGCGATGCTCTCCCCCCGTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGTATGCATTTGCTAACCTGCAAAA‑AGGTCGGTaa                                                                                                                           <  2:241716/129‑1 (MQ=255)
                 aTGCCCTGTACACGGCCAGGCTCTACCCCCTTGCCACGCGTGAGAACGTAAAAACAGCACCCATACTCACCAGCACTCTCAATTATGTATAAGAATGCATTTCCTAACCTGCAAAA‑AGGTCCGTAAATCGCTAAtgctgcagc                                                                                                          >  2:254246/1‑140 (MQ=255)
                       tttACACGGCGAGGCTCTCCCCCCTTTCCACTCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCg                                                                                                    <  1:78181/141‑1 (MQ=255)
                                 gAGGGTCTCCCCCCGTGCCACGCGTGAGAATGTAACAAAAGCACCCATACTCAAGAGCACTCTCAAGTATGTTTAAGAATGCATTTGCTGACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTAt                                                                                          <  2:57874/143‑1 (MQ=255)
                                                          tataaCGTCAAGAAAGCACCCGTTCTCCGGAGCACTCTCAATTATGCTTAAGAATGCATTTGCTAACCTGC‑AAATAGGTCGGTAAAGCGCTGATGCTGCCGGTATCCGTATTGTCTATCCCAGGTATTCTGCTGGGCGTCGGt                                                                 <  2:198222/140‑1 (MQ=255)
                                                               cgtaacagtaGCCCCCATCCTTCGGAGCACTCCCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCcgc                                                            <  2:149243/134‑1 (MQ=255)
                                                                                                cAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTCTCCGTACTGCCTATCGCATGTATTCTGCTGGGAAGCGAGAACGCGAGTTTCAGCTGGCTGCCCGCCGTTGTATCGCa                           >  1:295206/1‑143 (MQ=255)
                                                                                                          tAAGAATGCATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCTGTATCCGTACTGCCTATCGCAGGTATTCTTCTGTGCGTCGGTTCCACGAATTTCAGCTAGCTGACCGCCGTTGTATCGCATGTTATGTCa                 >  1:190141/1‑143 (MQ=255)
                                                                                                                  cATTTGCTAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGACGGTATCCGTACTGCCTATCGCAGGTATTCTGCTCGGCGTCGGTTCCGCGAATTTCAGCTCGCTGCCCGACGTGGTATCGCAGGTTCTGGCGGAAGCAgg         >  2:8826/1‑143 (MQ=255)
                                                                                                                         tAACCTGCAAAA‑AGGTCGGTACATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCCTTCTATCTCATGTTCTGGCAGAACCAGGCCGTTcc  >  2:77488/1‑143 (MQ=255)
                                                                                                                         tAACCTGCAAAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATAGCAAAAGCAGGCggaagc  >  1:292791/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                       
AAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGC‑AAA‑AGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCC  >  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: ‑