Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 933,585 | 0 | C | A | 72.7% | 47.8 / 26.5 | 33 | A121E (GCG→GAG) | ftsK | cell division DNA translocase FtsK |
Reads supporting (aligned to +/- strand): ref base C (5/4); new base A (17/7); total (22/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.38e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTTGCCTGGCGTCATCAGTCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGCGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATT > NC_000913/933508‑933661 | gTTTGCCTGGCGTCATCAGTCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGCGCt < 2:178134/81‑1 (MQ=255) tGCCTGGCGTCATCAGTCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGCGCTCa < 1:178135/80‑1 (MQ=255) ccTGGCGTCATCAGTCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGAGCTCATcc > 1:178162/1‑81 (MQ=255) cTGGCGTCATCAGTCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGAGCTCATCCt > 2:178167/1‑81 (MQ=255) ggCGTCATCAGTCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGAGCTCATCCTTa > 2:178170/1‑81 (MQ=255) gTCATCAGTCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGCGCTCATCCTTAcct < 1:178136/81‑1 (MQ=255) gTCATCAGTCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGAGCTCATCCTTAcct > 2:178159/1‑81 (MQ=255) tcatcaGTCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGAGCTCATCCTTAcctc > 1:178203/1‑81 (MQ=255) atcaGTCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGCGCTCATCCTTAcctcct > 2:178158/1‑81 (MQ=255) gTCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGAGCTCATCCTTACCTCCTGTgg > 2:178176/1‑81 (MQ=255) gTCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGAGTTTTGGAGCTCATCCTTACCTCCTGTgg > 2:178221/1‑81 (MQ=255) tCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGt < 1:178137/81‑1 (MQ=255) tACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATc < 1:178138/81‑1 (MQ=255) aCATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATc > 1:178173/1‑80 (MQ=255) gATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGc < 2:178139/80‑1 (MQ=255) attattTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGCGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCt > 1:178165/1‑80 (MQ=255) ttttGCCGTTTCGCTACGCATCATTGGCGTTTTGGCGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGa > 1:178166/1‑81 (MQ=255) tttGCCGTTTCGCTACGCATCATTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGa > 2:178164/1‑80 (MQ=255) tttGCCGTTTCGCTACGCATCATTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGa > 2:178163/1‑80 (MQ=255) cGTTTCGCTACGCATCATTGGCGTTTTGGCGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTg > 2:178179/1‑81 (MQ=255) gTTTCGCTACGCATCATTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTgg > 1:178190/1‑81 (MQ=255) tttCGCTACGCATCATTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTgg > 1:178182/1‑80 (MQ=255) ttCGCTACGCATCATTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTGGTa > 2:178168/1‑81 (MQ=255) gCTACGCATCATTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTGGTAttt > 2:178174/1‑81 (MQ=255) tACGCATCATTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTGGTATTTTg > 1:178183/1‑81 (MQ=255) atcatTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTGGTATTTTGCCTcc > 1:178191/1‑81 (MQ=255) atcatTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTGGTATTTTGCCTcc > 1:178180/1‑81 (MQ=255) atTGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTAACGATATCTGGTATTTTGCCTCCGGt < 1:178140/81‑1 (MQ=255) ttGGCGTTTTGGCGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTGGTATTTTGCCTCCGGTg < 1:178142/81‑1 (MQ=255) ttGGCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTGGTATTTTGCCTCCGGTg < 2:178141/81‑1 (MQ=255) ggCGTTTTGGCGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGc > 1:178177/1‑81 (MQ=255) gCGTTTTGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCg < 2:178143/81‑1 (MQ=255) tGGAGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCAtt < 1:178144/80‑1 (MQ=255) | GTTTGCCTGGCGTCATCAGTCCAGCGACGAATACATTGATTATTTTGCCGTTTCGCTACGCATCATTGGCGTTTTGGCGCTCATCCTTACCTCCTGTGGTCTGGCGGCAATCAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATT > NC_000913/933508‑933661 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |