Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,783,565 | 0 | A | C | 34.6% | 178.1 / 45.0 | 101 | T744P (ACA→CCA) | ECB_03521 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (43/23); new base C (0/35); total (43/58) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.24e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 229 total aligned reads displayed.
CTGTTTGAAGATGATGGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTGATAA > REL606/3783490‑3783640 | cTGTTTGAAGATGATGGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGGGGGGGGCTGca > 1:3346844/1‑76 (MQ=255) gTTTGAAGATGATGGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGcaca < 2:810850/76‑1 (MQ=255) gTTTGAAGATGATGGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGcaca < 2:4958345/76‑1 (MQ=255) tGAAGATGATGGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGcc < 2:3437956/76‑1 (MQ=255) tGAAGATGATGGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGcc < 1:3541284/76‑1 (MQ=255) gAAGATGATGGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCa < 1:6696183/76‑1 (MQ=255) aaGATGATGGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGGGGTGGGGTGCACAGcccc > 1:738514/1‑74 (MQ=255) aaGATGATGGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTGTATGGGAGGGGGAATGCccacccac > 2:3397121/1‑76 (MQ=255) aGATGATGGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGGGGTGGACTGCGCAGcccca > 2:2769111/1‑73 (MQ=255) aGATGATGGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTTAATGGGGGGTGGGCTGCACCGCCAcc > 1:6520213/1‑75 (MQ=255) aGATGATGGC‑AAACTTGGGGGATATACTGAGGGTAATGCTCTGTGTCTTGAAGGGGAGCGGGCCTCCCCAGCCACt < 2:6660018/76‑1 (MQ=255) atgatgGCAAAATTG‑GGGTATATACTGAGGGTAATTCTCTGTGGCTTGAAGGGGAGCGGGCCCCCCCAGCCACTAc < 2:4405190/76‑1 (MQ=255) tgatgGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGGGGTGGGCTGCACAGGCACcaca > 1:5293633/1‑76 (MQ=255) tgatgGCAAAATTT‑GGGGATATACTGAAGGTAATCCTCTGTGGCTTGAATGGGGGCGGGCCGCCCCAGCCACTACa < 2:5528833/76‑1 (MQ=255) tgatgGCAAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTGGAAGGGGGGCTGGACTCCCCAGCCACTACa < 2:3647108/76‑1 (MQ=255) atgGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAAt < 1:1599049/76‑1 (MQ=255) tgGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGGGGTGGGCTGCACAGCCCCTACaaat > 1:4944019/1‑74 (MQ=255) tgGCGAAACTG‑GGGGATATACTGAAGGTATTGCTCTGTGTCTTGAAGGGGAGCTGGACTCCCCAGCCACTACAAtt < 2:5321372/76‑1 (MQ=255) ggCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGGGGTGGGCTGCACCGCCCCCACAATTg > 1:6315019/1‑76 (MQ=255) gCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGa < 1:3774178/76‑1 (MQ=255) aaaaaCTT‑GGGAATATACTGAAGGTAGTTCTCTGTGGCTTGAAGGGGAGCGGGACGCCCCAGCCACTACAATTGAg < 2:638097/74‑1 (MQ=255) aaaCTT‑GGGTATATACTGAAGGTAATGCTCTGTGGCTGGAAGGGGAGCTGGCCTGCCCAGCCACTACAATTGAGtt < 2:5052005/76‑1 (MQ=255) aaCTG‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTGAAGGGGAGCTGGACTGCACAGCCACTACAATTGAGTTa < 1:488161/76‑1 (MQ=255) aCTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGGGGCGGGCTGGACCGCCCCCACAAATTAGttta > 2:168775/1‑74 (MQ=255) aCTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTTAATGGGGGCTGGGCTGCACAGCCACTACAATTGAGttta > 1:1378362/1‑74 (MQ=255) cTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAg < 1:2555385/76‑1 (MQ=255) tt‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTTAATGGGGGCTGGACTGCACAGCCACTACAATTTAGTTAAGa > 2:5161384/1‑76 (MQ=255) t‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTTAATTGGAGCTGGAGTGCGCCACCACTACAAATTAGGTAAGaa > 1:4260373/1‑76 (MQ=255) tgaggATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCGAAACTCCCCAGCCACTACAATTGAGTTAAGaa < 2:5864890/73‑1 (MQ=255) ggggATATACTGAAGGTAATGCTCTGTGGCTTTAATGGGAGCTGGACTGCACAGCCACTACAAATGAGGTAAGaaa > 1:6407387/1‑75 (MQ=255) ggggATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAAt < 1:257932/76‑1 (MQ=255) gggATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATa < 1:1336431/76‑1 (MQ=255) gggATATACTGAAGGTAATGCTCTGTGGCTTGAAGGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATa < 1:6084265/76‑1 (MQ=255) ggATATACTGAAGGTAATGCTCTGTGGCTTTAATGGGGGCTGGGCTGCACAGCCACCACAAATGAGGTAAGAaaaa > 2:112062/1‑76 (MQ=255) gATATAGGGAAGGTAATGCTGTGTGGTTGGAAGGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATaaa < 2:3028400/76‑1 (MQ=255) gATATACTGAAGGTAATGCTCTGTGGCTTTAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATaaa > 1:6937907/1‑76 (MQ=255) atatACTGAAGGTAATGCTCTGTGGCTTTAATGGGGGGTGGGCTGGCCAGCCCCCACAATTTAGTTAaaaaaaaaa > 1:4707938/1‑75 (MQ=255) tataCTGAAGGTAATGCTGTGTGGCTTGAAGGGGAGCGGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAaca < 2:2599759/76‑1 (MQ=255) ataCTGAAGGTAATGCTCTGTGGCTTTAATGGGAGCTGGGCTGCACAGCCACTACAAATGAGTTAaaaaaaaaaaa > 2:6263334/1‑73 (MQ=255) ataCTGAAGGTAATGCTCTGTGGCTTTAATGGGAGCTGGACTGCACAGACACTACAATTGAGTTAAGAATAAaaac > 2:14764/1‑76 (MQ=255) ataCTGAAGGTAATGCTCTGGGGCTTGAAGGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAacac < 1:1272685/76‑1 (MQ=255) taCGGAAGGTAATGCTCTGTGGCTTGAATGGGAGCCGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACt < 1:1394500/76‑1 (MQ=255) cTGAAGGTATTGCTCTGTGGCTGGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCa < 2:4587922/76‑1 (MQ=255) cTGAAGGTAATGCTCGGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCa < 1:6478724/76‑1 (MQ=255) aaGGTAATGCTCTGTGGCTTTAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATgg > 1:1010681/1‑76 (MQ=255) tAATGCTCTGTGGCTTTTATGGGGGCTGGGCTGCACAGCCCCTACAATTGAGTTAAGAATAAACACTCATGgaggt > 1:633443/1‑74 (MQ=255) tAATGCTCTGTGGCTTTAATGGGGGCTGGGCTGCACAGCCACTACAAATGAGTTAAGAATAAACACTCATGGAGAt > 2:4242855/1‑76 (MQ=255) tAATGCTCTGTGGCTTGAAGGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGAt < 1:5975951/76‑1 (MQ=255) aaTGCTCTGTGGCTTTTATGGGGGCTGGGCTGCACCGCCCCCACAATTGAGTTAAGAAAAAAAAATCATGGAGGtt > 2:4079673/1‑76 (MQ=255) aaTGCTCTGTGGCTTTTATGGGGGCTGGGCTGCACAGCCCCTACAAATGAGTTAAGAATAAAAACTCCTTGAGAtt > 2:2545007/1‑76 (MQ=255) aaTGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGAtt < 1:4128396/76‑1 (MQ=255) aTGCTCTGTGGCTTGAATGGGAGCGGGACTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTa < 2:6146010/76‑1 (MQ=255) tGCTCTGTGGTTTGAAGGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTAt < 1:3104389/76‑1 (MQ=255) gCTCTGTGGTTTGAAGGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATc < 2:2493403/76‑1 (MQ=255) gCTCTGTGGCTTTTATGGGGGGTGGGCTGCACAGGCCCTACAAATTAGTTAAAAATAAAAAATCCCGGGGATTAtt > 2:1228265/1‑75 (MQ=255) ctctGTGGCTTTAATGGGGGCTGGGCTGCACAGCCACTACAATTGAGTTTAGAATAAACACTCATGGAGATTATCg > 2:4868293/1‑76 (MQ=255) tctGTGGCTGGAAGGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGc < 1:2886501/76‑1 (MQ=255) gtgGCTTTAATGGGGGCTGGACTGCACAGCCACCACAATTGAGTTAAGAATAAAAACTCATGGAGATTATCGCCCt > 2:1765328/1‑76 (MQ=255) gtgGCTTGAAGGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCt < 1:1121219/76‑1 (MQ=255) tgGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTg > 2:1800531/1‑76 (MQ=255) tgGCTGGAAGGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTg < 1:5515177/76‑1 (MQ=255) ggCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGc > 1:1992887/1‑76 (MQ=255) gCTTGAATGGGAGCTGCACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCa < 2:1758853/76‑1 (MQ=255) gCTTAAAGGGGAGCTGGCCTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCa < 2:3116831/76‑1 (MQ=255) cTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCAt > 2:4531303/1‑76 (MQ=255) ttGAATGGGAGCTGGACTGCACAGCCACTACAACTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATg > 2:4358684/1‑76 (MQ=255) ttGAAGGGGAGCTGGCCTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATg < 2:6506549/76‑1 (MQ=255) ttGAAGGGGAGCTGGACTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATg < 1:6147645/76‑1 (MQ=255) tGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATgg > 2:2923817/1‑76 (MQ=255) tGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATgg > 1:5179547/1‑76 (MQ=255) tGAAGGGGAGCGGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATgg < 2:1873369/76‑1 (MQ=255) gAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATggg > 2:4714035/1‑76 (MQ=255) gAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATggg > 1:50757/1‑76 (MQ=255) gAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATggg > 2:2180699/1‑76 (MQ=255) gAATGGGAGCGGGACTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATggg < 1:5471663/76‑1 (MQ=255) gAATGGGAACTGGCCTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATggg < 1:3263032/76‑1 (MQ=255) aaTGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGa > 1:2283254/1‑76 (MQ=255) aaTGGGAGCGGGCCTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGa < 2:3856771/76‑1 (MQ=255) aaTGGGAGCGGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGa < 1:2654039/76‑1 (MQ=255) gggggAGCGGGCCTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGaa < 2:3989661/74‑1 (MQ=255) aggggAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGaa < 2:5467221/74‑1 (MQ=255) ggCTGGACTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACAtt < 2:3220529/75‑1 (MQ=255) gCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTg > 1:2219878/1‑76 (MQ=255) cTGGCCTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGa < 2:3520044/76‑1 (MQ=255) cTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGa > 1:3652223/1‑76 (MQ=255) cTGAACTACACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGa < 1:1128486/76‑1 (MQ=255) tGGACTCCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGaa < 2:2432442/76‑1 (MQ=255) ggACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGaaa > 2:5603910/1‑76 (MQ=255) gccTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAg < 2:6597485/74‑1 (MQ=255) gACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAg > 1:6536314/1‑76 (MQ=255) aCTCCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGt < 2:5926623/76‑1 (MQ=255) ctccacaGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTg < 2:5185142/73‑1 (MQ=255) tGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTGa > 1:3947832/1‑76 (MQ=255) tGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTGa > 1:6362762/1‑76 (MQ=255) cacaGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTGATa < 2:4883423/76‑1 (MQ=255) cacaGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTGATa > 2:1377412/1‑76 (MQ=255) cacaGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTGATa > 2:2322633/1‑76 (MQ=255) acaGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTGATaa > 2:6132566/1‑76 (MQ=255) acaGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTGATaa > 2:2046074/1‑76 (MQ=255) acaGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTGATaa > 2:2128253/1‑76 (MQ=255) | CTGTTTGAAGATGATGGCGAAACTT‑GGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTGATAA > REL606/3783490‑3783640 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |