| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 4,013,020 | C→T | intergenic (‑204/‑33) | metR ← / → metE | DNA‑binding transcriptional dual regulator MetR/cobalamin‑independent homocysteine transmethylase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,013,020 | 0 | C | T | 96.4% | 88.1 / ‑3.7 | 28 | intergenic (‑204/‑33) | metR/metE | DNA‑binding transcriptional dual regulator MetR/cobalamin‑independent homocysteine transmethylase |
| Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (15/12); total (16/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01 | |||||||||||
AACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATCAACAAAAGCAAGCGGGTATCGACCTGCTGCCGGTGGGCGATTTTGCCTGGTACGATCAT > NC_000913/4012779‑4013259 | aaCGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGAttt > 1:97254/1‑251 (MQ=255) tgtggtgttttctgtttttttACTGCTATCATGAAAGGCCTTCCCTGCGGATTGAATATTTTGCGCTTGAAGAATACACATGAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGTATTTACaga < 2:200478/233‑3 (MQ=255) tGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACatat > 2:91876/1‑251 (MQ=255) gTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACatat < 1:91876/250‑1 (MQ=255) gcGATCATGAAAGTCCTTCACTTCGGAATGAATAATTTGCGCTGGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTGCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACaa < 1:198411/250‑1 (MQ=255) aTGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCATTCCCTAAATTCAAAATCCATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTg > 1:110897/1‑250 (MQ=255) ccTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACa > 2:178940/1‑210 (MQ=255) ccTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACa < 1:178940/210‑1 (MQ=255) atTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTCTTAAATTTATGGAGCGTTGGCTGCGTTCCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGCTGAATAAACTCACCGCCTTCCCTAAATTCAAAATCTATAGGTTTTACATATAATTAGAGGACGAAAAAATGACACTATTGAATCACACCCTCGGTTTcccctgcgttcg > 2:136585/1‑242 (MQ=255) gggggataataaaGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTcg < 2:228274/239‑1 (MQ=255) aCAGTAACCGACAATTATGGATGTGTAAACATCTGGCCGGCTACAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTAAATATAATTAGAGGAAGAAAAAATGACAATATTGACTCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGaa > 1:98201/1‑250 (MQ=255) aTTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATTAATAAACTTTCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACTCCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGaaa < 1:36170/251‑1 (MQ=255) tATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAg < 2:184247/250‑1 (MQ=255) tANGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCGTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATTAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGTATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCTGCGCGAGCTGAAAAAAGCGCAAGAAAg < 2:259789/250‑1 (MQ=255) aaCATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCgggg < 1:137793/251‑1 (MQ=255) cATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGTTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTAAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGACGCGAGCGGAAAAAAGAGCAAGTAAGTTATTGCTCGGGGa > 2:200328/1‑250 (MQ=255) aTCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCGGGGAACTCCACGCGTgaaga > 1:35190/1‑250 (MQ=255) cTTTTCAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAGTTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCATGCGTCGCGAGCTGAAAAAAGCGCACGAGAGTTATTGGGCGGGGAACTCCACGCGTGAAGAACTGCTGGCg > 2:227166/1‑251 (MQ=255) tAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGATTCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCGGGGAACTCCACGCGTGAAGAATTGCTGGCGGt > 1:235818/1‑249 (MQ=255) gCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCGGGGAACTCCACGCGTGAAGAACTGCTGGCTGGAGGGCGTGAATTg < 2:35190/250‑1 (MQ=255) gTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCGGGGAACTCCACGCGTGACGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCAc > 2:42674/1‑251 (MQ=255) gTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCAc > 1:119717/1‑251 (MQ=255) agTCACTTACTTCAGTAAGCTCCCGGGGATGAATCAACTTGCCGCCTTCCCTAAATTCAAAATCTATCGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATcaacaa > 2:141075/2‑250 (MQ=255) acttacttCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATCAACAAaagc > 1:188752/1‑250 (MQ=255) aGTACGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAACTCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGGTATTGGGCGGGGAACTCCACGCGTGAAGAACTGATGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATCAACAAAAGCAAGCGGGTAt > 2:44213/1‑251 (MQ=255) aaGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATCAACAAAAGCAAGCGGGTATCg < 1:42674/250‑1 (MQ=255) aTAAACTTGCCGCCTTCCCTAAATTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATCAACAAAAGCAAGCGGGTATCGACCTGCTGCCGGTgggg > 1:164118/1‑250 (MQ=255) aaTTCAAAATCTATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGGGGGTAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATCAACAAAAGCAAGCGGGTATCGACCTGCTGCCGGTGGGCGATTTTGCCTGGTACGATCAt < 2:265087/251‑1 (MQ=255) | AACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACATATAATTAGAGGAAGAAAAAATGACAATATTGAATCACACCCTCGGTTTCCCTCGCGTTGGCCTGCGTCGCGAGCTGAAAAAAGCGCAAGAAAGTTATTGGGCGGGGAACTCCACGCGTGAAGAACTGCTGGCGGTAGGGCGTGAATTGCGTGCTCGTCACTGGGATCAACAAAAGCAAGCGGGTATCGACCTGCTGCCGGTGGGCGATTTTGCCTGGTACGATCAT > NC_000913/4012779‑4013259 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |