Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,296,329 | 0 | A | C | 26.2% | 141.0 / 15.5 | 61 | intergenic (‑149/‑66) | inaA/yfaH | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (44/1); new base C (0/16); total (44/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.38e-14 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTA > REL606/2296254‑2296404 | tttgTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca > 2:373934/4‑76 (MQ=255) gtCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca < 1:195396/76‑1 (MQ=255) gtCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca > 1:58936/1‑76 (MQ=255) tCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccac < 2:1675296/76‑1 (MQ=255) cGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacc > 2:1244341/1‑76 (MQ=255) cGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacc > 1:1810953/1‑76 (MQ=255) gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca < 1:757569/76‑1 (MQ=255) ttAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcc > 1:1114294/1‑76 (MQ=255) aaTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccac > 1:202990/1‑76 (MQ=255) aaTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccac < 1:1244341/76‑1 (MQ=255) aTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccact > 1:1664659/1‑76 (MQ=255) gAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactac > 1:1582731/1‑76 (MQ=255) aattttCTTTTGCCACTACTAGCTTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGccactac < 2:776068/73‑1 (MQ=255) tttCTTTTGCCCCTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATCCCCCCAGccactaccac < 1:2548/76‑1 (MQ=255) ttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccact > 1:1674585/1‑76 (MQ=255) ttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccact > 1:1121386/1‑76 (MQ=255) ntCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccact > 1:1344126/2‑76 (MQ=255) ttttttCCCACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGccactaccacta < 2:1154088/74‑1 (MQ=255) ttttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTGGCAGCGCATAGATTCCCCCAGCCACTACCACTAAt < 1:1253041/76‑1 (MQ=255) ttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAg > 2:900803/1‑76 (MQ=255) cACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgcg > 2:195993/1‑76 (MQ=255) actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa > 2:1801304/1‑76 (MQ=255) ctaGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAg < 1:517233/76‑1 (MQ=255) ctaGCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAg < 1:272207/76‑1 (MQ=255) taGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGt > 2:546069/1‑76 (MQ=255) taGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGt > 2:1720374/1‑76 (MQ=255) aGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTa < 1:825912/76‑1 (MQ=255) aGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTa < 2:1428099/76‑1 (MQ=255) aGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTa > 1:55036/1‑76 (MQ=255) cGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGa > 1:23764/1‑76 (MQ=255) tGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATa > 2:546838/1‑76 (MQ=255) gACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaa > 2:726694/1‑76 (MQ=255) gACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaa > 1:829240/1‑76 (MQ=255) ccccGTTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaaa < 2:1357256/74‑1 (MQ=255) cccGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAt < 2:1683962/75‑1 (MQ=255) cGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtgt > 2:245724/1‑76 (MQ=255) cGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtgt > 2:723047/1‑76 (MQ=255) gtttttCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtgtg < 1:1554841/74‑1 (MQ=255) ttttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGaa > 1:1255975/1‑76 (MQ=255) tttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt > 2:1361053/1‑76 (MQ=255) ttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAtt > 2:375016/1‑76 (MQ=255) tACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt > 2:1407995/1‑76 (MQ=255) tACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt > 1:1851034/1‑76 (MQ=255) aCCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTa > 2:187530/1‑76 (MQ=255) cccTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTaa > 2:836804/1‑76 (MQ=255) cccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTaa < 2:1674585/76‑1 (MQ=255) ccTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt > 2:378647/1‑76 (MQ=255) ccTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt > 2:378657/1‑76 (MQ=255) ccTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt > 2:378674/1‑76 (MQ=255) ttCATTTGCAGCGCAAAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCa < 1:732694/76‑1 (MQ=255) aTTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTc > 1:1071881/1‑76 (MQ=255) ttGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCgg > 2:1497264/1‑76 (MQ=255) gCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGAt > 2:1069921/1‑76 (MQ=255) cAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATa > 2:303775/1‑76 (MQ=255) cAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGGGAATTTAATCAGTCGGATa > 2:1036261/1‑76 (MQ=255) aGCGCATAAATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATAc < 1:919576/76‑1 (MQ=255) cgcATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAg < 1:1037740/76‑1 (MQ=255) cgcATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAg < 2:826708/76‑1 (MQ=255) gcATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGc > 1:425928/1‑76 (MQ=255) gcATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGc < 1:721936/76‑1 (MQ=255) cATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCg > 2:1527207/1‑76 (MQ=255) aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa > 2:668795/1‑76 (MQ=255) aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg > 2:1693314/1‑76 (MQ=255) aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg > 1:403414/1‑76 (MQ=255) aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg > 1:1116794/1‑76 (MQ=255) aGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg < 2:1121386/76‑1 (MQ=255) gATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa > 1:447513/1‑76 (MQ=255) gATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa > 2:1166071/1‑76 (MQ=255) ttCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAg > 1:1610523/1‑76 (MQ=255) tCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGt > 2:1490409/1‑76 (MQ=255) accaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTa > 1:385348/1‑76 (MQ=255) accaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTa > 1:1241531/1‑76 (MQ=255) | GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTA > REL606/2296254‑2296404 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |