Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,296,329 | 0 | A | C | 28.3% | 123.0 / 13.0 | 53 | intergenic (‑149/‑66) | inaA/yfaH | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (34/4); new base C (0/15); total (34/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.20e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTT > REL606/2296256‑2296403 | cGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacc < 1:1412924/76‑1 (MQ=255) gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca < 2:1636250/76‑1 (MQ=255) aTTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGc < 2:1674772/76‑1 (MQ=255) ttAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcc > 2:885733/1‑76 (MQ=255) tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTGCCCTTCATTTGCAGCGCATAGATTCCACCAGcca < 2:790498/76‑1 (MQ=255) aaTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCCCCAGccac < 2:448554/76‑1 (MQ=255) aaTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccac > 2:561910/1‑76 (MQ=255) gAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactac > 2:981089/1‑76 (MQ=255) cTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTaa > 1:287892/1‑76 (MQ=255) ttttCCCACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAAt < 2:1645568/76‑1 (MQ=255) tttGCCACTACTAGCGTGACCCCGCTTTTACCCTTCATTTCCAGCGAATAGATTCCCCCAGCCACTACCACTAATa < 2:1509107/76‑1 (MQ=255) tttCCCACTACTAGCGTCACACCGCTTTTACCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATa < 2:1509103/76‑1 (MQ=255) gCCGCTCCTAGCGTGACCCCGCTTTTTCCCTTCATTTGCACCGCAGAGATTCCCCCAGCCACTACCACTAATAGTg < 2:408350/76‑1 (MQ=255) ccACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgc > 2:1377470/1‑76 (MQ=255) actactAGCGTGACACCTCTTTTACCCTTCATTTGCAGCGCATAGTTTCCACCAGCCACTACCACTAATAGTgcgc < 1:811655/76‑1 (MQ=255) actactAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAgcgcgc > 1:1008700/1‑76 (MQ=255) ctactaGCGTGACCCCTCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGcc < 1:1374175/76‑1 (MQ=255) ctactaGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGcc < 1:291254/76‑1 (MQ=255) ctactaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGcc > 2:124148/1‑76 (MQ=255) tactgGCGTGCCCCCGCTTTTCCCCTCCATTGGCCGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGccc < 1:539360/76‑1 (MQ=255) actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa > 1:372587/1‑76 (MQ=255) ctaGCGTGACCCCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAg < 1:1303127/76‑1 (MQ=255) ctaGCGGGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAg < 2:159596/76‑1 (MQ=255) taGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGt > 1:841520/1‑76 (MQ=255) taGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGt > 1:1769306/1‑76 (MQ=255) aGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTa > 2:1149169/1‑76 (MQ=255) gCGTGACCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg < 1:1469867/76‑1 (MQ=255) gCGGGACACCCGTTTTCCCCTTGATTTGCAGCGCATAGTTTCCTCCAGCCACTACCACTAATAGTGCGCCCAGTAg < 2:76973/76‑1 (MQ=255) cggTGACCCCCCTTTTCCCCTTCATTGGCAGCGCATGTATCCCACCAGCCACTACCACTAATAGTGCGCCCAGTAg < 1:679774/74‑1 (MQ=255) tGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATa > 1:145863/1‑76 (MQ=255) gACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaa > 1:160196/1‑76 (MQ=255) acacCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaa > 1:1065241/1‑76 (MQ=255) cacCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaaa < 1:738114/76‑1 (MQ=255) cGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtgt > 2:1768610/1‑76 (MQ=255) ttttACCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGaa < 1:199090/76‑1 (MQ=255) tttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt > 1:421679/1‑76 (MQ=255) tttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt > 1:1166765/1‑76 (MQ=255) tttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt > 1:164358/1‑76 (MQ=255) tccccTTCATTGGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt < 2:162793/74‑1 (MQ=255) tACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt > 1:437416/1‑76 (MQ=255) aCCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTa > 2:1131086/1‑76 (MQ=255) cccTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTaa < 2:206550/76‑1 (MQ=255) ccTTCTTTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTCGATAAAATGTGAATTTAAt < 1:20019/76‑1 (MQ=255) ccTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt > 2:1207757/1‑76 (MQ=255) ccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt < 2:1079696/76‑1 (MQ=255) cTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGCGCGCCCAGCAGATAAAATGTGAATTTAATc > 1:83684/1‑76 (MQ=255) tCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg > 2:1376200/1‑76 (MQ=255) atgtGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTc < 2:1500987/73‑1 (MQ=255) aTTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTc > 1:1101731/1‑76 (MQ=255) aTTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTc > 1:1595730/1‑76 (MQ=255) ttGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCgg < 2:539274/76‑1 (MQ=255) tGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGa > 1:626128/1‑76 (MQ=255) cAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATa > 2:1660743/1‑76 (MQ=255) aGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATAc > 2:708347/1‑76 (MQ=255) cATGGACTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCg < 2:390151/76‑1 (MQ=255) aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa > 2:1320252/1‑76 (MQ=255) aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa > 1:398768/1‑76 (MQ=255) aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa > 1:482151/1‑76 (MQ=255) tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt > 1:190977/1‑76 (MQ=255) aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg > 1:1125215/1‑76 (MQ=255) aTTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGaa > 1:856032/1‑76 (MQ=255) caccaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGttt > 1:137619/1‑76 (MQ=255) | CGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTT > REL606/2296256‑2296403 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |