Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 127,123 | G→C | 12.6% | intergenic (‑33/+37) | thrB ← / ← PP_0122 | homoserine kinase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 127,123 | 0 | G | C | 12.6% | 118.6 / 10.5 | 48 | intergenic (‑33/+37) | thrB/PP_0122 | homoserine kinase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (20/21); new base C (3/3); total (24/24) | |||||||||||
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.96e-01 |
AAATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCT > NC_002947/126975‑127260 | aaaTTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTccgccg > 7:210562/1‑149 (MQ=255) gCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAg > 2:80131/1‑149 (MQ=255) ttCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCCGGCGGAAGTGcc < 2:100732/149‑1 (MQ=255) ttCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGAGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGcc > 7:239810/1‑149 (MQ=255) ctcGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCgg > 3:241500/1‑149 (MQ=255) tcGGTACCGGCGGCAATGCCCTTGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACTGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCgg < 8:210562/148‑1 (MQ=255) ggTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGtt < 4:241500/149‑1 (MQ=255) gctACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGATGACTGACAAGATGAAAAATTGCCCATACCGCCACTTCCCGGCGGAAGTGcccgtat < 1:207923/147‑6 (MQ=255) ggcggcAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGAtgtg < 1:131563/149‑1 (MQ=255) gCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCTGCCGGGAAGTGCCGGGTTGATGTGAACGGTgcg > 5:155974/1‑149 (MQ=255) ggAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGAt > 6:70854/1‑132 (MQ=255) ggAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGAt < 5:70854/132‑1 (MQ=255) gAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGt > 3:24633/1‑146 (MQ=255) gAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCa < 4:24633/148‑1 (MQ=255) gAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAg > 5:27415/1‑149 (MQ=255) aaGTCGAGCAGACGGCCCAGCTCGTAGGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGc < 7:199463/149‑1 (MQ=255) aGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCGGCt < 1:80131/149‑1 (MQ=255) tCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTa < 3:118102/149‑1 (MQ=255) gAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTAc < 6:176933/148‑1 (MQ=255) gACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTc > 8:244961/1‑148 (MQ=255) aCGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCg > 7:117380/1‑148 (MQ=255) gCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTc < 7:204576/149‑1 (MQ=255) gCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTc < 5:253112/149‑1 (MQ=255) cGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAg > 7:137089/1‑149 (MQ=255) gTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCa < 8:233181/147‑1 (MQ=255) gTACGGCGCCAGAAAGGTTTCAAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCAATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTTCGCTTCAGCTTACAAATCGAAGATCTTCCAg > 8:190155/1‑148 (MQ=255) cGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCAGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACg < 6:155974/149‑1 (MQ=255) gcgcCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGa < 8:137089/149‑1 (MQ=255) aaGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCa < 6:27415/148‑1 (MQ=255) cAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCt > 5:272336/1‑149 (MQ=255) gCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGAGGGAATCAGCATGTCAGGCTGGTCGGAt > 7:91107/1‑149 (MQ=255) gAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCAt > 5:143344/1‑149 (MQ=255) cTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTg < 8:117380/149‑1 (MQ=255) atgaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAg < 7:244961/148‑1 (MQ=255) tgaAAAATTGCCCATACGGGCACTTACGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAgc < 7:190155/148‑1 (MQ=255) gaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAgcgc > 8:234493/1‑149 (MQ=255) cccATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAAtt > 4:164254/1‑120 (MQ=255) cccATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAAtt < 3:164254/120‑1 (MQ=255) ccATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTACGGCTGGTCGGATCGAATGAAATTGGCATCACTGCAAt > 8:190068/1‑133 (MQ=255) ccATACGGGCACTTCCGCCAGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGaa > 3:258731/1‑148 (MQ=255) aTACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAAt > 5:62071/1‑111 (MQ=255) aTACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAAt < 6:62071/111‑1 (MQ=255) aTACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGaaaa > 1:76296/1‑148 (MQ=255) ccggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAg > 7:57542/3‑149 (MQ=255) cGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATaa < 8:239810/148‑1 (MQ=255) cggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGaaa < 7:251813/142‑1 (MQ=255) cggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGaaa > 8:251813/2‑143 (MQ=255) cggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAgg > 6:53780/2‑149 (MQ=255) gCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCt > 5:36830/1‑149 (MQ=255) | AAATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCT > NC_002947/126975‑127260 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |