| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 500,213 | 0 | A | C | 31.6% | 62.3 / 10.5 | 38 | H33P (CAC→CCC) | ylbH | conserved protein, rhs‑like |
| Reads supporting (aligned to +/- strand): ref base A (15/11); new base C (0/12); total (15/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.92e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTAC > REL606/500140‑500283 | ggAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCAcc > 1:1087921/1‑76 (MQ=255)ggAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCAcc < 1:348768/76‑1 (MQ=255) aaaCATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACTAGCTGAATGAAGAGAACCCGCATCACCTGCACCAg < 2:417682/76‑1 (MQ=255) aCATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGcc < 1:70587/76‑1 (MQ=255) cATTGCGTGGAGCGTGTAGTATGATGAGTGGGGCAACCAGCTGAATTAAGAGAACCCGCATCACCTGCACCAGCCg < 2:196844/76‑1 (MQ=17) cATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCg > 2:1153409/1‑76 (MQ=255) tGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTAc > 1:1161207/1‑76 (MQ=255) tGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTAc < 2:917589/76‑1 (MQ=255) gTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCgg > 2:1047282/1‑76 (MQ=255) ggAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCCCCTGCACCAGCCGTACCGGCt < 2:1368853/76‑1 (MQ=255) gagaATGAGGGGGGGGGCACCCGGCTGAATGAAGAGACCCCGCATCCCCGGCCCCAGCCGTACCGGCTGCCGGGgc < 2:786343/76‑1 (MQ=16) gAGTATGATGGGTGGGGCAACCAGCTGAATGAAGAGACCCCGCATCCCCTGCCCCAGCCGTACCGGCTGCCGGGgc < 2:1688683/76‑1 (MQ=17) gAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGACCCCGCATCCCCTGCCCCAGCCGTACCGGCTGCCGGGgc < 1:1411921/76‑1 (MQ=17) gAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGACCCCGCATCACCTGCCCCAGCCGTACCGGCTGCCGGGgc < 2:1710418/76‑1 (MQ=255) gAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGACCCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgc < 2:1698805/76‑1 (MQ=255) gAGTATGATGAGTGGGCCACCCAGCTAAATGAAGAGACCCCGCATCCCCTGCCCCAGCCGTACCGGCTGCCGGGgc < 2:1352257/76‑1 (MQ=21) gAGTAGGATGAGGGGGGCAACCAGCTAAATGAAGAGACCCCGCATCACCTGCCCCAGCCGTACCGGCTGCCGGGgc < 2:759817/76‑1 (MQ=17) aGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGACCCCGCATCACCTGCCCCAGCCGTACCGGCTGCCGGGgca < 2:408357/76‑1 (MQ=255) aGTATGATGAGGGGGGCAACCAGCTAAATGAAGAGACCCCGCATCACCTGCCCCAGCCGTACCGGCTGCCGGGgca < 1:1081543/76‑1 (MQ=17) tATGAGGAGTGGGGCAACCAGCTGAATGAAGAGACCCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgcagc < 2:1572204/76‑1 (MQ=255) atgatgAGGGGGGCACCCAGCTGAATGAAGAGACCCCGCATCCCCTGCCCCAGCCGTACCGGCTGCCGGGgcagca < 1:1479304/76‑1 (MQ=17) atgatgAGGGGGGCAACCAGCTAAAGGAAGAGACCCCGCATCCCCTGCCCCAGCCGTACCGGCTGCCGGGgcagca < 1:1633355/76‑1 (MQ=17) gaGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACTAGCCGTACCGGCGGCCGGGGCAGGAGTTTg > 1:502510/1‑76 (MQ=17) gTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGAACCAGCCGTACCGGCTGCCGGGGCAGCAGgatgat > 2:1092026/1‑76 (MQ=255) tGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATa > 2:1087291/1‑76 (MQ=255) ggCAACCAGCTGAATGAAGAGACCCCGCATCCCCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAgg < 1:333038/76‑1 (MQ=255) ggCAACCAGCTGAAGGAAGAAACCCCGCATCACCTCCCCCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAgg < 1:1265685/76‑1 (MQ=21) ggCAACCAGCTAAATGAAGAGACCCCGCATCCCCTGCCCCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAgg < 1:1606651/76‑1 (MQ=21) cAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAg > 2:1505099/1‑76 (MQ=255) cAACCAGCAGAATAAAGAGAACCCGCATCCCCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAg < 2:126734/76‑1 (MQ=18) aCCAGCTGAATGAAGAGACCCCGCATCCCCTGCCCCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTc < 1:415758/76‑1 (MQ=18) aCCAGCTGAATGAAGAGACCCCGCATCACCTGCAACAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTc < 1:1430551/76‑1 (MQ=255) cAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGGATGATAAGGAGTCgg > 2:1329971/1‑76 (MQ=255) cAGCGGAATGAAGAGACCCCGCATCCCCTGCCCCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCgg < 2:573152/76‑1 (MQ=21) aGCTGAATGAAGAGACCCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCggg < 2:1384160/76‑1 (MQ=255) cTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGGCGGGGc > 1:1134709/1‑76 (MQ=255) aaatgaaGAAACCCCGCACCCCCTCCCCCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTg < 2:519290/75‑1 (MQ=17) gaaGAGACCCCGCATCACCTGCACCAGCCATACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGtact < 2:286624/76‑1 (MQ=255) aaGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTCCCGGGGCAGCACTATGNTAAGGAGTCGGGGCTGtacta > 2:278537/1‑76 (MQ=255) agagACCCCGCATCCCCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGtactac < 1:1695867/76‑1 (MQ=255) agaACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGGTGgactacca > 2:475203/1‑74 (MQ=17) gaccccGCATCCCCTGCCCCAGCCGTACCTGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACaac < 2:298842/73‑1 (MQ=17) gaccccGCATCACCTGCCCCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACaac < 1:706928/73‑1 (MQ=255) gaccccGCACCCCCTGCCCCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACaac < 1:1159758/73‑1 (MQ=17) ccccGCACCCCCTGCCCCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACaaccg < 1:1715069/75‑1 (MQ=17) aCCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGGATGATAAGGAGGCGGGGCTGGACTACaaacg > 1:1315138/1‑76 (MQ=17) gCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACaaccggaac > 1:1179167/1‑76 (MQ=255) aTCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGGATGATAAGGAGGCGGGGGTGGAATACcaccgggaccg > 2:480289/1‑76 (MQ=17) tCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACaaccggaaccgg > 2:94803/1‑76 (MQ=255) ccTGCCCCAGCCGTGCCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGtac < 1:511687/76‑1 (MQ=255) | GGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTAC > REL606/500140‑500283 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |