Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,011,323:1 | +T | 100% | intergenic (‑258/+192) | PP_0871 ← / ← SBRG_Ecoli_galM | glycine betaine/carnitine/choline ABC transporter permease/SBRG_Ecoli_galM |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,011,323 | 1 | . | T | 82.4% | 35.8 / 5.9 | 17 | intergenic (‑258/+192) | PP_0871/SBRG_Ecoli_galM | glycine betaine/carnitine/choline ABC transporter permease/SBRG_Ecoli_galM |
Reads supporting (aligned to +/- strand): ref base . (0/3); new base T (8/6); total (8/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.06e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
AGTCCCCAACCGAAGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATG‑CTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGT > NC_002947/1011192‑1011427 | aGTCCCCAACCGAAGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGtctc < 2:267882‑M1/136‑5 (MQ=255) aaGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCtttt < 2:198993/136‑1 (MQ=255) gTTGCTGGGGAGTGGTGGGCTCGTTGGGCGTATAGGTTCCCATCTGAGGGGGCATTTGGCCTGTGTTAATCGGTGGGTTGCGGGTTTTTGGGGGAGATAAACAGCGTCTTGATGTCTCAAATCATCTTTTAATGat < 2:424161/136‑1 (MQ=255) cAGGGTAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGTGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGataata < 2:469066/136‑1 (MQ=255) aCGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATtat > 1:293187/1‑136 (MQ=255) aCGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATtat > 1:327455/1‑136 (MQ=255) gTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTatat > 1:377083/1‑136 (MQ=255) cGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGcgacga > 2:312674/1‑136 (MQ=255) ggTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAAc < 2:266072/136‑1 (MQ=255) gcggcATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGtgctgcaggcgttgaaataactgcgttttggttatttgtcgtgataaagcgattgcctgttgg < 2:181618‑M1/129‑64 (MQ=255) gcggcATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGtgctgcaggcgttgaaataactgcgttttggttatttgtcgtgataaagcgattgcctgttgg > 1:181618‑M1/1‑66 (MQ=255) gcggcATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGc > 1:513879/1‑136 (MQ=255) tttGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAAcc < 1:312674/136‑1 (MQ=255) tAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCTATGCTTACGCATAACCATAGCGAAAATa > 2:407056/1‑136 (MQ=255) gCTGGGTTGCGGGTTTTTGGGCCATAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGgc < 1:9016/136‑1 (MQ=255) gTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAgt > 1:364492/1‑136 (MQ=255) ttttttgGGCGAGAAAAACAGCGACTTGATGtgctgcaggcgttgaaataactgcgttttggttatttgtcgtgataaagcgattgcctgttgggtttgttttgtttgtgccggcctcttcgcgggtaaactcgct < 1:303928‑M1/135‑106 (MQ=255) gTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTAAATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCCAAAATAGTGGCGCAGTGTAAGgt > 2:108054/1‑136 (MQ=255) | AGTCCCCAACCGAAGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATG‑CTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGT > NC_002947/1011192‑1011427 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |