Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 2,866,401 | 0 | G | A | 67.5% | 44.2 / 35.2 | 40 | Q337* (CAA→TAA) | nlpD | predicted outer membrane lipoprotein |
Reads supporting (aligned to +/- strand): ref base G (0/13); new base A (14/13); total (14/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.09e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGAT > W3110S.gb/2866354‑2866467 | ggTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt < 1:776121/71‑1 (MQ=255) ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATtg > 1:1097116/1‑45 (MQ=255) ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATtg > 1:1112268/1‑45 (MQ=255) ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATtg > 1:3552129/1‑45 (MQ=255) ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATtg > 1:3734980/1‑45 (MQ=255) ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATtg > 1:1583806/1‑45 (MQ=255) ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt < 1:850629/51‑1 (MQ=255) ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt < 1:3110721/51‑1 (MQ=255) ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt < 1:1992884/51‑1 (MQ=255) ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt < 1:1729207/51‑1 (MQ=255) ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt < 1:1725412/51‑1 (MQ=255) ttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCa > 1:1388661/1‑40 (MQ=255) ttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt > 1:1478827/1‑49 (MQ=255) ttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt > 1:209454/1‑49 (MQ=255) ttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt > 1:2201844/1‑49 (MQ=255) ttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTa > 1:1390484/1‑51 (MQ=255) ttttGCCCCGCCTAAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt > 1:1632646/1‑49 (MQ=255) tAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt < 1:3062960/36‑1 (MQ=255) aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt < 1:4085314/35‑1 (MQ=255) aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt < 1:3574827/35‑1 (MQ=255) aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt < 1:2491828/35‑1 (MQ=255) aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt < 1:3749229/35‑1 (MQ=255) aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt < 1:1752391/35‑1 (MQ=255) aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt < 1:3337507/35‑1 (MQ=255) aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt < 1:419988/35‑1 (MQ=255) aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt < 1:689165/35‑1 (MQ=255) aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt < 1:1956167/35‑1 (MQ=255) aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt < 1:1569855/35‑1 (MQ=255) aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt < 1:141030/35‑1 (MQ=255) aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt < 1:136057/35‑1 (MQ=255) aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt < 1:106816/35‑1 (MQ=255) aaCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATgg > 1:2660763/1‑40 (MQ=255) aaCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt > 1:3293165/1‑71 (MQ=255) aaCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt > 1:1282330/1‑71 (MQ=255) aCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGttt < 1:1086840/71‑1 (MQ=255) aCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGttt < 1:2978763/71‑1 (MQ=255) aCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGttt < 1:2826075/71‑1 (MQ=255) aCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGttt < 1:268498/71‑1 (MQ=255) aCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGttt < 1:1355826/71‑1 (MQ=255) ttcttGTTATTCCCGGACCAGCATTGTGTCGTTATg < 1:1526964/36‑1 (MQ=255) gttattCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat < 1:3303581/70‑1 (MQ=255) gttattCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat < 1:2728014/70‑1 (MQ=255) gttattCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat < 1:4108541/70‑1 (MQ=255) gttattCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat < 1:2394376/70‑1 (MQ=255) gttattCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat < 1:2225912/70‑1 (MQ=255) | GGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGAT > W3110S.gb/2866354‑2866467 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |