Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,987,442 | Δ1 bp | 15.1% | intergenic (+14/+65) | ftnB → / ← yecJ | ferritin B, putative ferrous iron reservoir/DUF2766 family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,987,442 | 0 | T | . | 15.1% | 204.9 / 33.5 | 73 | intergenic (+14/+65) | ftnB/yecJ | ferritin B, putative ferrous iron reservoir/DUF2766 family protein |
Reads supporting (aligned to +/- strand): ref base T (34/28); new base . (5/6); total (39/34) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.45e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01 |
AAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTT > NC_000913/1987343‑1987541 | aaGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCt < 1:2131377/100‑1 (MQ=255) aGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTa > 1:926520/1‑100 (MQ=255) gTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa > 1:1094646/1‑100 (MQ=255) gTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa < 1:286903/100‑1 (MQ=255) ctcAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGc > 1:1613023/3‑100 (MQ=255) gTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCAtt > 1:1783470/1‑99 (MQ=255) gTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGC‑AATGCAttag < 1:1818316/100‑3 (MQ=255) aaaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCg < 2:360797/100‑1 (MQ=255) aaaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTAACACCAGCTGCATTGATCATCATATGCGCTAATGCATT‑tcgc > 2:229744/1‑98 (MQ=255) aaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGC‑AATGCATTAGCGCCGa < 2:1016284/100‑1 (MQ=39) aaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCGGCATTGATCATCATCGGCGC‑AATGCATTAGCGCCGa < 2:1700142/100‑1 (MQ=39) aCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGAt < 1:733418/100‑1 (MQ=255) aCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGC‑AATGCATTAGCGCCGAt < 1:1555669/100‑1 (MQ=39) aCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCAATGATCATCATCGGCGC‑AATGCATTAGCGCCGAt < 2:2262770/100‑1 (MQ=38) cTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATg > 2:1961514/1‑100 (MQ=255) cTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATg > 2:292065/1‑100 (MQ=255) cGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAgg > 1:386730/1‑100 (MQ=255) cGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAgg > 2:1318362/1‑100 (MQ=255) cGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAgg < 1:813652/100‑1 (MQ=255) tgtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTga < 1:1421966/100‑1 (MQ=255) tgtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGATTTga > 1:900650/1‑100 (MQ=255) tgtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGATTTga > 1:61116/1‑100 (MQ=255) tgtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑ACGCCGATGAAAGTTTTga > 1:809454/1‑100 (MQ=255) gtgCTCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTgag < 1:1016784/100‑1 (MQ=255) gtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTTGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTgag < 1:1688036/100‑1 (MQ=255) gtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTgag < 1:2169547/100‑1 (MQ=255) gCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGaa > 1:312971/1‑100 (MQ=255) gCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGaa > 1:1098460/1‑100 (MQ=255) ccTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAAc < 1:1696754/100‑1 (MQ=255) gtgCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATAATCGGCGCTAATGCATT‑GCACCGATGAAGGTATTGAGAAACCg > 1:904566/1‑99 (MQ=255) gaccgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCt < 2:1947851/100‑1 (MQ=255) ccgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCa > 1:344572/1‑100 (MQ=255) gaccAACATGGTCTGTATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGCGCAGCCGCTGCctcctc > 2:1639910/1‑100 (MQ=255) ccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATc > 2:1346712/1‑98 (MQ=255) aaCATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCAGCATCGGCGCTAATGCATT‑GAGCCGATGAAGGTTTTGAGAAACCGCTGCCGCATCTGt > 2:1303662/1‑99 (MQ=255) cATGTTCTGAATGTCGTGTCACACCCGCTGCATTGCTCATCATCGGCGCTCATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTg < 2:705563/100‑1 (MQ=255) cATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTg < 2:31905/100‑1 (MQ=255) cATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTg > 2:379651/1‑100 (MQ=255) tGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCTTCaaacagatgag < 2:386727/100‑11 (MQ=255) tGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGaa > 1:84135/1‑100 (MQ=255) tGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAagcagc < 1:110615/100‑1 (MQ=255) aTGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGt > 1:1384379/1‑100 (MQ=255) gTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttt > 1:2136055/1‑100 (MQ=255) tCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttt < 1:306756/100‑1 (MQ=255) gtgtCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttttt < 2:867889/100‑1 (MQ=255) tgtCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttttttt < 1:993519/100‑1 (MQ=255) cacaCCAGGTGCATTGATAATCATCGGAGCTAATGCATT‑GCGTCGATGAAGGTATTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGGTGTTTTAAt < 2:34622/100‑1 (MQ=255) cacaCCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt > 1:2398213/1‑100 (MQ=255) cacaCCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt > 1:2398061/1‑100 (MQ=255) acacCAGCTGCATTGATCATCATCGGCGCTAATGCGTT‑GCGTCGATGAAGGTTTTGTGAATTCGTTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt < 2:499339/99‑1 (MQ=255) acCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATggg < 2:413678/100‑1 (MQ=255) acCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATgg < 1:610770/99‑1 (MQ=255) ccAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGTGa > 1:1945966/1‑100 (MQ=255) ccAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGa > 1:1685702/1‑100 (MQ=255) cAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGAt > 1:1391980/1‑100 (MQ=255) aGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGAtt < 1:220691/100‑1 (MQ=255) gCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTc > 2:2316062/1‑100 (MQ=255) gCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAcc > 2:553649/1‑100 (MQ=255) ctcaaaaccttcatcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAcc > 1:1724361/11‑100 (MQ=25) cATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAccc < 2:210714/100‑1 (MQ=255) gATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtgtg > 1:1251392/1‑100 (MQ=255) atcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTTTTTGAAGCAGCGGTTTTTTTAATGGGATACACCCTGTgggct < 1:1622513/100‑3 (MQ=255) atcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGt < 2:1603446/100‑1 (MQ=255) catcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTaa < 2:639244/100‑1 (MQ=255) catcGGCGC‑AATGCATTAGCGCCGCTGAAGGTTCTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTa > 1:2322685/1‑99 (MQ=38) atcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTaaa > 1:2260168/1‑100 (MQ=39) atcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTaaa > 1:790350/1‑100 (MQ=39) atcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTaaa > 2:493136/1‑100 (MQ=39) tcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAAc < 2:1724340/100‑1 (MQ=39) ggCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGAGTAAAc > 1:18527/1‑98 (MQ=255) cgcTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGa > 1:2038869/1‑100 (MQ=255) gcTAATTCATT‑GTGGCAATGAAGATTTTGAGAAACCGCTGCCTCATCTGTTAGAAGAAGCGTGTTATTTAATGGGTGTCAGCCGCTGGGGTAAAc < 2:1597504/95‑1 (MQ=255) gcTAATGCATT‑GCGCCGATGAAGGTTTTGTGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAg < 2:1444455/100‑1 (MQ=255) cTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAg > 1:412531/1‑99 (MQ=255) tAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTTATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGt < 2:409071/99‑1 (MQ=255) tAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGtt > 1:1484912/1‑100 (MQ=255) tAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGtt < 1:118671/100‑1 (MQ=255) | AAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTT > NC_000913/1987343‑1987541 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |