Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,903,077 | T→A | 11.2% | intergenic (+42/+73) | PP_3445 → / ← ilvA‑I | hypothetical protein/threonine deaminase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,903,077 | 0 | T | A | 11.2% | 116.5 / 5.8 | 45 | intergenic (+42/+73) | PP_3445/ilvA‑I | hypothetical protein/threonine deaminase |
Reads supporting (aligned to +/- strand): ref base T (16/23); new base A (3/2); total (20/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.38e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01 |
CCCCGAGCCGGATACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCATCTCGTC > NC_002947/3902936‑3903215 | ccccGAGCCGGATACCGAGAACTTCAAGCAGCTCCGCCGACTTCAAACGCCTGTAAACGTTCCGCAAGTGCGTTGCCGCCTTTTTCATCCGCGAGCGTTGAAATTGGCTCCCCCAGTTCCCTCCAGGAGATCGCCATCCCCTGcgggcg > 1:216055/1‑149 (MQ=255) cccGAGCCGGATACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTcgggcggg < 1:80829/149‑1 (MQ=255) aGCCGGATACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGggcggggcgg < 2:133632/149‑1 (MQ=255) gATACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACg > 5:310097/1‑149 (MQ=255) aTACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCCGCCCGAAgggcagggcga < 5:293427/143‑2 (MQ=255) aTACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCCGCCCGAAgggcagggcga > 6:293427/1‑142 (MQ=255) tACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACg > 7:182632/1‑147 (MQ=255) aCCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAg < 7:312639/149‑1 (MQ=255) cctagcaCTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGACACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGc < 4:212232/143‑1 (MQ=255) ccGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGc > 2:224967/1‑149 (MQ=255) ccGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGc > 5:288284/1‑149 (MQ=255) cGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCg < 1:224967/149‑1 (MQ=255) cGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCg < 6:288284/149‑1 (MQ=255) gaACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAac < 8:182632/149‑1 (MQ=255) cAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcaggc > 8:191844/1‑149 (MQ=255) cAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcagg < 1:81457/148‑1 (MQ=255) cAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcagg < 8:337084/148‑1 (MQ=255) cGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcaggcggcaggcgg > 5:63674/1‑149 (MQ=255) cGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcaggcggcaggcgg < 5:348467/149‑1 (MQ=255) gACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAacac > 2:233179/1‑132 (MQ=255) gACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAacac < 1:233179/132‑1 (MQ=255) aCTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGggcggggcgg > 2:75441/1‑116 (MQ=255) aCTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGggcggggcgg < 1:75441/116‑1 (MQ=255) cTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGAc < 6:94295/149‑1 (MQ=255) tCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTa > 5:182918/1‑144 (MQ=255) tCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTa < 6:182918/144‑1 (MQ=255) aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc > 8:181365/1‑149 (MQ=255) aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc > 4:14914/1‑149 (MQ=255) aaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGt < 6:187330/149‑1 (MQ=255) aaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGt > 6:55082/1‑149 (MQ=255) ttCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGcac < 6:467/149‑1 (MQ=255) aGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGcgcg < 3:14914/149‑1 (MQ=255) gTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGc < 7:181365/149‑1 (MQ=255) ccGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAg > 1:134000/1‑149 (MQ=255) tttCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAg < 6:63674/148‑1 (MQ=255) cgAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGcc < 2:216055/149‑1 (MQ=255) tCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCGTCGGGCGGGGCGGTCACGCAGCGAACACCGGCGGCAGGCGGCTAGCCAGGCGTGCATTGGCCGCACAGGCGCGTCCGCCCAGGAACCGCGTGTACGCCGGGTTGGCAGTCGCGTc > 5:300293/1‑149 (MQ=255) gcgcAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTcc < 6:310097/148‑1 (MQ=255) ggAGATCGCCCTGCCCTTCGGGCGTGGCGGTCACGCAGCGAACACAGTCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCa < 4:24143/149‑1 (MQ=255) gaccgccccgcccGAAGGGCAGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCAtct > 8:273687/4‑149 (MQ=255) accgccccgcccGAAGGGCAGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTAc < 6:119353/133‑1 (MQ=255) accgccccgcccGAAGGGCAGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTAc > 5:119353/3‑135 (MQ=255) ccctgccctTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCg > 2:334730/1‑108 (MQ=255) ccctgccctTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCg < 1:334730/108‑1 (MQ=255) ccctgcccgTAGGGCCGGGCGGTAACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCCTGCATTTGCCGCCCAGGCCCGTCCGCCCAGGAACAGCTTGTACCCCGGGTTGTGAGTCTCGTACCAGGACCCGGACCCCaacaagta > 4:2289/1‑141 (MQ=255) | CCCCGAGCCGGATACCGAGAACTTCAAGC‑GCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCATCTCGTC > NC_002947/3902936‑3903215 |
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: ‑ |