Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 718,461 | G→A | 9.1% | intergenic (+31/+46) | thiO → / ← PP_0613 | FAD‑dependent glycine/D‑amino acid oxidase/amidase family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 718,461 | 0 | G | A | 9.1% | 115.0 / 2.6 | 49 | intergenic (+31/+46) | thiO/PP_0613 | FAD‑dependent glycine/D‑amino acid oxidase/amidase family protein |
Reads supporting (aligned to +/- strand): ref base G (16/26); new base A (2/2); total (18/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.39e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-02 |
ATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCAC > NC_002947/718315‑718605 | aTCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTAcgc > 2:296376/1‑148 (MQ=255) tCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTAcgc > 7:65360/1‑147 (MQ=255) cGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAg > 1:82666/1‑149 (MQ=255) aaTGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCgg < 3:135674/149‑1 (MQ=255) ggCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTTAACCCGCGCCGAAGCGGGGCGGGtt > 7:142363/1‑149 (MQ=255) gCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTa < 1:302548/149‑1 (MQ=255) cTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑c < 1:296376/149‑1 (MQ=255) gTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑ccc < 5:59977/148‑1 (MQ=255) cTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGttt < 4:106776/131‑1 (MQ=255) cTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGttt > 3:106776/1‑131 (MQ=255) tttACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCGGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGa < 7:315460/149‑1 (MQ=255) tttACCGATTTGCTGACAGGCGCCGACCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGCACCTCTCCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTAA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGa < 1:123977/149‑1 (MQ=39) cGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCTTATGCGCCAGAAGGGCGTTTGGGCTGAGTACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGTTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa < 7:28754/149‑1 (MQ=255) cGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCTTATGCGCCAGAAGGGCGTTTGGGCTGAGCACCTCTTCGCGTTTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa < 1:95891/149‑1 (MQ=255) cGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa < 8:65360/149‑1 (MQ=255) cGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa < 5:186243/149‑1 (MQ=255) cGATTTGCTCACAGGGGCCCAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGCACCTCTTCGCGGGTAACCCCGCTCCAACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa < 7:195812/149‑1 (MQ=255) gATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGGTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGAtttacg > 6:183808/1‑146 (MQ=255) aGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCGCCTACGC‑GAGCGGGTTTT‑CCCGCGAAAGGGCCGCAGTAGCCAATAAAGATTTAATCCAGTAATCgg > 8:75121/1‑149 (MQ=255) ggggCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCACCTACGC‑GAGCGGGGTTA‑CCCCCGAAAAGGCCGCACAACCCAAAAAAGATTGAAGCCAGTAAACGGt > 3:75685/1‑149 (MQ=255) ggCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCGCTACAACGC‑GAGCGGGTTTA‑ACCGCGAAAGGGCCGCAGAAGCCAAAACAGATTTAAGCCAGTAAAAGGTGc > 6:144671/1‑149 (MQ=255) ggccgaTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAG‑GAGCGGGTTTA‑c < 2:89735/93‑1 (MQ=25) ggccgaTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCGCGTAG‑GAGCGGGTTTA‑c > 1:89735/2‑94 (MQ=25) cgaTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGt < 5:183808/149‑1 (MQ=255) aTTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTCCGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGc < 2:260830/149‑1 (MQ=255) gACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCTCCAACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATAAAGATTTGAGCCAGTAATAGGTGGTTGCAGGTGCCGGcgc > 6:119289/1‑147 (MQ=255) gACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCGCCGAAGC‑GAGCGGGGTTT‑ACCGCGAAAGGGGCGCCGAAGGCCATAAAGAGTTCAGCCAGTAATAGGTGGTTGCAGGTggcggc > 4:79136/1‑147 (MQ=255) cgggggcCTTATGCGCCAGAAGGGCTTTTGGGCGGTGGACCTCTTCGCGGGTAAACCCGCTCCTCCGC‑TAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGa < 3:228175/144‑1 (MQ=38) gggCGCCTTATCCCCCAGACGCGGGTTTGCGCGCAGCACGGCTTCGCGGTTAAACCCGCCCCACCGC‑GAGCGGTTTAA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGAt < 3:29266/148‑1 (MQ=255) gcCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAAtg < 7:247735/149‑1 (MQ=255) tATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAAACCCCTCCTACGC‑GAGCGGGTTTA‑ACCGCGAAAGGGCCGCAGAAGCCAAAACAGATTTAAGCCAGTAATAGGTGCTTGCAGGTGCCGGCTCGGCTAAtgctgc > 2:228366/1‑149 (MQ=255) atactggggGTTTGGGCTGCGAACGCCTGCG‑GGGTAAAACCGCCACTACCC‑CCGGGGGTTTACCCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATaca < 1:367396/144‑1 (MQ=255) gggCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTg < 2:82666/149‑1 (MQ=255) ggCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCgc < 1:228366/149‑1 (MQ=255) cTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCgcag < 5:144671/149‑1 (MQ=255) ccgctgcgcgGGTAAACCCGCCCCACCCC‑CAGCCGGATTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGc < 8:142363/143‑1 (MQ=255) ccTCTTCGCGGGTAACCCCGCTCCAACGC‑TCGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGc < 3:79136/149‑1 (MQ=255) tctTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGTTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCa < 3:125841/74‑1 (MQ=255) gcgGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAg > 3:211197/1‑149 (MQ=255) gagtgaaaaCACGCTCATACGC‑CAGCGGTTTTA‑CCCGCGAAAGTGCCGCAGAATCCAATACATATTTCAGCCATTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAg < 6:117778/143‑1 (MQ=255) ggTAAACCCGCTCGCGTAG‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCtcatca > 4:208397/1‑145 (MQ=39) ggTAAACCCGCTCGCGTAG‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCtcatca < 3:208397/145‑1 (MQ=39) aaaCCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAAtg < 4:141095/102‑1 (MQ=255) aaaCCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAAtg > 3:141095/1‑102 (MQ=255) aCCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGc < 6:68743/117‑1 (MQ=255) aCCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGc > 5:68743/1‑117 (MQ=255) ccGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGccc > 6:227171/1‑149 (MQ=255) ccGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGccc > 2:272339/1‑149 (MQ=255) ccGCTCCTACCC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGccc < 7:75121/149‑1 (MQ=255) cTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGc < 5:119289/149‑1 (MQ=255) cTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGccac > 8:80549/1‑149 (MQ=255) | ATCGCAATGGGCTGGTGCTGGCGCCGGCTTCGTGCCAGCTGTTTACCGATTTGCTGACAGGGGCCGAGCCGATTATCGACCCGGCGCCGTATGCGCCAGAAGGGCGTTTGGGCTGAGGACCTCTTCGCGGGTAAACCCGCTCCTACGC‑GAGCGGGTTTA‑CCCGCGAAAGGGCCGCAGAAGCCAATACAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCAC > NC_002947/718315‑718605 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |