Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,852,477 | 0 | G | A | 52.6% | 6.3 / 18.4 | 19 | E676K (GAA→AAA) | mutS | DNA mismatch repair protein MutS |
Reads supporting (aligned to +/- strand): ref base G (6/3); new base A (0/10); total (6/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.10e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACC‑TTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTA > NZ_CP009273/2852396‑2852556 | ccTATCGATCGCATCTTTCCCCGCGTAGGCGCGGCAGATGCCCTGGCGTCCGGGCGCTCAACC‑TTTATGGTGGAGATGACTGAAACCGcc < 1:391449/90‑1 (MQ=255) tttCCCCGCGTGGGCGCGGCAGATGCCCTGGCGTCCGGGCGCTCAACC‑TTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAAc < 1:261582/90‑1 (MQ=255) tttCCCCGCGGGGGCGGGGCAGATGCCCGGGCCCCGGGGCCCCCACCT‑TTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAAc < 2:206572/90‑1 (MQ=255) tccccGGGGAGGCGCGGCAGATGCCCGGGCG‑CCGGGGCCCCAACCTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGc < 2:700874/88‑1 (MQ=255) tccccGCGTAGGCGCGGCAGATGCCCTGGCGTCCGGGCGCTCACCC‑TTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGc < 1:460334/88‑1 (MQ=255) tccccGCGAGGGCGCGGAAGATGCCCGGGCGCCCGGGCCCTCAACC‑TTTATGGGGGAGATGACTAAACCCGCCAATATTTTACATAACGc < 2:130178/88‑1 (MQ=255) tACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACC‑TTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGc > 1:213970/1‑90 (MQ=255) cgcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACC‑TTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCAcc > 2:676073/1‑90 (MQ=255) cgcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACC‑TTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCAcc > 2:781789/1‑90 (MQ=255) ggCGCGGCAGTTGCCCGGGGGTCCGGCCGCTCACCC‑TTTATGGGGGGGTTGACTAAAACCGCCAATATTTTACATAACGCCACCGAATAc < 2:673551/90‑1 (MQ=255) ggCGCGGAAGATGCCCTGGCGCCGGGGCCCCCAACT‑TTTAGGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATAc < 2:605487/90‑1 (MQ=255) ggCGCGGAAGATGCCCGGGCGTCGGGGCCCCAAACT‑TTTATGGGGGAGAGGACGAAACCCGCCAATATTTTACATAACGCCACCGAATAc < 2:730874/90‑1 (MQ=255) cgGCAGAGGCCCGGGCGTCGGGGCCCTAACCT‑TTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTc < 2:264750/90‑1 (MQ=255) cgGAAGATGCCCGGGCGTCGGGGCGC‑CAACCTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTc < 1:883177/90‑1 (MQ=255) gggCCCTCACCC‑TTTGTGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATc < 2:1109997/90‑1 (MQ=255) ttttATGGTGGAGATGACAAAACCCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAAc < 2:802213/89‑1 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 1:691450/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 2:486407/1‑90 (MQ=255) ggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTa > 1:802213/1‑90 (MQ=255) | CCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACC‑TTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTA > NZ_CP009273/2852396‑2852556 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |