Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,747,903 | 0 | G | T | 20.4% | 129.8 / 22.1 | 44 | A8S (GCA→TCA) | ydiL | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (28/7); new base T (9/0); total (37/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.14e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGC > REL606/1747829‑1747977 | agGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAgg > 1:293791/1‑75 (MQ=255) agGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAgg < 1:1291735/75‑1 (MQ=255) gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa > 1:352321/1‑75 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc > 1:224585/1‑74 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca > 1:793274/1‑74 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcct > 1:523228/1‑73 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtg > 1:1580225/1‑74 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACta > 1:79765/1‑74 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACta > 1:51115/1‑74 (MQ=255) tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtagc > 1:1326920/1‑72 (MQ=255) tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtagc > 1:179872/1‑72 (MQ=255) aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc < 1:533954/75‑1 (MQ=255) aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc > 1:1265847/1‑75 (MQ=255) gTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat > 1:1241893/1‑75 (MQ=255) gTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat < 1:494671/75‑1 (MQ=255) tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat > 2:812533/1‑75 (MQ=255) atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt > 2:229730/1‑75 (MQ=255) ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt > 2:1411210/1‑75 (MQ=255) ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt > 1:617868/1‑75 (MQ=255) ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt > 2:422925/1‑75 (MQ=255) ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt > 2:707854/1‑75 (MQ=255) tataTTTTATTGGTTATTTCTTGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt > 2:799473/1‑75 (MQ=255) tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc > 2:81869/1‑75 (MQ=255) atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc > 1:949315/1‑75 (MQ=255) tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt > 1:635500/1‑75 (MQ=255) tttATTGTTTATTTCATCCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt < 1:876272/75‑1 (MQ=255) ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg > 2:55111/1‑75 (MQ=255) aTTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc > 2:958907/1‑75 (MQ=255) ttGTTTATTTCATGCAGCCATTACACGATGAACGCTTATCAACTCCAGGCACTCCGCCATATTTTTGCCATGACt > 2:238604/1‑75 (MQ=255) tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCCCTCCGCCATATTTTTGCCATGACTa < 1:748812/75‑1 (MQ=255) ttttttTCATGCAGCCATTACACGATGAACGCTTATGACCTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt < 1:1068338/75‑1 (MQ=255) ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg > 2:882851/1‑75 (MQ=255) tATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga < 2:616654/75‑1 (MQ=255) aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat < 1:100629/75‑1 (MQ=255) aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat > 1:481394/1‑75 (MQ=255) tttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg > 1:1549269/1‑75 (MQ=255) tCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa > 2:502175/1‑75 (MQ=255) cAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcgc > 1:192139/1‑75 (MQ=255) aGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg > 1:1476008/1‑75 (MQ=255) aGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg > 2:583417/1‑75 (MQ=255) ccATTACACGATGAACGCTTATGACCTCCAGCAACGCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAc < 1:619961/75‑1 (MQ=255) acacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGa > 2:1487093/1‑75 (MQ=255) cacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAt > 2:720902/1‑75 (MQ=255) aTGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGcc < 1:607670/75‑1 (MQ=255) gCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAcc > 2:544069/1‑75 (MQ=255) cTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCg > 2:812049/1‑75 (MQ=255) ttATGAACTCCAGGCCCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCgg < 1:155901/75‑1 (MQ=255) ttATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCgg < 1:211921/75‑1 (MQ=255) aTGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTg > 2:707025/1‑75 (MQ=255) aTGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAACGGTg > 2:138344/1‑75 (MQ=255) accTCCAGGCCCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACa < 1:1172241/73‑1 (MQ=255) ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa > 1:962183/1‑75 (MQ=255) gCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGc < 1:1370150/75‑1 (MQ=255) | AGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGC > REL606/1747829‑1747977 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |