Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,973,535 | 0 | T | C | 54.5% | 4.7 / 12.7 | 11 | D597G (GAT→GGT) | cheA | fused chemotactic sensory histidine kinase in two‑component regulatory system with CheB and CheY: sensory histidine kinase/signal sensing protein |
Reads supporting (aligned to +/- strand): ref base T (2/3); new base C (0/6); total (2/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.82e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03 |
TTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGC‑AGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAA‑CCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCC > NC_000913/1973387‑1973669 | ttCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGC‑AGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAA‑CCTGGTGTTGACCAATTAATTGATc < 2:315805/150‑1 (MQ=255) cAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGC‑AGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAA‑CCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAgcgg < 2:245470/151‑1 (MQ=255) cAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGC‑AGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAA‑CCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAgcg < 2:383884/150‑1 (MQ=255) gTGCCACGCTGCCGTCGCCAAGAATGGTCGC‑AGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAA‑CCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATca > 1:369860/1‑150 (MQ=255) cgcgAAGATGGGTAGC‑AGCAAAAATGCGGGGTCTTTTGGGATA‑TTTCATTCAAAGTTTTTAACCCCAACCCGGGTGTGAAAGAATAATTTGCCCCCCCAGAAGGGGGTGGGGGGGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGt < 1:350488/151‑1 (MQ=255) gcaaaGAATGG‑CGCAAGCAAAAACGCGGGGGAGTTGGCAATAGTTTCTTTCAAGGTTTTAACCCACAC‑CCGGTTGTTGCCAAATTAATGGACCCCCCAGCAGGGGGTAGGGGGGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGt < 1:218927/147‑1 (MQ=255) atatttttttaaaaGTTTTTAAACCACAC‑CCGGTTGTTAACCAATTAATACACCCCCCAACAAGGCGTGGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCcac < 2:215543/139‑1 (MQ=255) aTAGTTACTTTCAAGGTTTTTAACCACAA‑CCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCca > 2:131455/1‑150 (MQ=255) ccACAC‑CCGTGTGTTGACAAATTAATTAACCCCCCAGCAGGGGGTGGGGGCGCCCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTc < 2:76164/151‑1 (MQ=255) aa‑CTTGGTGTTAACAAATAAATTCACCCACCAGCAGGGCGGAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCAccc < 1:14534/151‑1 (MQ=255) ttaacaaataaatgGACCCCCCGGAAAGGGGGAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTcc < 2:283816/145‑1 (MQ=255) | TTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGC‑AGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAA‑CCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCC > NC_000913/1973387‑1973669 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |