Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 4,297,115 | 0 | G | A | 21.9% | 56.4 / 12.7 | 33 | P168S (CCC→TCC) | alsK | D‑allose kinase |
Reads supporting (aligned to +/- strand): ref base G (12/13); new base A (2/5); total (14/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.26e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTT > CP009273/4296974‑4297261 | tttCCGCATGGACGAAAAGATCGCCCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGa < 2:110099/149‑1 (MQ=255) gatctgCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGa < 2:176128/144‑1 (MQ=255) ttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAc < 2:570548/149‑1 (MQ=255) cGCATGGACGAAAAGATCGGGCAAGGGGTAATTTCGGAGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCa < 2:231295/149‑1 (MQ=255) cATGGACGAAACGATCGCGCAATGGCTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAc < 1:353237/144‑1 (MQ=255) cATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCAATTCCAGAGCACTTGGTTTCCAGGCACCCAGGAATGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAc > 1:160627/1‑144 (MQ=255) tgggACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCGGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGtt < 2:496174/147‑1 (MQ=255) tGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAccc < 1:11943/144‑1 (MQ=255) ggccGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCg < 2:256157/146‑1 (MQ=255) gAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg > 2:328837/1‑149 (MQ=255) tAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACTCCGTGTGCAccaca > 2:720698/1‑146 (MQ=255) tCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCAcc > 2:346650/1‑148 (MQ=255) ggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCAc > 1:59984/1‑144 (MQ=255) cTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACgg > 1:492960/1‑144 (MQ=255) gttgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGCTATGACCCAGTTCGCCTGCCCCACCGTGTGCACCCGTCCgcggcg > 1:412111/1‑144 (MQ=255) cAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGCCATATCTCCCAGGGGCATATGACCCAGTTCGCCTACCACACCGTGTGCACCCGTCCACGGCGCACCGTTCGTc > 1:34393/1‑144 (MQ=255) tCTTAGCCCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGACATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGa < 2:50821/149‑1 (MQ=255) ttAGCGCCATTCCAGAGCAATTGGTTCCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAc < 2:14117/149‑1 (MQ=255) gcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCAt < 2:564329/149‑1 (MQ=255) agagCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTAc > 2:312200/1‑149 (MQ=255) gagCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATAGCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTAcc < 2:338667/149‑1 (MQ=255) cagttGGTTTCGCGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCCTGTCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTa < 1:183819/141‑1 (MQ=255) cAATTGGTTTCCAGGCACCCAGGATGGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTa < 1:604152/144‑1 (MQ=255) ttGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCAGCGTTCATCCACACTGCGAACCCCATCCCCGTACCg > 1:399181/1‑144 (MQ=255) attCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGACATCTCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgcc > 2:659503/2‑149 (MQ=255) aGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATACCTCCCAGGGGGATGTGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACTGCGCACCGCTCCTCCACACGGCGGACCCCACCCCCGTACCCAGCGAGgcg > 1:392647/1‑143 (MQ=255) cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGGGATAGGCCGCCAGAACCAgt > 2:151919/1‑149 (MQ=255) tgtTGGGGCATATCCCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTc < 2:84858/149‑1 (MQ=255) gggTCATATCTCCCAGGGCGATATGACCCAGTTCGCCTGCCACACCGGGTGCACCCGTCCACGGCGCACCGGTCATCCACCCTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGtttt < 1:401908/144‑1 (MQ=255) tCATATCGCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTcta < 1:295516/144‑1 (MQ=255) gcctacctcCCAGGGGGATATGCCCCAGTTCGCCGGCCACACCGTGTGCCCCCGTCCACGGCGCACCGTTCAGCCACACTGCGAACCCCACCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTcta < 1:11429/138‑1 (MQ=255) atatCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCATGTGCACCCATCCCCGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTc > 2:177819/1‑149 (MQ=255) ccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGgag < 2:66523/149‑1 (MQ=255) ggggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGtt > 2:595779/1‑149 (MQ=255) ggggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGtt > 2:94783/1‑149 (MQ=255) | TTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTT > CP009273/4296974‑4297261 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |