Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 273,581 | 0 | T | G | 37.5% | 80.9 / 17.0 | 56 | G136G (GGT→GGG) | yagU | conserved inner membrane protein |
Reads supporting (aligned to +/- strand): ref base T (15/20); new base G (21/0); total (36/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.54e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTCC > REL606/273508‑273655 | tgtTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGc < 2:1026704/76‑1 (MQ=255) ttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGGttt > 2:1371015/1‑76 (MQ=255) gCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg < 2:976508/76‑1 (MQ=255) gCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg < 2:723418/76‑1 (MQ=255) tGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGcc > 1:1088028/1‑76 (MQ=255) gTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGccc > 1:101861/1‑76 (MQ=255) gTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCCTTTGccc > 2:503088/1‑76 (MQ=255) tCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCa < 1:817692/76‑1 (MQ=255) tCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCa < 2:213590/76‑1 (MQ=255) cGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCaa < 2:559780/76‑1 (MQ=255) cGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCaa > 1:680860/1‑76 (MQ=255) cGGTTATTGTGTGGTCGCTGAAGTATTCCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTTCCCaa > 2:435875/1‑76 (MQ=255) ttATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACtt > 2:1034743/1‑76 (MQ=255) ttATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCCAcct > 2:950181/1‑74 (MQ=255) tATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACttt > 1:195801/1‑76 (MQ=255) aTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCCTTAGCCCAACtttt > 1:1527400/1‑76 (MQ=255) ttGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCCACttttt > 2:221177/1‑76 (MQ=255) tgtgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGGTTTTGCCCAACtttttt > 2:928432/1‑75 (MQ=255) tgtgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACtttttt > 1:1117242/1‑75 (MQ=255) tgtgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACcttttt > 1:1046374/1‑75 (MQ=255) gtgtgGTCGCTGAAGTATTTCCAAAACTTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACttttttt > 1:988300/1‑74 (MQ=255) gtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCCACttttttttt > 2:652757/1‑75 (MQ=255) gtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGGTTTTGCCCACCTTTTTGttt > 2:598180/1‑75 (MQ=255) gtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTGttt > 2:706329/1‑75 (MQ=255) tgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTGttta > 1:565777/1‑74 (MQ=255) ggTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACTTTTTGTTTat > 2:1159821/1‑76 (MQ=255) ggTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTGTTCat > 1:197529/1‑76 (MQ=255) cGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTCGCCCAACTTTTTGTTCATATg > 2:1401129/1‑76 (MQ=255) gCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTTTTCATATGa > 1:333162/1‑76 (MQ=255) tGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACCTTTTTTTCATATGAtt > 1:449873/1‑76 (MQ=255) gAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTTTTCATATGAttt > 2:1354357/1‑76 (MQ=255) gAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAttt < 2:146551/76‑1 (MQ=255) gAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCCACTTTTTTTTTATATGAttt > 1:830495/1‑76 (MQ=255) aaGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttc < 1:984852/76‑1 (MQ=255) aGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttct < 2:870161/76‑1 (MQ=255) gTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttctt < 1:1467874/76‑1 (MQ=255) aTTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTTTTCATATGAtttctttc > 1:96969/1‑76 (MQ=255) tttCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTTTTTTTTTGAttttttttc > 2:805108/1‑74 (MQ=255) ccAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCctc > 1:1406107/1‑76 (MQ=255) cAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTTTTCATATGATTTCTTTTCctct > 2:190044/1‑76 (MQ=255) aaaTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCAACTTTTTTTTTATATGATTTTTTTTCCtctcct > 1:883179/1‑74 (MQ=255) aaaTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCAt < 2:924101/76‑1 (MQ=255) aaaTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCAt < 2:1498905/76‑1 (MQ=255) aaTTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTCATATGATTTTTTTTCCTCTCCTg > 2:950338/1‑76 (MQ=255) aaaCTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTGTTTATATGATTTTTTTTCCTCTCATggggc > 2:777201/1‑74 (MQ=255) aaaCTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTCATATGATTTTTTTTCCTCTCATGGGAc > 1:1341734/1‑76 (MQ=255) aaCTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACTTTTTGTTCATATGATTTCTTTTCCTCTCATGGGACt > 1:1273043/1‑76 (MQ=255) ctctGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTTTTTTCCTCTCATGGGACTGa > 1:209532/1‑76 (MQ=255) ctctGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGa < 1:324241/76‑1 (MQ=255) ctctGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGa < 2:827392/76‑1 (MQ=255) tctGGCAGGGCTTACTGGCAGGGGCTTTAGCCCCACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGac > 1:1484944/1‑76 (MQ=255) ctGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTTATATGATTTTTTTTCCTCTCATGGGACTGaca > 2:1034615/1‑76 (MQ=255) ctGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTGTTCATATGATTTTTTTCCCTCTCATGGGGCTGaca > 1:306754/1‑76 (MQ=255) ggCAGGGCTTACTGGCAGGTGCTTTAGCCGAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGAcacc < 2:447161/76‑1 (MQ=255) ggCAGGGCTTACTGGCAGGTGCTTTAGCCCAACCTTTTTTTCATATGAATTCTTTCCCCCCCCCGGGACCGAcacc > 1:987324/1‑76 (MQ=255) aGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTTTTCTTATTATTTCTTTTCCTCTCCTGGGACTGACaccacc > 1:375090/1‑76 (MQ=255) gggCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTCATATGATTTCTTTTCCTCTCCTGGGACTGACaccaccc > 2:291387/1‑75 (MQ=255) cTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTg > 2:359927/1‑76 (MQ=255) ttACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGt < 2:1149257/76‑1 (MQ=255) tACTGGCAGGTGGTTTTGCCCAACTTTTTTTTTATATGATTTTTTTTCCTCTCATGGGACTGACACCACCTCTGtt > 2:914178/1‑76 (MQ=255) cTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTg > 1:642710/1‑76 (MQ=255) cTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTg < 2:404007/76‑1 (MQ=255) tGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGa < 2:90387/76‑1 (MQ=255) ggCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGAt < 2:1038961/76‑1 (MQ=255) cAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCt > 2:1047989/1‑76 (MQ=255) aGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCtt < 2:389418/76‑1 (MQ=255) ggTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTc > 2:491586/1‑76 (MQ=255) gTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTcc < 2:446600/76‑1 (MQ=255) | TGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTCC > REL606/273508‑273655 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |