Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 248,447 | C→T | pseudogene (1624/1713 nt) | lfhA ← | flagellar system protein, promoterless fragment |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 248,447 | 0 | C | T | 96.5% | 87.5 / ‑6.3 | 29 | pseudogene (1624/1713 nt) | lfhA | flagellar system protein, promoterless fragment |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (20/8); total (21/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 |
TCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATAT > NC_000913/248319‑248563 | tCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATaa > 2:22894/1‑133 (MQ=255) cACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAACCACGGCGCACACAAGATTGtatacc > 2:313677/1‑131 (MQ=255) gccgcATCCGGCAGTTGTGGACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAAcgcg < 1:516057/133‑1 (MQ=255) cacATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGcc > 1:80805/3‑133 (MQ=255) aTCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAACCCCGGCGCAAACAACCTTGTATAACCCGCCAGc > 2:640468/1‑133 (MQ=255) tacGGCAGTTGTGCAACAATGCTGCATACTTATATCAACGCTCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCGTGTATAACGCGCCAGCa < 2:291664/131‑1 (MQ=255) aTGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAACCAACCTTGCATAACGCGCCAGCAAATGGCGCAGCTGCGGTg > 2:216505/1‑133 (MQ=255) gCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGc < 1:36973/133‑1 (MQ=255) gCTGGATACTTACATCAACGCCCCCATAATATTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGc < 1:313677/133‑1 (MQ=255) tACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCTAGCAAAGAGCGTAGATGCGGTAGCACcagc > 1:299128/1‑133 (MQ=255) tACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACcagc < 2:405755/133‑1 (MQ=255) cTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCTAGCAACGGGCGCAGCTGCGGTAGCACcagcag > 2:302295/1‑133 (MQ=255) aCATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCAc < 2:453940/133‑1 (MQ=255) aaCGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACcag < 1:379003/122‑1 (MQ=255) aaCGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACcag > 2:379003/1‑122 (MQ=255) cGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGCTGGTACTAGCAGCTCCGGGTCt > 1:296859/1‑133 (MQ=255) aTCTTCAACTCCAGCCCACCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATACCGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCccgccg < 1:250717/133‑1 (MQ=255) cAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCGGCAGCGGGTCTTTCCCAGCCGCTTTCATCTGc > 2:250717/1‑133 (MQ=255) gCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCAGCGGGTCTATCCCCGTCGCTTTAATCTGctc > 1:211959/1‑133 (MQ=255) tCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTATCATCTGCTCTTTCa > 1:313585/1‑133 (MQ=255) gCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATACCGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCCCCGGGTTTTCCCCCGCCGCTTTCATCTGCTCTTTCACCtg < 2:324964/133‑1 (MQ=255) tCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGAGGCa > 2:325287/1‑133 (MQ=255) ttATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCg < 2:80805/133‑1 (MQ=255) tataCGACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCCCCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGt < 1:22894/133‑1 (MQ=255) gACAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTActct > 1:325744/1‑133 (MQ=255) cAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTgg > 2:406268/1‑133 (MQ=255) cAGCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTgg > 2:450004/1‑133 (MQ=255) gCACATGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGaa > 2:518299/1‑133 (MQ=255) gCAAATCCAGCCGCGGCCCAAACACCCTTTTAAAGCGCGCCAAAAAAGGGCGCAGCCGGGGCGCCACCGCCACCCCCGGTTCTTCCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGaa < 2:477276/133‑1 (MQ=255) cacaTGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAAc > 1:122449/1‑133 (MQ=255) cacaTGCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAAc > 1:56054/1‑133 (MQ=255) tGCAGCCCCGTCGCACACACCCTTCTATCACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCCCCGGGTCTTCCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGt < 2:181029/133‑1 (MQ=255) tGCACCCCGGGCGCAAACACCCTTGTATACCCCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCACCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGt < 2:325744/133‑1 (MQ=255) gCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGtt > 2:451281/1‑133 (MQ=255) gCAGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGtt > 2:590413/1‑133 (MQ=255) aGCCCCGGCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGAATCTTATTCTGGAACTGGTTGa > 1:213473/1‑133 (MQ=255) ggCGCAAACAACCTTGTATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCatat > 1:396169/1‑133 (MQ=255) | TCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCATAT > NC_000913/248319‑248563 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |