Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,383,157 | 0 | T | G | 53.8% | 4.7 / 32.7 | 26 | G224G (GGT→GGG) | ycjX | conserved protein with nucleoside triphosphate hydrolase domain |
Reads supporting (aligned to +/- strand): ref base T (3/9); new base G (14/0); total (17/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.04e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCACCACTGTAAACAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGCCGATAA > REL606/1383085‑1383231 | ccaccaCTGTAAACAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGgcg > 1:301532/1‑76 (MQ=255) ccaccaCTGTAAACAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGgcg > 1:1011117/1‑76 (MQ=255) cacTGTAAACAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCcg < 1:1192038/76‑1 (MQ=255) acTGTAAACAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCcgc < 2:1134330/76‑1 (MQ=255) tGTAAACAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCgcgc < 2:1064085/76‑1 (MQ=255) tAAACAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGGgcgcgcgccg > 2:295193/1‑74 (MQ=255) tAAACAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGGCCccgccg > 2:1052976/1‑74 (MQ=255) tAAACAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCccgcgg > 2:62592/1‑74 (MQ=255) aaaCAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCGCGCTGc > 1:1753146/1‑76 (MQ=255) gcagGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAAttc < 1:1841367/76‑1 (MQ=255) gcagGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAAttc < 1:340315/76‑1 (MQ=255) cagGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAAttct < 2:683498/76‑1 (MQ=255) gCACTTTATTCAGCCTGGGCGCTTTGTTTTGCCGGGAGATATGGCAGGTGCGCCCGCGCGGGAATTCTTCCCGTgg > 2:2190/1‑76 (MQ=255) gCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCCCCCTGCCATTCTTTCCCTgg > 2:1348851/1‑76 (MQ=255) ttATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGGCCCCGCGGCCAATCTTTCCccggccggc > 1:804458/1‑75 (MQ=255) aTTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCGCGCCGCCATTCTTCCCCTGGCCGGGTg > 2:1647863/1‑76 (MQ=255) ttCAGCCGGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGt < 2:808541/76‑1 (MQ=255) cAGCATGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCCCGCCGCAATTTTTCCCCTGGCCGGGTGTCg > 1:334762/1‑76 (MQ=255) aGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCCCGCTGCAATTCTTCCCGTGGCCGGGTGTCgg > 1:447578/1‑75 (MQ=255) gCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCCCCCTGCAATTCTTCCCCTGGCCGGCTGGCggt > 1:586737/1‑74 (MQ=255) cTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCGCGCCGCAATTCTTCCCGTGGCCGGATGTCGATAc > 2:1481307/1‑76 (MQ=255) tGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGGTGTCGATAcc > 1:1847071/1‑76 (MQ=255) tGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCGCGCTGCAATTCTTCCCCTGGCCGGATGTCGATAcc > 1:1554550/1‑76 (MQ=255) ttGTCTTGCCGGGAGATATGGCAGGGGCGCCCGCGCTGCAATTTTTTCCCTGGGCGGGTGTCGGTACCTggggggg > 1:416634/1‑73 (MQ=255) gTCTTGCCGGGAGATATGGCAGGGGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGGTGGCGATACCTGGGGCgggg > 1:52417/1‑73 (MQ=255) cgtGCCGGGAGATGTGGGAGGTGCGCCCGCGCTGCAATTCTTCCCTTGGCCGGATGTCGATACCTGGGGCGAGTcc < 1:742513/74‑1 (MQ=255) cGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACt < 2:918689/76‑1 (MQ=255) gaTATGGCAGGGGCGCCCCCGCCGCAATTCTTCCCGTGGCCGGGGGGCGATAACTGGGGGGGGGCCAAACAGgcgc > 2:1538204/1‑76 (MQ=255) atGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAgg > 1:1685573/1‑76 (MQ=255) tGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGc > 1:842395/1‑76 (MQ=255) cAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGCCGa > 2:997681/1‑76 (MQ=255) gTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGCCGATaa < 1:1496754/76‑1 (MQ=255) | CCACCACTGTAAACAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGCCGATAA > REL606/1383085‑1383231 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |