Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 3,845,703:1 | +G | 100% | intergenic (‑70/+28) | PP_3394 ← / ← PP_3395 | putative 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/Transcriptional regulator, LysR family |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 3,845,703 | 1 | . | G | 97.4% | 133.6 / ‑3.0 | 38 | intergenic (‑70/+28) | PP_3394/PP_3395 | putative 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/Transcriptional regulator, LysR family |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (23/14); total (23/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.95e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-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. |
GCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑G‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACAC > AE015451/3845471‑3845919 | gCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTCCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑g‑g > 1:785705/1‑235 (MQ=255) cACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGc > 1:59185/1‑236 (MQ=255) gCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGCAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACg > 2:353989/1‑207 (MQ=255) gCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGc > 2:501631/1‑233 (MQ=255) cGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATgg < 2:517097/209‑1 (MQ=255) cGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTg < 1:780541/236‑1 (MQ=255) cgcCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCa < 1:531046/210‑1 (MQ=255) ggcgTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCcg > 1:37721/1‑205 (MQ=255) aTCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTg < 1:119769/205‑1 (MQ=255) ctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTc < 1:402778/159‑1 (MQ=255) gCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGcc < 1:93221/182‑1 (MQ=255) gATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCg > 2:634525/1‑184 (MQ=255) cAGGCTCTGCCAGCCCTCGCGCAGGCCGGCTTCATTGATTGTTCTGGTTGTCATCGCGGTCACTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGGCGGGCAGAAAAGACGTGTTTggc < 2:319300/124‑1 (MQ=255) aaGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg > 2:502980/1‑233 (MQ=255) gcgcAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCa > 2:123304/1‑184 (MQ=255) gcAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGCTCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCccggccgg > 2:32059/1‑156 (MQ=255) tgttgtCATCGCGGTCTCTCTGTCTGAAGCGGTGCGCCGGACTAAGAGCGGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGAGGGATCAGCCCGCGCCCTGCCAGGCCGGTGc < 2:228758/134‑1 (MQ=255) gttgtCATTGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACGGCGAGAGCGGGTGCCg > 1:29620/1‑187 (MQ=255) ttgtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAg > 2:644806/1‑156 (MQ=255) tCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCGGCCACGAt < 1:476569/99‑1 (MQ=255) tCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGAt > 2:476567/1‑99 (MQ=255) aTCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAg > 1:502526/1‑218 (MQ=255) cgcgGTCTCGCTGTCTGAGGCGTTGCGCCAGACTAAGAGCTGGATGCCGGGCAGAAAAGACGTGTTTCGCGG‑GCGGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTg < 2:270139/124‑1 (MQ=255) gcgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGcc > 1:478833/1‑114 (MQ=255) gcgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGcc < 2:478831/114‑1 (MQ=255) gtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGa < 1:613135/235‑1 (MQ=255) ttGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGAccca > 2:276931/1‑216 (MQ=255) gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGGTCAGCCCGCGCCCTGCCAGGCCGGTGCAc > 2:695518/1‑101 (MQ=255) gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcag > 1:422054/1‑197 (MQ=255) gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCCATGACCCGCAGGTCGAACAAcgccgc > 1:730024/1‑231 (MQ=255) ccAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAACCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGAc < 2:761972/209‑1 (MQ=255) gACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GTGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTg < 1:738497/94‑1 (MQ=255) gACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCa > 2:799420/1‑99 (MQ=255) gCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCCGATGACGAATTCACCAAGCGCCCAgtc > 2:178508/1‑168 (MQ=255) gCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCaa > 2:329924/1‑158 (MQ=255) gCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGa > 1:697404/1‑233 (MQ=255) gCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAaca > 1:717061/1‑236 (MQ=255) gCAGAAAAGACGTGTTTAGCGG‑GCTGATCCATCGTCAGCCACGATGGATCACCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCCGGCGCCGAGCCGACAGCGAGAGCgt > 1:256551/1‑123 (MQ=255) cAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAacac < 1:364454/236‑1 (MQ=255) tttCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGc < 1:36250/176‑1 (MQ=255) | GCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑G‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCTCTTGCGGGAACAC > AE015451/3845471‑3845919 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |