Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | AE015451 | 5,674,751 | +G | intergenic (+113/+27) | rhlE‑II → / ← yceI | ATP‑dependent RNA helicase/polyisoprene‑binding protein, acidic stress response factor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 5,674,750 | 1 | . | G | 96.5% | 95.3 / ‑2.9 | 29 | intergenic (+112/+28) | rhlE‑II/yceI | ATP‑dependent RNA helicase/polyisoprene‑binding protein, acidic stress response factor |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (17/11); total (18/11) | |||||||||||
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.40e-01 |
CCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑GCG‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCAC > AE015451/5674554‑5674951 | ccGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCtt > 1:309490/1‑237 (MQ=255) ccTGAGCAACTGGATGACCTGCCCAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCACCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTa < 1:299449/215‑1 (MQ=255) ccACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGat < 1:56417/212‑1 (MQ=255) ctgctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACAACTTCGACACTGATAtcca > 1:242010/1‑196 (MQ=255) gggCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGat > 1:26593/1‑177 (MQ=255) gggCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGat < 2:26593/177‑1 (MQ=255) gATAGAAGTCCCGAGGGGGCTGCTTCGCAGCCCTAGCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGaca < 2:82685/167‑1 (MQ=255) ccGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGt > 1:35253/1‑212 (MQ=255) gctTCTCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGTCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttca > 1:434948/1‑207 (MQ=255) caAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGcc < 2:242008/174‑1 (MQ=255) ccGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCACCGCGCTCGCCACCCCATGGGTCCt > 1:145764/1‑234 (MQ=255) cGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGT‑CAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGcc > 2:434367/1‑236 (MQ=255) cTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCANCCGGGCCTTTg > 1:3046/1‑158 (MQ=255) cTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgtt > 2:154256/1‑172 (MQ=255) aTTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAgg > 2:214875/1‑237 (MQ=255) ttCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTa > 1:391184/1‑237 (MQ=255) gATCTCTACAGGAGCGGCCTCGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCATGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtg > 2:441248/1‑172 (MQ=255) ctctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCt < 1:73867/116‑1 (MQ=255) cAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGT‑CAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCAc < 1:434369/237‑1 (MQ=255) ccTTGTGTCGTGATAGGGCCGCAAAGC‑GGGCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCt > 2:73866/1‑66 (MQ=37) ttGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCg > 2:126707/1‑103 (MQ=255) ttGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATgg > 1:58909/1‑188 (MQ=255) ttGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTg > 2:217620/1‑207 (MQ=255) ttGTGTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGACGGTAGc > 2:241177/1‑215 (MQ=255) tgtgTCGTGATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCt < 1:55157/197‑1 (MQ=255) gATAGGGCCGCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCATGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCAc < 1:441250/218‑1 (MQ=255) gCAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTTCTTCTGCTTCGCACCTGCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCt < 2:58908/157‑1 (MQ=255) cAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccag > 2:132754/1‑76 (MQ=255) cAAAGCGGCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccag < 1:132755/76‑1 (MQ=255) | CCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑GCG‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCAC > AE015451/5674554‑5674951 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |