Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,568,967 | 0 | A | C | 20.0% | 155.2 / 10.6 | 62 | G1380G (GGT→GGG) | yfhM | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (33/15); new base C (0/12); total (34/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.24e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCA > REL606/2568892‑2569042 | cTCATCAGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCa < 2:741182/76‑1 (MQ=255) cATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTAGCGTAATCAcc > 2:1386655/1‑76 (MQ=255) cATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCAcc > 2:1055500/1‑76 (MQ=255) tCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCa < 2:473965/76‑1 (MQ=255) tCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCa < 2:1677191/76‑1 (MQ=255) ggTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGcc < 1:1514409/76‑1 (MQ=255) ggTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGcc < 1:676161/76‑1 (MQ=255) tAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCat > 2:1864116/1‑76 (MQ=255) cAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATc > 2:1361538/1‑76 (MQ=255) aGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATcc < 2:738977/76‑1 (MQ=255) aGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATcc < 2:307458/76‑1 (MQ=255) gTTTTTTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCg < 2:1888023/76‑1 (MQ=255) tttATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCGGc > 2:1404497/1‑76 (MQ=255) ttATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCt > 2:437533/1‑76 (MQ=255) ttATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCt > 1:630266/1‑76 (MQ=255) ttATTTTCTCCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCt < 1:509843/76‑1 (MQ=255) aTTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCtca > 1:517611/1‑76 (MQ=255) ttcttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCtcatca > 2:1092010/1‑76 (MQ=255) cttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACt > 1:937472/1‑76 (MQ=255) cttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACt > 1:1182210/1‑76 (MQ=255) cttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACt > 1:244755/1‑76 (MQ=255) cttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACt > 1:281577/1‑76 (MQ=255) tcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTAt < 2:42157/76‑1 (MQ=255) tcCAGCAAGGAGAGCATTAACGCGTTGGCGCGCGGGGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTAt < 1:352054/76‑1 (MQ=255) ccAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTAtt > 2:1871482/1‑76 (MQ=255) aGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTAc > 1:697992/1‑76 (MQ=255) gCAAGGAGAGCATTAACGCGTTGTCGCCCAGTAAACTACCGTAATCGCCCAGCCATATCCGCTCATCACTATTAcg < 2:1182828/76‑1 (MQ=255) gCAAGGAGAGCATCAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTAcg < 2:1706741/76‑1 (MQ=255) cAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTAcgc > 1:468027/1‑76 (MQ=255) aaGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTAcgcg < 1:91683/76‑1 (MQ=255) ggAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACgcggc > 1:1329604/1‑76 (MQ=255) agagCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGt > 1:1206322/1‑76 (MQ=255) gagCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGtt < 1:1642441/76‑1 (MQ=255) aagCATTAACGCGTTTTCGCGCAGTGAACTCCCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGtt < 1:985133/75‑1 (MQ=255) gCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGtttt > 1:1521848/1‑76 (MQ=255) cATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTc < 1:1155822/76‑1 (MQ=255) cATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTc < 1:1957845/76‑1 (MQ=255) aTTAACGCGTTGTCGCGCAGTGAACTACCGTAATCGCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTca > 2:1792461/1‑76 (MQ=255) ccgcgTTGTCGCGCAGTGACCTACCGTAATCCCCCAGCCATAGCCGCTCATCACTATTACGCGGCGTTTTcagcgc < 2:636966/75‑1 (MQ=255) cgcgTTGTCGTGCAGCGAACCACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgcc < 2:1386261/76‑1 (MQ=255) gcgTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgcca < 2:1494118/76‑1 (MQ=255) cgTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccag < 1:411691/76‑1 (MQ=255) cgTTGTCGCGCAGTAAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccag < 2:53372/76‑1 (MQ=255) gTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccagc > 1:1672865/1‑76 (MQ=255) gTTGTCGCGCAGTAAACTCCCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccagc < 1:580688/76‑1 (MQ=255) ttGTCGCGCAGTGAACTACCATAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccagcg < 1:898872/76‑1 (MQ=255) tGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccagcgc > 1:969529/1‑76 (MQ=255) cgcgCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAAt > 1:550753/1‑76 (MQ=255) gcgcAGTGACCCCCCCTATTCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATc < 2:19277/76‑1 (MQ=255) gcgcAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGAGGAGATTACAGCGCCAGCGCAATc > 1:928817/1‑76 (MQ=255) gcgcAGTGAACTACCGTAACCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATc < 1:860258/76‑1 (MQ=255) cgcAGCGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGCTTTCAGCGCCAGCGCCATCg > 2:356832/1‑76 (MQ=255) gcAGTAAACTACCGTAACCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGc < 1:1092010/76‑1 (MQ=255) gTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGcttc > 2:1067077/1‑76 (MQ=255) gAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGcttctt > 1:1522845/1‑76 (MQ=255) aCTCCCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCNCAATCGCTTCTTca < 1:1087043/76‑1 (MQ=255) aCTCCCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTca < 1:1598095/76‑1 (MQ=255) aCTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTca < 1:1827262/76‑1 (MQ=255) aCTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTca < 1:1397273/76‑1 (MQ=255) aCTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTca < 1:1987271/76‑1 (MQ=255) tACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTcacc > 2:1119080/1‑76 (MQ=255) tACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTcacc > 1:937746/1‑76 (MQ=255) tACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTCCTTcacc > 2:142452/1‑76 (MQ=255) aCCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTcacca < 1:357586/76‑1 (MQ=255) aCCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTcacca > 1:1496195/1‑76 (MQ=255) tAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCAcgcg < 1:1522144/76‑1 (MQ=255) tAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCAcgcg > 1:1700523/1‑76 (MQ=255) aaTCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGt > 2:459392/1‑76 (MQ=255) aaTCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGt > 1:1643881/1‑76 (MQ=255) tCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCg > 2:37179/1‑76 (MQ=255) aCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCa > 1:1464282/1‑76 (MQ=255) | CTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCA > REL606/2568892‑2569042 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |