Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,845,703:1 | +G | 100% | intergenic (‑70/+28) | PP_3394 ← / ← PP_3395 | 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,845,703 | 1 | . | G | 95.8% | 74.7 / ‑2.8 | 24 | intergenic (‑70/+28) | PP_3394/PP_3395 | 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (14/9); total (14/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.17e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
CTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGC > NC_002947/3845558‑3845842 | ctgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑ggg > 6:243762/1‑148 (MQ=255) gctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑gggct < 2:229305/149‑4 (MQ=255) tCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGAt > 3:176368/1‑148 (MQ=255) cATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATc < 5:35174/149‑1 (MQ=255) gATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGt > 8:188222/1‑148 (MQ=255) aTCACTTGCAGGCTCTGCTAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTc < 4:32250/148‑1 (MQ=255) tCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCa < 7:53594/148‑1 (MQ=255) aCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAggc > 2:144575/1‑147 (MQ=255) ggCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATgg < 8:250865/149‑1 (MQ=255) gCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGc > 1:80457/1‑149 (MQ=255) ttgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAg < 5:309013/122‑1 (MQ=255) ttgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAg > 6:309013/1‑122 (MQ=255) tattattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTTATACATGGTGAGCCACCATGGATGAGCCCCCCCCCTGCCCCGCCGGGGGCCgcagca > 2:149177/3‑148 (MQ=255) gtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCCGCACCAggtggtt > 5:17824/1‑147 (MQ=255) ctctGTCTTAAGCGTTGCGCCAGAATAATAGTTGTATGAGGGGAAGAAAAGACGTGTTTCGGG‑GGCTTATCCATCGTCAGCCACGTTGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccga < 2:80457/148‑1 (MQ=255) gCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc > 3:278246/1‑146 (MQ=255) gcCATACTAAGAGCTGTATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGAACCATCGTCACCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtgt < 7:34859/149‑1 (MQ=255) cAGACTAAGAGCTGGAGGACGGGCAGAAACGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCACCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGGCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGc < 7:210954/149‑1 (MQ=255) agagCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgt > 7:284611/1‑138 (MQ=255) gagCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgcg > 8:219816/1‑149 (MQ=255) ggATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGc < 5:243762/149‑1 (MQ=255) gATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg > 2:302363/1‑147 (MQ=255) gATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATACATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg > 2:297573/1‑147 (MQ=255) gACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTGAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGAt > 3:30895/1‑149 (MQ=255) cAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGCCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTATGCCGGTCGCGCATGCAGATGACGaa < 7:188222/149‑1 (MQ=255) gtgtTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgc > 7:42866/1‑149 (MQ=255) | CTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGC > NC_002947/3845558‑3845842 |
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: ‑ |