Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,783,565 | 0 | A | C | 36.0% | 150.9 / 34.4 | 89 | T744P (ACA→CCA) | ECB_03521 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (36/21); new base C (0/32); total (36/53) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.86e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 180 total aligned reads displayed.
TGTTTGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTGATA > REL606/3783491‑3783639 | tGTTTGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGGGGGGGGCTGcac > 1:3494935/1‑76 (MQ=255) tGTTTGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGcac < 1:4711730/76‑1 (MQ=255) tGTTTGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGcac < 2:5154222/76‑1 (MQ=255) gTTTGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGTCTTGAATGGGAGCTGGACTGcaca < 2:3750381/76‑1 (MQ=255) gTTTGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTTAATGGGAGGTGGGCTGcaca > 2:285702/1‑76 (MQ=255) gTTTGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGcaca < 2:176286/76‑1 (MQ=255) tttGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGGGGTGGACTGcccac > 1:3670399/1‑75 (MQ=255) tttGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAg < 2:4396270/76‑1 (MQ=255) tGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGGGGTGGGCTGCCCAAcc > 2:302477/1‑76 (MQ=255) tGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGGGGTGGGCTGCACAAcc > 1:10486/1‑76 (MQ=255) tGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGcc < 1:1454843/76‑1 (MQ=255) tGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGcc < 1:2508799/76‑1 (MQ=255) gAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCa < 2:3129214/76‑1 (MQ=255) aaGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATTGGGGGGGGGCTGCACAGacac > 1:217876/1‑76 (MQ=255) aaGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAAGGGGAGCTGGACTGCACAGCCAc < 2:3610293/76‑1 (MQ=255) atgatgGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACCGCCCAACcaccac > 2:5004373/1‑76 (MQ=255) gatgGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGAGCTGGGATGCCCAGCCACCACaa > 2:2435172/1‑76 (MQ=255) gatgGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACaa < 1:439323/76‑1 (MQ=255) tgGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGGGCTGGGCTGCACAGCCACTACaaat > 1:2650026/1‑74 (MQ=255) ggCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTg < 1:843624/76‑1 (MQ=255) ggCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTg < 1:157489/76‑1 (MQ=255) ggCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAAGGGGAGCGGGACTCCCCAGCCACTACAATTg < 2:1451565/76‑1 (MQ=255) gCGAAACTTGGGGATATACGGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCCCAGCCACTACAATTGa < 1:4667160/76‑1 (MQ=255) cGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAg < 1:2429073/76‑1 (MQ=255) cGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAg < 1:1885610/76‑1 (MQ=255) cGAAACTGGGGGATATACTGAAGGAAATGCTCTGTGGCTTGAAGGGGAGCGGGACTGCCCAGCCACTACAATTGAg < 2:3790022/76‑1 (MQ=255) aaCTTGGGGATATACTGAAGGTAATGCTCTGTGTCTTGAATGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTa < 2:4912395/76‑1 (MQ=255) aCTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTaa < 1:1515298/76‑1 (MQ=255) cTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAg < 1:2559944/76‑1 (MQ=255) ttGGGGATATACTGAAGGTAATGCTCTGTGGCTTTAATGGGAGGTGGGGTGCACAGCCACTACAAATGAGGTAAGa > 2:2297517/1‑76 (MQ=255) tGGGGATATACTGAAGGTAATGCTTTGTGGTTTGAATGGGAGCGGGACTGCCCAGCCACTACAATTGAGTTAAGaa < 1:372601/76‑1 (MQ=255) tGGGGATATACTGAAGGTAATGCTCTGTGGCTTTTATGGGGGCTGGACTGCACAGCCCCTACAATTTAGGTTAGaa > 1:404622/1‑76 (MQ=255) tGGGGATATACTGAAGGTAATGCTCTGTGGCTTTAATGGGGGGTGGGCTGCACAGCCCCCACAATTGAGTTAAGaa > 2:3351731/1‑76 (MQ=255) ggggATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAAt < 1:1484034/76‑1 (MQ=255) ggATATACTGAAGGTAATGCTCTGGGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATaa < 1:1698609/76‑1 (MQ=255) gaatatACTGAAGGTAATGCTCTGTGGCTTGAAGGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATaa < 2:2795464/74‑1 (MQ=255) atatACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAac < 1:1163849/76‑1 (MQ=255) tataCTGAGGGTAATGCTCTGTGGCTTGAAGGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAaca < 1:3286254/76‑1 (MQ=255) tataCTGAAGGTAATGCTCTGTGGCTTTTATGGGGGCTGGGCTGCACCGCCCCTACAAATGTGTTTAGaaaaaaac > 1:2825362/1‑74 (MQ=255) ataCTGAAGGTAATGCTCTGTGGCTTTTATGGGGGGTGGACCGGCCAGCCCCCACACATGAGTTTAGAAAAAacaa > 1:278723/1‑75 (MQ=255) ataCTGAAGGTAATGCTCTGTGGCTTTTATGGGGGCTGGGCTGCACAGCCCCCACAAATGAGTTAAGAATAAacaa > 1:2921546/1‑75 (MQ=255) ataCTGAAGGTAATGCTCTGTGGCTTGAATGGGAGGTGGACTGCACAGCCACTACAAATTAATTAAGAATAAacac > 1:388342/1‑76 (MQ=255) taCTGAAGGTAATGCTCTGTGGCTTGAAGGGGAGCTGGACTCCCCAGCCACTACAATTGAGTTAAGAATAAACACt < 2:1154413/76‑1 (MQ=255) aCTGAAGGTAATGCTCTGTGGCTTTAATGGGGGCTGGGCTGCACAACCACCACAAATGAGTTGAAAAAAAAAACTc > 1:421414/1‑76 (MQ=255) tGAAGGTAATGCTCTGTGGCTTTAATGGGGGCTGGACTGCACCGCCACTACAATTGAGTTAAGAATAAACAATCAt > 2:1263374/1‑76 (MQ=255) gAAGGTAATGCTCTGTGGCTTTTATGGGGGGCGGGCTGCACAACCACCACAAATGAGTTTAAAAAAAAAACCCATg > 2:563516/1‑76 (MQ=255) aGGTAATGCTTTGTGGCTTGAATGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGga < 1:5188179/76‑1 (MQ=255) aGGTAATGCTCTGTGGCTTTTATGGGGGGTGGGCTGCACAGCCACTACAAATGAGTTAAGAAAAAAAACCCATGga > 2:3425470/1‑76 (MQ=255) ggTAATGTTCTGGGGTTTGAATGGGAGCTGTCCTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGgag < 1:463562/76‑1 (MQ=255) ggTAATGCTCTGTGGCTTGAAGGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGgag < 1:4465372/76‑1 (MQ=255) gTAATGCTCTGTGGCTTGAATGGGGGCTGGGCTGCACCGCCACCACAAATGAGGTAAGAATAAAAACTCATGggga > 2:3292510/1‑76 (MQ=255) tAATGCTCTGTGGCTTTAATGGGGGGTGGGCTGGCCAGCCACTACAATTTAGTTAAGAAAAAACACTCACGCAGAt > 2:563143/1‑76 (MQ=255) tAATGCTCTGTGGCTTTAATGGGGGCTGGACTGCACAGCCCCTACAAATGAGTTAAGAATAAACACTCATCGAGAt > 2:2218164/1‑76 (MQ=255) tAATGCTCTGTGGCTTTAAAGGGGGCTGGACTGCCCAGCCCCTACAATTTAGTTAAAAAAAAAAACCCATGGAGAt > 2:943664/1‑76 (MQ=255) tAATGCTCTGTGGCTTGAAGGGGAGCTGGCCCGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGAt < 2:3226450/76‑1 (MQ=255) aaTGCTCTGTGGCTTTTATGGGGGCTGGGCTGCACAGCCCCTACAATTTAGTTTAGAAAAAACACCCATGGAGAtt > 2:4032285/1‑76 (MQ=255) aaTGCTCTGTGGCTTTAATGGGAGGTGGGCTGCACAGCCACTACAATTTAGTTAAGAAAAAAAACTCCTGGAGAtt > 2:4912478/1‑76 (MQ=255) aaTGCTCTGGGGCTTGAATGGGAGCGGGACTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGAtt < 2:832915/76‑1 (MQ=255) tgtctcgGTGGCTGGAATGGGAGCTGGCCCGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTa < 2:745867/73‑1 (MQ=255) tGCTCTGTGGCTTTAATGGGGGCTGGGATGCACAGCCACCACAATTTAGGTTAGAATAAAAACCACAGGAGATTAt > 2:2183968/1‑76 (MQ=255) gCTCTGTGGCTTTAATGGGGGGTGGGCTGCACAGCCACTACAAATGAGTTAAGAAAAAACAATCATGGAGATTATc > 2:4295092/1‑76 (MQ=255) tctGTGGCTTGAATGGGAGCTGGACTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGc < 1:2564418/76‑1 (MQ=255) ctGTGGCTTGAAGGGGAGCTGGCCTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGcc < 1:4416142/76‑1 (MQ=255) gtgGCTTGAATGGGAGCGGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCt < 1:3816773/76‑1 (MQ=255) gtgGCTTGAAGGGGATCTGCACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCt < 2:1672259/76‑1 (MQ=255) tgGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTg > 2:5021889/1‑76 (MQ=255) tgGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTg > 1:2023498/1‑76 (MQ=255) tgGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTg > 2:1580135/1‑76 (MQ=255) ggCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGc > 1:1253696/1‑76 (MQ=255) ggCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGc > 1:5068989/1‑76 (MQ=255) tttGAAGGGGAGCGGGACTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCAt < 2:2015228/75‑1 (MQ=255) cTTGAAGGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCAt < 1:943152/76‑1 (MQ=255) tGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATgg > 1:5032055/1‑76 (MQ=255) gAATGGGAGCTGGACTGCACAGCCACTACAATTGTGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATggg > 2:2319560/1‑76 (MQ=255) gAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATggg > 2:3196060/1‑76 (MQ=255) gAAGGGGAGCTGTACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCCCCCTGCATggg < 2:3526517/76‑1 (MQ=255) aaggggAGCGGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGa < 1:4118789/73‑1 (MQ=255) aTGGGAGCTGGACTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGaa < 2:2723106/76‑1 (MQ=255) aTGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCCTGGGaa > 1:362949/1‑76 (MQ=255) aTGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGaa > 2:3672589/1‑76 (MQ=255) aTGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGaa > 1:1959802/1‑76 (MQ=255) tGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGaaa > 2:3471744/1‑76 (MQ=255) ggggAGCTGGCCTGCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGaaa < 1:1019093/75‑1 (MQ=255) gggAGCGGGACTACACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAAc < 1:3429310/76‑1 (MQ=255) ggAGCTGGCCTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACa < 1:466656/76‑1 (MQ=255) ggAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACa > 1:2833964/1‑76 (MQ=255) ggAGCGGGCCTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACa < 2:2178861/76‑1 (MQ=255) gAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACAt < 1:664666/76‑1 (MQ=255) aGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACAtt > 1:778969/1‑76 (MQ=255) aGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACAtt > 2:2427549/1‑76 (MQ=255) aGCTGGACCCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACAtt < 1:2883482/76‑1 (MQ=255) cTGGACTCCCCAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGa < 2:3904921/76‑1 (MQ=255) tGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGaa > 1:1725515/1‑76 (MQ=255) tGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGaa > 1:3756298/1‑76 (MQ=255) tGGACTGCACAACCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGaa > 1:4948285/1‑76 (MQ=255) ggACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGaaa > 2:4776775/1‑76 (MQ=255) ggACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGaaa > 2:4201003/1‑76 (MQ=255) aCTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGt > 2:1785851/1‑76 (MQ=255) cTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTg > 2:1639379/1‑76 (MQ=255) gCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTGAt > 2:3096126/1‑76 (MQ=255) | TGTTTGAAGATGATGGCGAAACTTGGGGATATACTGAAGGTAATGCTCTGTGGCTTGAATGGGAGCTGGACTGCACAGCCACTACAATTGAGTTAAGAATAAACACTCATGGAGATTATCGCCCTGCATGGGAAACATTGAAAGTGATA > REL606/3783491‑3783639 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |