Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,020,300 | 0 | A | C | 53.8% | 9.3 / 16.6 | 13 | E391D (GAA→GAC) | guaD | guanine deaminase |
Reads supporting (aligned to +/- strand): ref base A (6/0); new base C (0/7); total (6/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.83e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTG > NZ_CP009273/3020211‑3020381 | gCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGaa > 1:568151/1‑90 (MQ=255) gCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGaa > 2:536612/1‑90 (MQ=255) aTCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAca > 1:627230/1‑90 (MQ=255) tGGGCCTTGACGTTTTGATTGCCACTTTTTTCCCTGGAAAAGAGGCTGTTTTGGTGGTGAGGAACCCCACCGCCACTCCGCTACagcagc < 2:660243/90‑1 (MQ=255) gggCCTTGACGTTTTGATTGGCATTTTTTTCCCGGGAAAAGAGGTTGTTTTGGTGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCt < 2:208393/90‑1 (MQ=255) cTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTgcgc > 1:660243/1‑90 (MQ=255) ttGACGTTTTGTTGGGCATTTTTTTCCCTGGAAAAGGGGCTGTTTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCt < 2:58608/90‑1 (MQ=255) ttttGATTGGCACTTTTTTCCCTGGAAAAGAGGCTGTTTTCGGGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACa < 2:204675/89‑1 (MQ=255) aTTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACa > 1:494646/1‑90 (MQ=255) tGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGAACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGt < 1:509995/90‑1 (MQ=255) ttttACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAg > 1:242056/1‑90 (MQ=255) tttCCCTGGAAAAGAGGTTGTTTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGt < 2:251865/90‑1 (MQ=255) ggCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACaaa > 2:119942/1‑90 (MQ=255) gtggtgATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACg > 2:455207/1‑90 (MQ=255) gtgATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTg < 1:119942/90‑1 (MQ=255) | GCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTG > NZ_CP009273/3020211‑3020381 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |