Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 2,797,548 | 1 | . | C | 71.0% | 99.7 / 34.2 | 38 | intergenic (‑397/‑272) | stpA/ygaW | DNA binding protein, nucleoid‑associated/predicted inner membrane protein |
Reads supporting (aligned to +/- strand): ref base . (11/0); new base C (27/0); total (38/0) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACACCGAATTCTCTACGCTTATAGCGTTACCTCA‑CCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAATT > W3110S.gb/2797515‑2797595 | acacCGAATTCTCTACGCTTATAGCGTTACCTCA‑CCCCCAAACTTTCATTAGTGa > 1:2141206/1‑55 (MQ=255) acacCGAATTCTCTACGCTTATAGCGTTACCTCA‑CCCCCAAACTTTCATTAGTGATATTAAAg > 1:672038/1‑63 (MQ=255) acacCGAATTCTCTACGCTTATAGCGTTACCTCA‑CCCCCAAACTTTCATTAGTGATATTAAAg > 1:661098/1‑63 (MQ=255) acacCGAATTCTCTACGCTTATAGCGTTACCTCA‑CCCCCAAACTTTCATTAGTGATATTAAAg > 1:2140752/1‑63 (MQ=255) acacCGAATTCTCTACGCTTATAGCGTTACCTCA‑CCCCCAAACTTTCATTAGTGATATTAAAg > 1:339666/1‑63 (MQ=255) aCGCTTATAGCGTTACCTCA‑CCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCat > 1:222417/1‑59 (MQ=255) aCGCTTATAGCGTTACCTCA‑CCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCat > 1:133751/1‑59 (MQ=255) aCGCTTATAGCGTTACCTCA‑CCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCat > 1:375220/1‑59 (MQ=255) aCGCTTATAGCGTTACCTCA‑CCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCat > 1:367982/1‑59 (MQ=255) aCGCTTATAGCGTTACCTCA‑CCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCat > 1:1748428/1‑59 (MQ=255) aCGCTAATAGCGTTACCTCA‑CCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCat > 1:1350127/1‑59 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:360121/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:2236774/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:2278141/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:229111/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:282304/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:340585/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:767673/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:442843/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:550513/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:751549/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:1077269/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:2216814/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:2199075/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:2067036/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:2065135/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:1920046/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:1908854/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:17714/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:1712135/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:162155/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:1614745/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:1596039/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:1577109/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:1523282/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:1473142/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:1338769/1‑65 (MQ=255) cTTATAGCGTTACCTCACCCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAAtt > 1:1138089/1‑65 (MQ=255) | ACACCGAATTCTCTACGCTTATAGCGTTACCTCA‑CCCCCAAACTTTCATTAGTGATATTAAAGGTAAGGCCATATCTAATT > W3110S.gb/2797515‑2797595 |
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: ‑ |