Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,352,751 | C→T | 41.9% | R293H (CGT→CAT) | sapD ← | putrescine ABC exporter ATP binding protein SapD |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,352,751 | 0 | C | T | 41.9% | 9.6 / 24.7 | 30 | R293H (CGT→CAT) | sapD | putrescine ABC exporter ATP binding protein SapD |
Reads supporting (aligned to +/- strand): ref base C (8/9); new base T (6/6); total (15/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01 |
AGCAGCGTTTCGATCATCTCACTCTTTCTCCATGTTCAGCGGGAAATGACAGGCATAGAGATGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGACGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTC > NC_000913/1352618‑1352888 | agcagcGTTTCGATCATCTCACTCTTTCTCCATGTTCAGCGGGAAATGACAGGCATAGAGATGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGATGACAACCaa < 1:233477/143‑1 (MQ=255) gcagcGTTTCGATCATCTCACTCTTTCTCCATGTTCAGCGGGAAATGACAGGCATAGAGATGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCACTGTGGCCCCAGATGACAACCCAt > 2:70219/1‑143 (MQ=255) gcagcGTTTCGATCATCTCACTCTTTCTCCATGTTCAGCGGGAAATGACAGGCATAGAGATGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGATGACAACCAAt > 2:601236/1‑143 (MQ=255) tCGATCATCTCACTCTTTCTCCATGTTCAGCGGGAAATGACAGGCATAGAGATGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGACGACAACCAATCGGTAACt < 2:187309/143‑1 (MQ=255) gATCATCTCACTCTTTCTCCATGTTCAGCGGGAAACGCCAGGCATAGAGATGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGATGACAACCAATCGGTAACTGt < 1:140588/143‑1 (MQ=255) ttCTCCATGTTCAGCGGGAAATGACAGGCATAGAGATGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGACGACAACCAATCGGTAACTGGTCCCGCAGTGGGAt > 2:623447/1‑143 (MQ=255) cATGTTCAGCGGGAAATGACAGGCATAGAGATGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGACGACAACCCATCGGTAGCTGTTTCAGCAGTGGGATAgcgc > 2:3560/1‑143 (MQ=255) cATGTGCAGGGGGAAATGACAGGCATAGAGATGATTTTTCGCCCCCGTCAAGCGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGACGACAACCAATCTGTAACTGTTCCAGCAGTGGGATAgcgc < 1:285851/143‑1 (MQ=255) gTTCAGCGGGAAATGCCAGGCATAGAGATGATTTTTCGCCCCCGTCAAACTTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGATGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCgg < 1:506279/143‑1 (MQ=255) cGGGAAATGACAGGCATAGAGATGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGATGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCg > 1:568076/1‑143 (MQ=255) caaaTGACAGGCAGAGAGGTGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGACGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTAt < 2:1801/142‑1 (MQ=255) aaaTGACAGGCATAGAGATGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGACGACAACCAATCGGTAACTGTGCCAGCAGTGGGCTAGCGCCGGGCACCGTCtt > 2:651369/1‑143 (MQ=255) tGACAGGCATAGAGATGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAAACTGGAAACAGAGCGCGACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAg > 1:150324/1‑143 (MQ=255) aTAGAGATGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGACGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCAGCCACCGTATTGAGGCGACttt > 1:117106/1‑143 (MQ=255) gatTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGATGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGc < 1:601236/143‑1 (MQ=255) gatTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGACGACAACCCATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGCCTTTTATGCggg > 1:589840/1‑142 (MQ=255) ttttttCGCCGCCCTCAAGCTTGGCGTCACATTGCATTCTCGTTGTGCGTACGGGCAACGTGGCCCCAGACGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCa < 2:359181/142‑1 (MQ=255) aCGTGGCCTCACCATGCATTCTCGTTGTGCATCCGGGCAACGTGGCCCCAGACGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTg < 1:578263/143‑1 (MQ=255) cgccgCGTCACACTGCAGTCGCGCTGTGCATACGGGCAACGGGGCCCCACACGACAACCAGTCGGTAACTGTTCCAGCAGTGGGAGAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTgg < 1:484092/139‑1 (MQ=255) gTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGATGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGACCACGTTTTTGCGGCATCGCGCTGCCGAAGTCTGGTATcgc > 1:649834/1‑143 (MQ=255) gtGCACTCTCGTTGTGCATACGGGCAACGTGGCCACCGGCGACTGGCTATCGGTAACGGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGAt < 2:637128/142‑1 (MQ=255) gCATTCTCGTGGTGCATACGGGCACCGTGGCCCCAGATGACAACCAATCGGCAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCa < 2:514912/143‑1 (MQ=255) cATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGACGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCGGGCGCAGTGTGGTATCGCGCGGATcgg > 1:250418/1‑141 (MQ=255) tctcGTTGTGCATACGGGCAACGTGGCCCCAGACGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGCATCGAGGCGACTTTTGTGCGGCATCGCGCTGCCGACGTCTGGTATCGCGCGGATgcgcgc > 2:395786/1‑143 (MQ=255) tcGTTGTGCATACGGGCAACGTGGCCCCAGACGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCt < 2:104976/143‑1 (MQ=255) ggggtgCATACGGGCAACGGGGCCCCAGACGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTgg < 2:547600/140‑1 (MQ=255) gggCAACGGGGCCCCAGATCACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGgatg < 1:603389/143‑1 (MQ=255) ggCAACGTGGCCCCAGATGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCCGCGCCTGGTTAGAAtgatga > 1:86831/1‑143 (MQ=255) ccccAGATGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTCAGGCGACTTTTATGCGGCATCGCGCTGCCGCCGTCTGGTATGGCGCGGATCAGCGCCTCGGTATAAGGATGATGTGGCATCg > 2:356236/1‑143 (MQ=255) ccAGACGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGACGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTc > 1:301877/1‑143 (MQ=255) | AGCAGCGTTTCGATCATCTCACTCTTTCTCCATGTTCAGCGGGAAATGACAGGCATAGAGATGATTTTTCGCCCCCGTCAAACGTGGCGTCACAATGCATTCTCGTTGTGCATACGGGCAACGTGGCCCCAGACGACAACCAATCGGTAACTGTTCCAGCAGTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTC > NC_000913/1352618‑1352888 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |