Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 273,581 | 0 | T | G | 31.5% | 146.2 / 19.7 | 73 | G136G (GGT→GGG) | yagU | conserved inner membrane protein |
Reads supporting (aligned to +/- strand): ref base T (14/36); new base G (23/0); total (37/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.48e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTCCG > REL606/273506‑273656 | aGTGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGt < 2:878870/76‑1 (MQ=255) gtgtTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTg < 2:1468337/76‑1 (MQ=255) tgtTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGc > 2:1092214/1‑76 (MQ=255) ttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCttt < 2:1183938/76‑1 (MQ=255) ttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCgtt > 1:908703/1‑76 (MQ=255) tGCTGTCGGTTATTCTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTa < 2:622278/76‑1 (MQ=255) gCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg < 1:923380/76‑1 (MQ=255) gCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg < 2:1511647/76‑1 (MQ=255) gCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg < 2:1265638/76‑1 (MQ=255) gCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCgtttg > 2:370833/1‑74 (MQ=255) cTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGGtttttc > 2:790744/1‑73 (MQ=255) cTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGc < 1:1148935/76‑1 (MQ=255) tGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGGTTTTGcc > 2:1271900/1‑76 (MQ=255) tCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCa < 1:883054/76‑1 (MQ=255) tCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCa > 1:808844/1‑76 (MQ=255) cGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCaa < 2:893285/76‑1 (MQ=255) cGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCaa > 2:637797/1‑76 (MQ=255) ggTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCAc > 2:287992/1‑76 (MQ=255) ggTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCAc > 2:947207/1‑76 (MQ=255) ttATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACtt > 1:995024/1‑76 (MQ=255) aTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACcttt > 1:287023/1‑76 (MQ=255) ttGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACttttt > 1:614092/1‑76 (MQ=255) tgtgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACtttttt > 2:1551012/1‑75 (MQ=255) tgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTGtt > 1:226606/1‑76 (MQ=255) tgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTGtt > 2:378888/1‑76 (MQ=255) gtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCACCCTttttttt > 2:1149387/1‑75 (MQ=255) gtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACttttttttt > 2:586483/1‑75 (MQ=255) tgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTCa > 2:441176/1‑76 (MQ=255) ggTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACttttttttttt > 1:996446/1‑73 (MQ=255) gTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCCACTTTTTTTTtata > 1:337007/1‑76 (MQ=255) cGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTTGCAGGTGCTTTAGCCCAACTTTTTGTTCATATg < 2:1460640/76‑1 (MQ=255) cGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATg < 1:342177/76‑1 (MQ=255) ctcaaGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTTATATGAt > 1:278133/4‑76 (MQ=255) gAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTCATATGAttt > 2:316026/1‑76 (MQ=255) aaGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttc > 1:202632/1‑76 (MQ=255) aaGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttc < 2:663422/76‑1 (MQ=255) tATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTGATATGAtttcttt < 1:704367/76‑1 (MQ=255) tATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttcttt < 2:514330/76‑1 (MQ=255) tATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttcttt < 2:1223563/76‑1 (MQ=255) tATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttcttt < 2:889149/76‑1 (MQ=255) aTTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCAACTTTTTTTTCATATGAtttctttc > 2:1631967/1‑76 (MQ=255) tttCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTcc < 2:1142670/76‑1 (MQ=255) ccAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTTATATGATTTCTTTTCctc > 1:19406/1‑76 (MQ=255) cAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCctct < 1:673724/76‑1 (MQ=255) aaaaaTTAAAATCTTGCAGGGCTTACTGGCAGGGGGGTTTGCCCCACTCTTTTTTTTTTTTATTTTTTTTCCtctc > 2:1581301/1‑76 (MQ=255) aaaaTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCa < 1:896078/76‑1 (MQ=255) aaaTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCAt < 2:1387989/76‑1 (MQ=255) aaTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATg < 1:1635763/76‑1 (MQ=255) aaTTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTTTTTATATTATTTTTTTTCCTCTCATg > 2:1125493/1‑76 (MQ=255) aaTTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTTATATGATTTTTTTTCCTCTCATg > 2:1505449/1‑76 (MQ=255) aaTTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTCATATGATTTCTTTCCCTCTCCTg > 1:227753/1‑76 (MQ=255) aTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATgg < 2:1633191/76‑1 (MQ=255) aTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATgg < 2:1595735/76‑1 (MQ=255) aTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATgg < 1:1508857/76‑1 (MQ=255) ttAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTTTTCATATGATTTCTTTCCCTCTCATggg > 1:184816/1‑76 (MQ=255) aaaCTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACTTTTTTTTTATATGATTTTTTTTCCTCTCCTggggc > 1:328610/1‑74 (MQ=255) aaaCTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACTTTTTTTTTATATGATTTTTTTTCCTCTCATGGGAc > 2:1279502/1‑76 (MQ=255) aaaCTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTTATATGATTTTTTTTCCTCTCATggggc > 1:436204/1‑74 (MQ=255) aCTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTg < 1:116912/76‑1 (MQ=255) aCTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTg < 1:1166874/76‑1 (MQ=255) ctctGGCAGGGCTTACTGGCAGGTGCTTTAGGCTCACCCTTTGTTTATATTGTTTCTTTTCCTCTCCTTGGAccgc > 1:660488/1‑73 (MQ=255) ctctGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGa < 2:704541/76‑1 (MQ=255) ctctGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTTTTTATATTATTTTTTTTCCTCTCCTTGGGCTGa > 1:366311/1‑76 (MQ=255) ctctGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTTTTTATATGATTTTTTTTCCCCTCCTGGGACTGa > 1:1332176/1‑76 (MQ=255) ctctGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTTATATGATTTTTTTTCCTCTCATTGGGCTGa > 1:660472/1‑76 (MQ=255) ctGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTGTTCATATGATTTTTTTCCCTCTCTTGGGGCTGaca > 2:1099253/1‑76 (MQ=255) ggCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGAcacc < 1:106843/76‑1 (MQ=255) gCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGAcacca < 1:441176/76‑1 (MQ=255) aGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACaccccc > 2:1554830/1‑76 (MQ=255) aGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACaccacc < 2:1206428/76‑1 (MQ=255) gggCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTTATATGATTTTTTTTCCTCTCATGGGGCTGACaccaccc > 2:710908/1‑75 (MQ=255) gggCTTACTGGCAGGAGCTTTTGACCCACTTTTTGTTTCTTTTTTTTTTTTCCCTCTCTTGGGAGCGACacccccc > 2:1146686/1‑75 (MQ=255) gCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCGCTCATGGGACTGACACCACctct < 1:41823/76‑1 (MQ=255) gCTTACTGGCAGGGGCTTTAGCCCCACTTTTTTTTTATATGATTTCTTTCCCTCTCATGGGACTGACACCACctct > 2:129756/1‑76 (MQ=255) cTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCTCCTCTg < 2:1282332/76‑1 (MQ=255) ttACTGGCAGGGGCTTTTGCCCAACTTTTTTTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCCCTGt > 2:1280733/1‑76 (MQ=255) aCTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGttt < 2:272197/76‑1 (MQ=255) aCTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGttt < 1:1643152/76‑1 (MQ=255) cTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTg > 2:105071/1‑76 (MQ=255) tGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGa > 2:1666219/1‑76 (MQ=255) ggCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGAt > 1:1343767/1‑76 (MQ=255) ggCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGAt < 2:1484679/76‑1 (MQ=255) ggCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGAt > 1:701063/1‑76 (MQ=255) gCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATc > 2:815949/1‑76 (MQ=255) cAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCt > 2:1423739/1‑76 (MQ=255) aGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCtt > 1:1671051/1‑76 (MQ=255) ggTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTc < 1:751540/76‑1 (MQ=255) tGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTCCg < 2:695162/76‑1 (MQ=255) | AGTGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTCCG > REL606/273506‑273656 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |