Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,519,907 | 0 | C | G | 35.6% | 280.9 / 123.6 | 146 | S91S (TCG→TCC) | yddL | predicted lipoprotein |
Reads supporting (aligned to +/- strand): ref base C (87/7); new base G (0/52); total (87/59) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.55e-32 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 198 total aligned reads displayed.
TGGGATGATTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGATCATTGATT > REL606/1519832‑1519982 | tGGGATGATTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCcgcgc < 2:1226641/76‑1 (MQ=255) tGGGATGATTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCcgcgc > 2:5364388/1‑76 (MQ=255) tGGGATGATTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCcgcgc > 1:5158622/1‑76 (MQ=255) gggATGATTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCcgcgcg > 2:4163011/1‑76 (MQ=255) gggATGATTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCcgcgcg > 1:3379066/1‑76 (MQ=255) gggATGATTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCcgcgcg > 1:2271377/1‑76 (MQ=255) gggATGATTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCcgcgcg > 2:3961962/1‑76 (MQ=255) ggATGATTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGa < 1:3249111/76‑1 (MQ=255) atgatTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGc < 2:5867678/76‑1 (MQ=255) tgatTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGcc < 2:2848333/76‑1 (MQ=255) gatTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCg > 2:1825938/1‑76 (MQ=255) gatTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCg > 2:546094/1‑76 (MQ=255) gatTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCg > 2:6223653/1‑76 (MQ=255) tgtgTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTg > 1:1542695/1‑76 (MQ=255) tAATATAAGGTGAGTGATCATATTGTGATGGGCACTTAACTAAATCAGAGTTTCTTGCCGGGGGAGCCGTTGGCtt < 2:2572539/76‑1 (MQ=255) tAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGGGAGCCGTTGGCtt < 2:6352280/76‑1 (MQ=255) tAATATAAGGTCAGTGATAATATTGTGAGGGGCCCTTAACTCAATCAGAGTTTCTTGCCGCGGGAGCCGTTGGCtt < 2:4930403/76‑1 (MQ=255) aaTATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTc > 1:4896726/1‑76 (MQ=255) aaTATAAGGTCAGTGATCATATTGTGATGGGCACTTAACCAAATAAGAGTTTCTTGCCGGGGGAGCCGTTGGCTTc < 2:3789243/76‑1 (MQ=255) ataAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGCGATTCTTGCCGCGCGAGCCGTTGGCTTCATc > 1:1518796/1‑76 (MQ=255) ataAGGTCAGTGATCATATGGTGAGGGGCACTTAACTCAATAAGAGTTTCTTGCCGCGGGAGCCGTTGGCTTCATc < 2:2906833/76‑1 (MQ=255) aaaGGTCAGTGTTCATATTGTGAGGGCCATTTAACTAAATCAGAGATTCTGCCCGCGGGAGCCGTTGGCTTCAtag < 2:3539025/75‑3 (MQ=255) agttCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGGGAGCCGTTGGCTTCATCGtt < 2:2054034/73‑1 (MQ=255) aGGTCAGTGATCATATTTTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGGGAGCCGTTGGCTTCATCGtt < 1:5812753/76‑1 (MQ=255) aGGTCAGTGATAATTTTGTGAGGGGCATTTAACTCAATCAGAGTTTCTTGCCGGGGGAGCCGTTGGCTTCATCGtt < 1:4954077/76‑1 (MQ=255) gTCATTGATCATTTTGTGATGGGCACTTAACTCAATCAGAGTTTCTTGCCGGGGGAGCCGTTGGCTTCATCGTtat < 1:4524246/76‑1 (MQ=255) gTCAGTGATCATTTTGTGATGGGCATTAAACTAAATCAGAGTTTCTGCCCGCGGGAGCCGTTGGCTTCATCGTtat < 2:4009671/76‑1 (MQ=255) gTCAGATATCATATTGTGATGGGCACTTAACACAATCAGATTTTCTGGCCGCGGGAGCCGTTGGCTTCATCGTtat < 2:4627196/76‑1 (MQ=255) tCAGTGATCTGATTGTGAGGGGCACTAAACACAATCAGAGATTCTTGCCGGGGGAGCCGTTGGCTTCATCGTtata < 2:5849989/76‑1 (MQ=255) tCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTtata > 1:2736748/1‑76 (MQ=255) tCAGTGATCATATTGTGAGGGGCACTTAACTCAATCAGAGTTTCTTGCCGCGGGAGCCGTTGGCTTCATCGTtata < 1:1470870/76‑1 (MQ=255) taGTGATCATATTGTGATGGGCACTTACCTAAATCAGAGTTTCGTGCCGCGGGAGCCGTTGGCTTCATCGTTATAg < 2:835688/75‑1 (MQ=255) aGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGc > 1:4294284/1‑76 (MQ=255) gtaaTCATTTGGTGATGGGCATTTAACTCAATCAGAGATTCTTGCCGCGGGAGCCGTTGGCTTCATCGTTATAGcc < 2:1581278/73‑1 (MQ=255) tGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCt > 1:2773044/1‑76 (MQ=255) ggATCTTTTTGTGATGGGCCTTTAACTCAAACGGAGTTTCTTGCCGGGGGAGCCGTTTGCTTCATCGTTATAGCCt < 2:4877843/75‑1 (MQ=255) gATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCtt > 1:1164509/1‑76 (MQ=255) gATCATATTGTGAGGGGCACTAAACTAAATCAGAGTTTCTTGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCtt < 2:6304586/76‑1 (MQ=255) aTCATATTGGGAGGGGCATTTAACTCAATAAGATTTTCTTGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCttt < 2:4738899/76‑1 (MQ=255) cATATTTGGATGGTCACTTAACTCAATCAGAGTTTTGTCCGGGGGGAGCCGTTGGCTTCATCGTTATAGCCTTTaa < 2:5022296/76‑1 (MQ=255) cATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTaa > 1:105875/1‑76 (MQ=255) cATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTaa > 1:2499362/1‑76 (MQ=255) atatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTaaa > 1:6244113/1‑76 (MQ=255) atatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTaaa > 1:3678556/1‑76 (MQ=255) atatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTaaa > 2:1767197/1‑76 (MQ=255) tatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTaaaa > 1:3014890/1‑76 (MQ=255) tatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTaaaa > 2:685922/1‑76 (MQ=255) ttggtgAGGGGCACTTACCTAAATCAGAGTTTCTTGCCGGGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAAt < 2:6026729/73‑1 (MQ=255) atTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAAt > 1:6057227/1‑76 (MQ=255) atTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAAt > 1:3487775/1‑76 (MQ=255) atTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAAt > 1:2056831/1‑76 (MQ=255) ttGTGATGGGCACTTAACTCAAAAAGAGTTTCTTGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATc < 1:3146036/76‑1 (MQ=255) gtgATGGGCCCTTAACCAAATCAGAGTTTCTGGCCGGGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCat < 1:1722912/76‑1 (MQ=255) gtgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCat > 2:2343244/1‑76 (MQ=255) gtgATGGGCACTTAACTAAATCAGAGATTCTTGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCat < 2:2320823/76‑1 (MQ=255) tGGCCACTTACCTAAATCAGAGTTTCTGCCCGGGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTc < 1:2857700/76‑1 (MQ=255) gccACTTACCTCAATCAGATTTTCTTGCGGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTccc < 2:3205502/74‑1 (MQ=255) gCATTTAACTCAATAAAAGTTTCGTGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCa < 2:6249409/76‑1 (MQ=255) gCACTAAACAAAATCAGATTTTCTTGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCa < 2:780044/76‑1 (MQ=255) gCAATTAACACAATCAGAGTTTCTGGCCGGGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCa < 2:57593/76‑1 (MQ=255) cACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAg > 2:2710921/1‑76 (MQ=255) cACTTAACTAAAAAAGAGTTTCTTGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAg < 1:2552864/76‑1 (MQ=255) cACTTAACAAAATAAGATGTTCTTGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAg < 2:2270957/76‑1 (MQ=255) atttAACAAAATCAGAGTTTCTTGCGGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGt < 2:2514947/74‑1 (MQ=255) aCTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGt > 1:5336570/1‑76 (MQ=255) aCTTAACTCAATCAGAAATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGt > 2:2944407/1‑76 (MQ=255) ttAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGa > 2:5360821/1‑76 (MQ=255) ttAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGa > 2:6167588/1‑76 (MQ=255) aCTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCa > 2:3320285/1‑76 (MQ=255) aCTCAATCAAAGGTTCTCGCCGGGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCa < 1:6340699/76‑1 (MQ=255) aCTCAAACAGAGTTTCGTCCCGGGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCa < 1:1291742/76‑1 (MQ=255) cacAAACAGTGATTCTCGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCaa < 2:1029329/74‑1 (MQ=255) cTCAATAAGAGTTTCTCGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCaa < 1:3967959/76‑1 (MQ=255) tCAAGCAGAGTTTCTGGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCaaa < 2:2075949/76‑1 (MQ=255) cAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCaaaa > 2:104557/1‑76 (MQ=255) cAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCaaaa > 1:6134867/1‑76 (MQ=255) cAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCaaaa > 2:5154033/1‑76 (MQ=255) aaTCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAAc > 1:5755125/1‑76 (MQ=255) tCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCa > 2:4605197/1‑76 (MQ=255) cAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCaa > 1:5619875/1‑76 (MQ=255) cAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCaa > 2:5115767/1‑76 (MQ=255) cAGAGATTCTGGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCaa < 2:5678146/76‑1 (MQ=255) aagagTTTCTGGCCGGGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCaa < 2:909869/75‑1 (MQ=255) agagATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAAt > 1:2383261/1‑76 (MQ=255) gagtTTCTGGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATc < 2:268841/76‑1 (MQ=255) gagaTTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATc > 2:1488390/1‑76 (MQ=255) gagaTTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATc > 1:2176897/1‑76 (MQ=255) agaTTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCa > 2:3134840/1‑76 (MQ=255) agaTTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCa > 1:392713/1‑76 (MQ=255) aTTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTTACCAAAACCAATCAtt > 2:2395729/1‑76 (MQ=255) aTTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCAtt > 1:653513/1‑76 (MQ=255) ttCTGGCCGGGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCAttt < 1:599717/76‑1 (MQ=255) tctcGCCGGGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTg < 2:3488550/76‑1 (MQ=255) tCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTg > 2:6211404/1‑76 (MQ=255) tCTGGCCGCGGGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTg < 1:5513956/76‑1 (MQ=255) ttGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGAt > 2:426470/1‑76 (MQ=255) tGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGATc > 2:4337153/1‑76 (MQ=255) gcggcgggAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGATCa < 2:3661905/73‑1 (MQ=255) gCCGCGCGAGTCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGATCa < 2:546944/76‑1 (MQ=255) cgcgcgAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGATCAtt > 1:2317634/1‑76 (MQ=255) | TGGGATGATTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGATCATTGATT > REL606/1519832‑1519982 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |