| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 1,571,721 | G→A | 8.0% | H109Y (CAT→TAT) | gadB ← | glutamate decarboxylase B, PLP‑dependent |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,571,721 | 0 | G | A | 8.0% | 142.0 / 2.7 | 50 | H109Y (CAT→TAT) | gadB | glutamate decarboxylase B, PLP‑dependent |
| Reads supporting (aligned to +/- strand): ref base G (26/20); new base A (2/2); total (28/22) | |||||||||||
| 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.00e+00 | |||||||||||
CAGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGC > NC_000913/1571657‑1571855 | cAGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGt < 3:132609/148‑1 (MQ=31) aGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGtt > 3:86585/1‑148 (MQ=32) aGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGtt > 3:3964/1‑148 (MQ=32) aGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGtt < 4:3964/148‑1 (MQ=32) gCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAAt < 8:396663/149‑1 (MQ=34) gCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAAt > 2:457794/1‑149 (MQ=34) ccTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATg < 7:106546/149‑1 (MQ=35) cTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATg > 1:341831/1‑148 (MQ=35) gaagaaCCAATGGTGTTTGTGCCAAAGGCATGACCATTTTTCGGCGCAGTCGCATGCCACAGATCGGCAACCATATTTACGGAACGCAGGTCGATGGGTGCGGATTGCGGATATTCTTCTTGGTCGCTCCAGTTTTTGGTAATGGATa > 6:308632/1‑148 (MQ=11) aagaaCCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATaa > 3:232086/1‑148 (MQ=39) aagaaCCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATaa < 4:232086/148‑1 (MQ=39) aagaaCCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATACCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATaaa > 8:437598/1‑149 (MQ=255) aagaaCCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATACCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATaaa > 8:437603/1‑149 (MQ=255) ccAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCAt < 6:115828/149‑1 (MQ=39) ccAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCAt > 5:436671/1‑149 (MQ=39) ccAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCAt < 2:341831/149‑1 (MQ=39) ccAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCAt < 4:86585/149‑1 (MQ=39) cAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCAt > 7:114680/1‑148 (MQ=39) aaTGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGGCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCa < 4:418829/149‑1 (MQ=255) aaTGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATc > 2:442681/1‑148 (MQ=39) ttGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTtgtg > 1:256727/1‑149 (MQ=255) ggTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTgg < 1:442681/148‑1 (MQ=255) gTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGa > 4:365523/1‑148 (MQ=255) tGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACa > 4:371230/1‑149 (MQ=255) gCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTg > 4:198643/1‑118 (MQ=11) gCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTg < 3:198643/118‑1 (MQ=11) gCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACAt > 2:34632/1‑149 (MQ=255) cAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTGCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACAttt > 6:155216/1‑149 (MQ=255) aaCGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACAtttt > 1:434179/1‑149 (MQ=255) aaCGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACAtttt > 5:84508/1‑149 (MQ=255) ggCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAAt > 3:97155/1‑117 (MQ=35) ggCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAAt < 4:97155/117‑1 (MQ=35) cTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTcgtcgt > 1:13383/1‑149 (MQ=255) cTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTcgtcgt < 6:431368/149‑1 (MQ=255) cTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTcgtcgt > 4:390558/1‑149 (MQ=255) tGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATa > 3:266050/1‑118 (MQ=40) tGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATa < 4:266050/118‑1 (MQ=40) aCCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCgtcgtc < 2:65516/147‑1 (MQ=255) aCCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCgtcgtc < 2:13383/147‑1 (MQ=255) aCCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTcc > 7:304191/1‑148 (MQ=255) aCCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTcc > 7:49502/1‑148 (MQ=255) cATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAg < 8:114680/149‑1 (MQ=255) aTTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAgg > 7:208793/1‑149 (MQ=255) tttttCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGt > 2:63098/1‑149 (MQ=255) ttCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCt < 3:371230/148‑1 (MQ=255) tCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTg < 3:390558/148‑1 (MQ=255) tCGGCGCAGGCGCATACCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTg < 7:437598/148‑1 (MQ=255) tCGGCGCAGGCGCATACCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTg < 7:437603/148‑1 (MQ=255) ggcgcaggcgcaTGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGc < 2:469826/148‑1 (MQ=255) ggcccaggcgcaTGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCATCTTTGTCGATGCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCATCGTCCCAGGTCTgg > 4:22112/1‑147 (MQ=255) | CAGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGC > NC_000913/1571657‑1571855 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |