Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,285,261 | Δ1 bp | 10.3% | coding (123/711 nt) | narJ → | molybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,285,261 | 0 | T | . | 10.3% | 197.0 / 17.4 | 58 | coding (123/711 nt) | narJ | molybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1 |
Reads supporting (aligned to +/- strand): ref base T (27/25); new base . (3/3); total (30/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01 |
AGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGC > NC_000913/1285119‑1285407 | aGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑t > 2:493365/1‑148 (MQ=255) gCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑Tg < 1:493365/148‑1 (MQ=255) gAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCAtttt < 3:1036419/148‑1 (MQ=255) tGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCAtt > 2:425247/1‑103 (MQ=255) tGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCAtt < 1:425247/103‑1 (MQ=255) atatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGa > 3:761129/1‑111 (MQ=255) atatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGa < 4:761129/111‑1 (MQ=255) tatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAt < 4:1061/98‑1 (MQ=255) tatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAt > 3:1061/1‑98 (MQ=255) aTGCTGCCATAAGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa > 3:1015131/1‑148 (MQ=255) gcCTTATGGCAGCATCAAGAAGAGATGTTTGAGGCGATTGCCGCGTCGAAATATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATgg > 4:717623/1‑112 (MQ=255) gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATgg < 3:717623/112‑1 (MQ=255) gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAtc > 3:380005/1‑86 (MQ=255) gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAtc < 4:380005/87‑2 (MQ=255) cTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCt > 3:108979/1‑117 (MQ=255) cTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCt < 4:108979/117‑1 (MQ=255) ttATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAt > 1:917031/1‑83 (MQ=255) ttATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAt < 2:917031/83‑1 (MQ=255) ggCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTaa < 3:111793/96‑1 (MQ=255) ggCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTaa > 4:111793/1‑96 (MQ=255) cagcaTCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGta < 4:1015131/148‑2 (MQ=255) acaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGcgcg > 1:1037423/1‑80 (MQ=255) acaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGcgcg < 2:1037423/80‑1 (MQ=255) acaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTAc < 4:448459/121‑1 (MQ=255) acaaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTAc > 3:448459/1‑121 (MQ=255) aaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCTTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcgc > 1:111013/1‑114 (MQ=255) aaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGAt < 4:934324/95‑1 (MQ=255) aaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcgc < 2:111013/114‑1 (MQ=255) aaGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCAAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGAt > 3:934324/1‑95 (MQ=255) agaTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCa < 4:1065852/106‑1 (MQ=255) agaTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCa > 3:1065852/1‑106 (MQ=255) agaTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCa < 2:311176/106‑1 (MQ=255) agaTGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCa > 1:311176/1‑106 (MQ=255) ggCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCcgcg < 4:849117/131‑1 (MQ=255) ggCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCcgcg > 3:849117/1‑131 (MQ=255) cGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGAc < 4:476322/119‑1 (MQ=255) cGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGAc > 3:476322/1‑119 (MQ=255) aTTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa > 4:333696/1‑107 (MQ=255) aTTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa < 3:333696/107‑1 (MQ=255) ttGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCTCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGtt < 3:954925/112‑1 (MQ=255) ttGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGtt > 4:954925/1‑112 (MQ=255) aaaTCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtg > 1:769884/1‑140 (MQ=255) aaaTCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtg < 2:769884/140‑1 (MQ=255) aaTCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGttt > 3:504594/1‑132 (MQ=255) aaTCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGttt < 4:504594/132‑1 (MQ=255) cTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa < 2:237957/89‑1 (MQ=255) cTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa > 1:237957/1‑89 (MQ=255) gCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACa > 4:109526/1‑132 (MQ=255) gCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACa < 3:109526/132‑1 (MQ=255) aaaaGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGaa < 2:938560/142‑1 (MQ=255) aaaaGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGaa > 1:938560/1‑142 (MQ=255) aaaGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgcg > 3:834586/1‑148 (MQ=255) aaGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCa > 2:858651/1‑67 (MQ=255) aaGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCa < 1:858651/67‑1 (MQ=255) agGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcg > 2:865786/1‑148 (MQ=255) aTGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGt > 2:632083/1‑81 (MQ=255) aTGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGt < 1:632083/81‑1 (MQ=255) aTGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc > 2:866763/1‑148 (MQ=255) | AGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGC > NC_000913/1285119‑1285407 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |