Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,568,967 | 0 | A | C | 22.0% | 135.5 / 10.9 | 60 | G1380G (GGT→GGG) | yfhM | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (34/12); new base C (0/13); total (34/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.28e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGC > REL606/2568892‑2569041 | cTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCa < 1:385809/76‑1 (MQ=255) cATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCAcc > 1:81864/1‑76 (MQ=255) tCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCa < 1:640005/76‑1 (MQ=255) cGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAg > 1:1062440/1‑76 (MQ=255) gttAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGcc > 1:494906/3‑76 (MQ=255) gTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCa < 1:66436/76‑1 (MQ=255) tAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCat > 2:1510489/1‑76 (MQ=255) tAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCat > 1:45484/1‑76 (MQ=255) aGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCata < 1:1485748/76‑1 (MQ=255) aGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCata > 2:415797/1‑76 (MQ=255) gCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCatat < 1:456509/76‑1 (MQ=255) cAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATc < 1:984118/76‑1 (MQ=255) cAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATc < 2:601730/76‑1 (MQ=255) cAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATc < 1:1680657/76‑1 (MQ=255) tttATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGc > 1:1559709/1‑76 (MQ=255) aTTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCtca < 1:447431/76‑1 (MQ=255) ttttCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCtcat < 2:1487049/76‑1 (MQ=255) tttCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCtcatc > 2:181252/1‑76 (MQ=255) tttCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCtcatc > 2:936204/1‑76 (MQ=255) ttcttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCtcatca > 2:976265/1‑76 (MQ=255) cttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACt < 1:223501/76‑1 (MQ=255) cttcCAGCAAGGAGAGCATAAACGCGTTGCCGCCCAGGAAACTACCGTAATCACCCAGCCATATCCGCTCATCACt < 2:21827/76‑1 (MQ=255) ctcccAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACt < 2:96906/73‑1 (MQ=255) ttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTa < 1:792748/76‑1 (MQ=255) ccAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCTCCCAGCCATATCCGCTCATCACTAtt < 1:1313320/76‑1 (MQ=255) ccAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTAtt < 2:593914/76‑1 (MQ=255) ccAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTAtt < 1:62563/76‑1 (MQ=255) gCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTAcg < 1:715822/76‑1 (MQ=255) cAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACACAGCCATATCCGCTCATCACTATTAcgc > 1:1616957/1‑76 (MQ=255) aaGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTAcgcg > 1:1034220/1‑76 (MQ=255) aaGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTAcgcg > 2:1691282/1‑76 (MQ=255) aaGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCAAGCCATATCCGCTCATCACTATTAcgcg > 1:1502773/1‑76 (MQ=255) ggAGAGCATTACCGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACgcggc < 1:1216180/76‑1 (MQ=255) ggAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACgcggc > 2:1452976/1‑76 (MQ=255) agCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGttt > 1:1001770/1‑76 (MQ=255) cATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTc > 2:543111/1‑76 (MQ=255) aTTACCGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTca < 2:760452/76‑1 (MQ=255) ttAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcag > 2:408551/1‑76 (MQ=255) tAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagc > 2:135348/1‑76 (MQ=255) aaCGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcg > 2:1401851/1‑76 (MQ=255) aCGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgc > 2:1033517/1‑76 (MQ=255) aCGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgc > 1:814706/1‑76 (MQ=255) cgcgTTGTCGCGCAGTGAACTACCGTAATCACCCACCCATATCCCCCCACAACAATTACCCGCCGTTTTcaccgca > 1:1192180/1‑75 (MQ=255) gcgTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgcca < 2:350313/76‑1 (MQ=255) gTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccagc < 2:690741/76‑1 (MQ=255) tGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccagcgc > 1:864309/1‑76 (MQ=255) gTCGCGCAGGAAACTCCCGTAACCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCa < 2:1377764/76‑1 (MQ=255) gcgcAGTGAACTCCCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATc < 2:415587/76‑1 (MQ=255) cgcAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCg < 2:1374072/76‑1 (MQ=255) cAGTGAACTCCCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGct < 1:23110/76‑1 (MQ=255) cAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGct > 2:464482/1‑76 (MQ=255) gAACTCCCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGcttctt < 2:1682805/76‑1 (MQ=255) aaCTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCttcttc > 2:8793/1‑76 (MQ=255) aaCTACCGTAACCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCttcttc < 1:1377426/76‑1 (MQ=255) ccGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTcaccac < 1:577736/76‑1 (MQ=255) ccGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTcaccac < 1:836792/76‑1 (MQ=255) ccaaaTCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCAcg < 2:524653/73‑1 (MQ=255) gTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCAcgc > 2:1508390/1‑76 (MQ=255) gTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCAcgc > 2:1106878/1‑76 (MQ=255) gTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGACATCGCTTCTTCACCAcgc > 2:105884/1‑76 (MQ=255) tAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCAcgcg > 2:1525460/1‑76 (MQ=255) tAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCAcgcg > 1:1085111/1‑76 (MQ=255) aaTCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGt > 1:498295/1‑76 (MQ=255) aaTCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGt > 2:575292/1‑76 (MQ=255) cACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGATTCTTCACCACGCGTCGc > 2:126304/1‑76 (MQ=255) | CTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGC > REL606/2568892‑2569041 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |