Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,207,576 | A→G | 12.4% | intergenic (‑44/‑571) | purH ← / → rrsE | bifunctional AICAR transformylase/IMP cyclohydrolase/16S ribosomal RNA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,207,576 | 0 | A | G | 12.4% | 83.8 / 3.7 | 39 | intergenic (‑44/‑571) | purH/rrsE | bifunctional AICAR transformylase/IMP cyclohydrolase/16S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base A (16/18); new base G (3/2); total (19/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.61e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.25e-01 |
CACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCG‑CGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAA > NC_000913/4207439‑4207718 | cgccGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCG‑CGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAAttttt < 1:377740/141‑1 (MQ=255) cGGCAAGTGCCTGGGCCCATTCGACGCTAACAGCATTGTCTGCCAAACTAAGCAGCGCGACGTCGTA‑TAGCAGATCTTGTTGCATGGTAAATCCCCAAGATTCAGCTATTTCACTTAGCCTTAGCTGAATTTTTCGCGaaaac < 2:8916/143‑2 (MQ=255) gggCGCATTCGACGATACCGGCTTTGCCAGAAACACTCAGCAGAGCG‑CGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGc < 2:297599/143‑1 (MQ=255) gaTACCGGCTTTGTGAGAAACACTGAGCAGAGCG‑CGGGGGGCTGGACGACGTTGTTGCATGGTCACTCCCCTGGAGTTGACTATTACAGAGAGCGTTAGCTGAATGTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCa < 2:255042/143‑1 (MQ=255) cGGCTTTGTCAGAAACACTGAGCAGAGCG‑CGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCtt > 1:56661/1‑143 (MQ=255) cGGCTTTGTCAGAAACACTGAGCAGAGCG‑CGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCtt < 1:370322/143‑1 (MQ=255) tgCTTTGTCAGAAACACTGAGCAGAGCG‑CGGCGGACTGGACGACGTTGGTGCATGGTAAAGCCCCTGGATTTGACTATTACAGAGGGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTa < 1:81572/142‑1 (MQ=255) ggCTTTGTCAGAAACACTGAGCAGAGCG‑CGACGGACTGGACGACGTAGTTGCATGGTAAATCCCCTGGATTTGAATATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACGCAGCTAACGACACTAACGGGGCAACCTTa > 2:164624/1‑143 (MQ=255) ggCTGTGTCAGAAACACTGAGCAGAGCG‑CGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTa < 2:241404/143‑1 (MQ=255) acTGAGCAGAGCG‑CGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTCTTTTTCGCCCGCAtt > 2:310349/1‑143 (MQ=255) acTGAGCAGAGCG‑CGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACCGGGCATCCTTCTTTTTCGCCCGCCtt > 1:459301/1‑143 (MQ=255) agagCG‑CGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGa < 1:141340/143‑1 (MQ=255) gagCG‑CGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTCTCGCCCGCATTGTAACGaa < 2:253211/143‑1 (MQ=255) cg‑cgGCGGACTGGACGGCTTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACGCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAAc < 1:456626/143‑1 (MQ=255) cg‑cgGCGGACTGGACGACGTTGGTGCATGGTACATCCCCTGGATTTGACTACTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAAc < 1:190716/143‑1 (MQ=255) gcggACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCCTCCTTATTGTTCGCCCGCAATGTCACGGACACGttt > 2:437033/1‑143 (MQ=255) acgacgTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAGCGTTTGCGCAACGc < 1:109402/143‑1 (MQ=255) acgTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCCTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCCACGCTcg > 1:269645/1‑143 (MQ=255) cATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCACTAAAGGGGCATCCTTATTTTTAGCCCGCATTGTAACGAAGACGATTGCGCAAAGCTGGGGAAttttt > 2:318436/1‑143 (MQ=255) aTGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGTCATCCTTTTTTTTCTCCCGCATTTTAACGACAACGTATTCGCTACTCTCACGAATTTTtc > 1:95733/1‑143 (MQ=255) tGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGCTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTtct > 1:303955/1‑143 (MQ=255) tagaTCCCCTTGGTTGGACTGTTGCAGAGAGCGTTAGCTGAGTTTTTCGCGAAAAACTCAGCTACCCCCCCTAACGGGGCATCCTTATTTTTAGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTAtt < 2:451084/140‑1 (MQ=255) aTCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCACCGCTCGCGAATTTTTCTCTTTCa > 2:262457/1‑143 (MQ=255) aTCCCCTGGATTTGACTATTACAGAAAGGGTTAGCTGAATTTTTCGCGAAAGACTCAGCTACAGCGCTTATCGTGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCa < 1:65664/143‑1 (MQ=255) ccccTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAAt < 1:132332/143‑1 (MQ=255) gATTTGACTATTCCAGAGAGCGTTAGCTGAATTTTTCGCTAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGAt < 2:368562/143‑1 (MQ=255) gATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGAt < 1:222035/143‑1 (MQ=255) ttGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGACCCTAACGGGGCATACTTATTTTTCGCACGCATTGTAAATAAACCATTTGCTCAATGCACGCGAATTTTTCTCTTTCAATgatgatgcg > 2:321366/1‑140 (MQ=255) aTTACGGAGTGCGTTAGCTGAAGTTTTCGCTAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAAtttt < 2:69526/143‑1 (MQ=255) ttACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATCGTCACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGCTGAACACAAATTTg > 2:373855/1‑143 (MQ=255) aCAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTCAAAAACGCTGGCGGATTTTTCTATTTCCATGGTGATCACAATTTTGAc > 2:406411/1‑143 (MQ=255) gagagCGTGAGGTGAATTTTGCGCGAAACACTTAGCTATCGCCCCTAACGGGGCATCCGTCTTTTTCGCCCGCATTGTAGCGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACtgt < 2:43481/143‑1 (MQ=255) gTTAGCTGAGTTTTTCGCGAAAAATTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTAc > 2:104869/1‑143 (MQ=255) ttAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTGcc > 1:261893/1‑143 (MQ=255) cTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGAGCACAATTTTGACTGTGGTTAACgggg > 2:263858/1‑143 (MQ=255) tGAGTTTTTCGCGAAAAATTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACtgtg < 2:134885/132‑1 (MQ=255) tGAGTTTTTCGCGAAAAATTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACtgtg > 1:134885/1‑132 (MQ=255) tGAGTTTTTCGCGAAAAATTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTggg > 1:137282/1‑143 (MQ=255) tGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGTAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGGTCACAATGTTGACTGTGGTTGCCGTggg > 1:325561/1‑143 (MQ=255) aTTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACTGTGGGCaa > 1:117170/1‑143 (MQ=255) | CACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCG‑CGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAA > NC_000913/4207439‑4207718 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |