Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 2,087,670 | C→T | 11.9% | intergenic (+483/+435) | PP_1865 → / ← aroH | transposase/Phospho‑2‑dehydro‑3‑deoxyheptonate aldolase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 2,087,670 | 0 | C | T | 11.9% | 108.2 / 9.0 | 42 | intergenic (+483/+435) | PP_1865/aroH | transposase/Phospho‑2‑dehydro‑3‑deoxyheptonate aldolase |
Reads supporting (aligned to +/- strand): ref base C (23/14); new base T (3/2); total (26/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.93e-01 |
GCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGC > AE015451/2087447‑2087879 | gCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCg > 1:773134/1‑236 (MQ=255) cAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCtt > 1:737008/1‑226 (MQ=255) cAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCtt < 2:737005/226‑1 (MQ=255) gAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAg > 1:214844/1‑236 (MQ=255) gAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGA‑GGGAGATGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGcc > 2:560856/1‑233 (MQ=21) aTCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAgcg > 1:649036/1‑150 (MQ=255) ggATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGA‑GGGAGATGACGGGCATGGATGCCCGTCAAgcg > 2:644623/1‑191 (MQ=21) ggATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGA‑GGGAGATGACGGGCATGGATGCCCGTCAAgcg < 1:644626/191‑1 (MQ=21) aGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGcc > 2:423558/1‑234 (MQ=255) ggATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGc > 2:159099/1‑146 (MQ=255) gCCGGCATGGCTGCCCTCGGCCTGGTGTAACTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCtt < 2:214416/141‑1 (MQ=255) ccGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATggg > 1:447854/1‑194 (MQ=255) cGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGc < 1:755774/182‑1 (MQ=255) cGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGc > 2:755771/1‑182 (MQ=255) gCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGAAGCCCGTCAAGCGCTG‑CGCCCCA‑g > 1:822229/1‑185 (MQ=255) tGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGc > 1:685424/1‑190 (MQ=255) ggCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGc > 1:652408/1‑176 (MQ=255) cTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGa < 1:423561/236‑1 (MQ=255) cccTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGGTGCCCGTCAAGCGCTG‑GGCCCCA‑g > 1:351577/1‑176 (MQ=255) tGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTg > 1:229964/1‑178 (MQ=255) tgtACCGGGAAGCACGCCGTTGCTCTTGCTCTTGCTCGTGCTTTGGCTTTGGCTTCTGGGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGccc < 2:21831/145‑1 (MQ=17) ttgatgaaCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGACGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTg < 1:308762/106‑1 (MQ=255) tGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGGAGGGGCCGgat > 1:248774/1‑235 (MQ=255) gaagaaCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTg < 2:773131/200‑1 (MQ=255) ttgctcttgctcttgctTTGGCTTTTGCTTCTAGGCGCGCGGTAGTTCAGGCGACGAAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGAt < 1:307536/142‑1 (MQ=255) ttgctcttgctcttgctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGa < 1:388688/170‑1 (MQ=37) cttgctcttgctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGc > 2:388685/1‑143 (MQ=18) gctcttgctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCgg < 2:795416/87‑1 (MQ=35) gctcttgctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCgg > 1:795419/1‑87 (MQ=35) tcttgctTTGGCTTTTGCTTCTAAGCGTGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTc > 1:574848/1‑113 (MQ=17) tcttgctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTg > 1:603821/1‑235 (MQ=34) cttgctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGCGt > 1:870569/1‑171 (MQ=14) cttgctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCACGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGc > 1:573484/1‑236 (MQ=21) cttgctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGACTGCGACTTCAGGAGGCCGAGCGCAGGGCTTg > 1:380017/1‑81 (MQ=18) cTTTTGCTTCTAGGCGCGCGGTAGTTCTGGCGACGCAGATTGCGACTGCAGGAGGCCGAGCGCAGGGCGTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAgcg < 2:176473/107‑1 (MQ=11) ttCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTg < 2:927401/122‑1 (MQ=32) tCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCAt > 2:663180/1‑83 (MQ=33) tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTg < 1:40947/61‑1 (MQ=35) tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTg > 2:40947/1‑61 (MQ=35) aGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGA‑GGGAGATGACGGGCATCCATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGtt > 2:529768/1‑176 (MQ=11) aGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGA‑GGGAGATGACGGGCATCCATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGtt < 1:529771/176‑1 (MQ=11) aGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGgc > 1:555602/1‑234 (MQ=32) aGCGCAGGGCTGGCGGA‑GGGAGGTGACGGGCATGGATGCCCGGCAAGCGCTG‑GGCCCCA‑GGAGGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGt < 1:25740/180‑1 (MQ=16) | GCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGA‑GGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGC > AE015451/2087447‑2087879 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |