Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,103,107 | T→G | 5.8% | pseudogene (283/446 nt) | wbbL ← | pseudogene, lipopolysaccharide biosynthesis protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,103,107 | 0 | T | G | 5.8% | 245.1 / 7.0 | 69 | pseudogene (283/446 nt) | wbbL | pseudogene, lipopolysaccharide biosynthesis protein |
Reads supporting (aligned to +/- strand): ref base T (33/32); new base G (2/2); total (35/34) | |||||||||||
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.85e-01 |
ACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTGG > NC_000913/2102965‑2103252 | acaATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatgat > 3:145851/1‑149 (MQ=255) tCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCg > 4:430320/1‑147 (MQ=255) tCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCg < 3:430320/147‑1 (MQ=255) aaaGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTc < 4:144947/149‑1 (MQ=255) aaaGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTc > 3:144947/1‑149 (MQ=255) aGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatga < 1:203602/128‑1 (MQ=255) aGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatga > 2:203602/1‑128 (MQ=255) aGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACa > 3:377405/1‑149 (MQ=255) ggAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatga < 1:941144/127‑1 (MQ=255) ggAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatga > 2:941144/1‑127 (MQ=255) gAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaaaa < 2:1062348/145‑1 (MQ=255) gAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaaaa > 1:1062348/1‑145 (MQ=255) ttCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa < 1:1060649/149‑1 (MQ=255) ttCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGACATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa < 4:145851/149‑1 (MQ=255) cTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatgat < 1:833902/119‑1 (MQ=255) cTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatgat > 2:833902/1‑119 (MQ=255) aataatCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgat < 2:1142402/115‑1 (MQ=255) aataatCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgat > 1:1142402/1‑115 (MQ=255) atCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCgg > 1:227608/1‑149 (MQ=255) atCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCgg < 4:377405/149‑1 (MQ=255) aTGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa > 1:443960/1‑133 (MQ=255) aTGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa < 2:443960/133‑1 (MQ=255) aTGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgatt > 3:18397/1‑108 (MQ=255) aTGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgatt < 4:18397/109‑2 (MQ=255) aaaGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCtata > 2:280100/1‑149 (MQ=255) aaGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTatat > 1:96153/1‑149 (MQ=255) gATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGt > 3:299926/1‑141 (MQ=255) gATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGt < 4:299926/141‑1 (MQ=255) ttCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc < 2:227608/149‑1 (MQ=255) tctcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTc > 2:39084/1‑81 (MQ=255) tctcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTc < 1:39084/81‑1 (MQ=255) aGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCgg > 1:112852/1‑118 (MQ=255) aGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCgg < 2:112852/118‑1 (MQ=255) cATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATAtttt < 2:1282153/126‑1 (MQ=255) cATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATAtttt > 1:1282153/1‑126 (MQ=255) tAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAatat < 2:96153/149‑1 (MQ=255) tAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATAATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAatat > 4:961368/1‑149 (MQ=255) aaTGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa > 4:357685/1‑99 (MQ=255) aaTGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa < 3:357685/99‑1 (MQ=255) gTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCt > 1:638012/1‑124 (MQ=255) gTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCt < 2:638012/124‑1 (MQ=255) aCTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattattat < 1:280100/149‑1 (MQ=255) tAAAAGCGTAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaaa > 2:1200250/1‑86 (MQ=255) tAAAAGCGTAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaaa < 1:1200250/86‑1 (MQ=255) tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc > 2:532387/1‑118 (MQ=255) tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc < 1:532387/118‑1 (MQ=255) tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattatt < 1:629187/145‑1 (MQ=255) tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattatt > 2:629187/1‑145 (MQ=255) aaaaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattattattat > 1:546639/1‑149 (MQ=255) aaGCATAACGCTTACTTTCGACATCTTTAATATATGTCATCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGa < 3:961368/149‑1 (MQ=255) aGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGAc < 2:546639/149‑1 (MQ=255) aTAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc < 4:839310/111‑1 (MQ=255) aTAACGCTTACGTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc > 3:839310/1‑111 (MQ=255) tttCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGc > 1:1094146/1‑141 (MQ=255) tttCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGc < 2:1094146/141‑1 (MQ=255) tttAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTcc < 2:413058/91‑1 (MQ=255) tttAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTcc > 1:413058/1‑91 (MQ=255) tttAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTaa > 4:11393/1‑149 (MQ=255) atatatGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAAtata < 2:1986/149‑1 (MQ=255) tatGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatc < 2:49685/59‑1 (MQ=255) tatGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatc > 1:49685/1‑59 (MQ=255) tCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCaaa > 3:366813/1‑118 (MQ=255) tCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCaaa < 4:366813/118‑1 (MQ=255) catcatGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCtat < 1:83201/63‑1 (MQ=255) catcatGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCtat > 2:83201/1‑63 (MQ=255) atcatGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCg > 1:183207/1‑112 (MQ=255) atcatGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCg < 2:183207/112‑1 (MQ=255) atGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCt > 2:1053255/1‑149 (MQ=255) gCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTgg > 3:812136/1‑149 (MQ=255) ttCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCt < 1:1133390/53‑1 (MQ=255) ttCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCt > 2:1133390/1‑53 (MQ=255) | ACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTGG > NC_000913/2102965‑2103252 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |