Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,349,621 | 0 | A | C | 31.8% | 24.6 / 19.4 | 22 | intergenic (‑149/‑66) | inaA/yfaH | putative lipopolysaccharide kinase InaA/putative uncharacterized protein YfaH |
Reads supporting (aligned to +/- strand): ref base A (15/0); new base C (0/7); total (15/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.86e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCCACTACTAGCTTGACACCGCTTTTACCCTTCATTTG‑CAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATACAGCGAT > NC_000913/2349567‑2349688 | tcccACTCCTAGCTTG‑CCCCCCTTTTTCCCTTTAATTGCCGGCGCAAAGATCCCCCCAGccactacca < 1:3849549/66‑1 (MQ=255) tcccACTACTATCTTCACCCCTTTTTTCCCCTTCTTTGC‑CAGCGAAAATATCCCCCCAGccactacca < 1:2841162/66‑1 (MQ=255) tcccACTACTACCTTGCCCCCGTTTTTCCCCTTCTTTTG‑CAGCGAATAGATCCCCCCAGccactacca < 1:1744622/66‑1 (MQ=255) ccACTACTAGCTTGACCCCGCTTTTCCCCTTCATGTG‑CAGCGCATAGATCCCCCCAGccactaccact < 1:332448/68‑1 (MQ=255) ccACTACGACCTTGCCCCCGTTTTTCCCTTTCTTTTG‑AAGAGCATAGCTCCCCCCAGccactaccact < 1:3226956/68‑1 (MQ=255) cACTACTAGCTTGACACCGCTTTTACCCTTCATTTG‑CAGCGCATAGATTCCACCAGccactaccacta > 1:2442539/1‑68 (MQ=255) actaGCTTGACCCCGTTTTTCCCCTTCATTTG‑CAGCGAATAGATCCCCCCAGCCACTACCACTAATAg < 1:3562889/68‑1 (MQ=255) taGCTTGACACCGCTTTTACCCTTCATTTG‑CAGCGCATAGATTCCACCAGCCACTACCACTAATAGTg > 1:1533094/1‑68 (MQ=255) gCTTGACACCGCTTTTACCCTTCATTTG‑CAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgcg > 1:2610229/1‑68 (MQ=255) cTTGCCCCCGTTTTTCCCCTTCATTGC‑CAGCGAATAATTCCCCCCAGCCACTACCACTAATAGTgcgc < 1:809743/68‑1 (MQ=255) cTTGACACCGCCTTTACCNTTCATTTG‑CAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgcgc > 1:3089380/1‑68 (MQ=255) acacCGCTTTTACCCTTCATTTG‑CAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCcag > 1:1964491/1‑68 (MQ=255) cacCGCTTTTACCCTTCATTTG‑CAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCcagc > 1:1988200/1‑68 (MQ=255) ctcTTTTCCCCTTCATTTG‑CAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGCAGa < 1:166488/66‑1 (MQ=255) cgtttttCCCCTTCATTGC‑CCGCGCATAGCTCCCCCCAGCCACTACCACTAATAGTGCGCCCAGCAGa < 1:2070423/65‑1 (MQ=255) ttttACCCTTCATTTG‑CAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATaa > 1:1477260/1‑68 (MQ=255) ccTTCATTTG‑CAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAAtgtg > 1:490351/1‑68 (MQ=255) ccTTCATTTG‑CAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAAtgtg > 1:2307421/1‑68 (MQ=255) aGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAg > 1:1202483/1‑68 (MQ=255) gcATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCg < 1:3615516/68‑1 (MQ=255) gATTCCACCAGCCACTACCACTAATAGTGCGCCCCGCAGATAAAATGTGAATTTAATCAGGCGGATAc > 1:1984011/1‑68 (MQ=255) gATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATAc > 1:2302751/1‑68 (MQ=255) ttCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATACAg > 1:881160/1‑68 (MQ=255) ttCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATACAg > 1:2931589/1‑68 (MQ=255) ttCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATACAg > 1:2592909/1‑68 (MQ=255) accaGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATACAGCGAt > 1:2251073/1‑68 (MQ=255) | TGCCACTACTAGCTTGACACCGCTTTTACCCTTCATTTG‑CAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATACAGCGAT > NC_000913/2349567‑2349688 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |