Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,534,904 | T→G | 8.5% | intergenic (+35/+61) | nhoA → / ← yddE | N‑hydroxyarylamine O‑acetyltransferase/PhzC‑PhzF family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,534,904 | 0 | T | G | 8.5% | 199.3 / 10.8 | 59 | intergenic (+35/+61) | nhoA/yddE | N‑hydroxyarylamine O‑acetyltransferase/PhzC‑PhzF family protein |
Reads supporting (aligned to +/- strand): ref base T (27/27); new base G (3/2); total (30/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01 |
AGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTC > NC_000913/1534764‑1535052 | aGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc > 1:172029/1‑149 (MQ=255) gAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCa < 3:1044196/149‑1 (MQ=255) cAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATAc < 2:172029/149‑1 (MQ=255) gatgatGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATcc > 3:911749/1‑124 (MQ=255) gatgatGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGACGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATcc < 4:911749/124‑1 (MQ=255) tgatgCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCt < 3:357036/149‑1 (MQ=255) atgCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTa > 3:38333/1‑117 (MQ=255) atgCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTa < 4:38333/117‑1 (MQ=255) gAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCAGACATTAGCCCGGTTGTAGACCTGACTgg > 3:385450/1‑149 (MQ=255) aaaCATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGtt < 4:1132798/134‑1 (MQ=255) aaaCATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGtt > 3:1132798/1‑134 (MQ=255) aaaCATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATAAGGCGTTTACGCCGCATCCgg > 3:445570/1‑117 (MQ=38) cTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACa > 3:372052/1‑115 (MQ=255) cTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACa < 4:372052/115‑1 (MQ=255) tACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCa > 3:1248786/1‑149 (MQ=255) cGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATcc < 3:1040657/101‑1 (MQ=255) cGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATcc > 4:1040657/1‑101 (MQ=255) gacgAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACTCGGAGGCAGGAAAATAATGTATGTCAGCTTGCCGAATGCGGCGTAAACGCCTTACCTGGCATACATTAGCCTGGTTGTAGACCTGACAGGCGAAGCGCATAAGATAACg < 4:385450/146‑1 (MQ=255) ttAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCa < 2:630486/95‑1 (MQ=255) ttAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCa > 1:630486/1‑95 (MQ=255) ttAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCa < 4:1248786/149‑1 (MQ=255) ttAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCa > 1:355242/1‑149 (MQ=255) gTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCaga > 1:634364/1‑142 (MQ=255) gTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCaga < 2:634364/142‑1 (MQ=255) gcgTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTc < 2:355242/149‑1 (MQ=255) aCCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCTGGCATACATTa < 1:1087735/71‑1 (MQ=255) aCCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTa > 2:1087735/1‑71 (MQ=255) aaaaTAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCt < 1:1138723/42‑1 (MQ=255) aaaaTAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCt > 2:1138723/1‑42 (MQ=255) aataatTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGc < 1:480744/115‑1 (MQ=255) aataatTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGc > 2:480744/1‑115 (MQ=255) tttATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaaa < 4:499059/130‑1 (MQ=255) tttATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaaa > 3:499059/1‑130 (MQ=255) tttATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAAt > 2:1348249/1‑149 (MQ=255) tttATGTCAGGCTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGGCCATTCGGCATGAATCAAAATCACCGCTGTGCCTGAAAt > 3:961408/1‑149 (MQ=255) tGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTg > 3:952820/1‑149 (MQ=255) gTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCaga > 1:1291911/1‑95 (MQ=255) gTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCaga < 2:1291911/95‑1 (MQ=255) gTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGa < 4:952820/149‑1 (MQ=255) tCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCa > 3:461352/1‑92 (MQ=255) tCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCa < 4:461352/92‑1 (MQ=255) tCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGAc < 1:1348249/149‑1 (MQ=255) gCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTcc < 4:961408/149‑1 (MQ=255) cGGATAAGGCGTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaaa > 1:419047/1‑115 (MQ=38) cGGATAAGGCGTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaaa < 2:419047/115‑1 (MQ=38) gataaggcgTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCt < 3:116578/122‑1 (MQ=255) gataaggcgTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCt > 4:116578/6‑127 (MQ=255) gATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaa < 3:1427413/112‑1 (MQ=255) gATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaa > 4:1427413/1‑112 (MQ=255) gcgTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTAt > 3:1392106/1‑149 (MQ=255) cgTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATc > 1:113160/1‑149 (MQ=255) tAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCAc > 3:522723/1‑149 (MQ=255) aaaCGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACg > 4:1314782/1‑149 (MQ=255) aaCGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGa < 4:1392106/149‑1 (MQ=255) aCGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACTGTTTCAGAGTTCaa < 3:802485/77‑1 (MQ=255) aCGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCaa > 4:802485/1‑77 (MQ=255) aCGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCTCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCggtt > 3:991604/1‑137 (MQ=255) aCGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCggtt < 4:991604/137‑1 (MQ=255) aCGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGaa < 4:522723/149‑1 (MQ=255) gCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAAtt < 4:1306711/149‑1 (MQ=255) gCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAAtt > 3:1306711/1‑149 (MQ=255) ttATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTc > 3:1192803/1‑149 (MQ=255) ttATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTc > 2:860210/1‑149 (MQ=255) | AGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTC > NC_000913/1534764‑1535052 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |