Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,555,706 | (G)6→7 | intergenic (‑37/+120) | maeA ← / ← sra | malate dehydrogenase, NAD‑requiring/30S ribosomal subunit protein S22 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,555,700 | 1 | . | G | 89.5% | 54.7 / 0.7 | 19 | intergenic (‑31/+126) | maeA/sra | malate dehydrogenase, NAD‑requiring/30S ribosomal subunit protein S22 |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base G (10/7); total (11/8) | |||||||||||
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.97e-01 |
TGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑T‑GGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGA > NC_000913/1555573‑1555824 | tGAAGGCACTGCCTTTATTCAAGAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑T‑ggggg < 1:192578/133‑1 (MQ=255) aaGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑T‑ggggggg > 2:263024/1‑132 (MQ=255) tGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑TGGGGGGG‑CTTTTAt > 1:211059/1‑133 (MQ=255) tGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑TGGGGGGG‑CTTTTAt > 1:31515/1‑133 (MQ=255) gCCTTTATGCAACAAAGGCAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑TGGGGGGG‑CTTTTATc < 1:217305/133‑1 (MQ=255) cGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCCAG‑GGGGGG‑ATTTGATAGTCTTTGCTTTACCg > 1:141350/1‑133 (MQ=255) gcagTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑T‑GGGGGGCCTTTTTTGGCCTTTGCTTTACCCCCAggcggccg > 1:438112/1‑127 (MQ=255) gcagTACAGGCCCAGCGTAAGGGATATATAGCGAACGCTGTCTTTTTGTTTTTGGTGCCATGTCCCTCACTCTTTTTTGAATATCCAT‑CCC‑TGGGGGGG‑CTTCTATCGTCTTTGCATTGCCGCCAGGGCGTCg < 2:211618/133‑1 (MQ=255) gTACAGGGCCCGCGTAAGGGATATCAAGCGAACGCTGTTTTCTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGcc < 2:211059/133‑1 (MQ=255) tcAGGGGCCGCGTAAAGGAAAAAAAGGCAACGCCGTTTTTTTTGTTTTGGTTCCATGTCCCTCCCTCTTTTTTGAATAACCCTCCCC‑G‑GGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTc < 2:438112/132‑1 (MQ=255) aCAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTc < 2:54841/133‑1 (MQ=255) aCAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTc < 1:263024/133‑1 (MQ=255) aCAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTc < 1:273827/133‑1 (MQ=255) gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATaaa > 1:348114/1‑133 (MQ=255) atatAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaa < 2:348114/133‑1 (MQ=255) gAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaaaacaaaaac > 1:313229/1‑133 (MQ=255) tctTTTTTGAATATCCAT‑CCC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAgg > 2:414622/1‑133 (MQ=255) ttttGAATATCCAT‑CCC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACGCCATTGCGAAGGTAATAGTGTTTAAACGTCTTTGAGGACTg > 2:160887/1‑133 (MQ=255) ttttGAATATCCAT‑CCC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTg > 1:19654/1‑133 (MQ=255) gAATATCCAT‑CCC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATgg < 2:266204/133‑1 (MQ=255) atccat‑ccC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTAt > 1:386495/1‑133 (MQ=255) atccat‑ccC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGGGTTTAATAGTCTTTGAGGACTTATGGTTAt > 1:150645/1‑133 (MQ=255) ccat‑ccC‑TGGGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGa > 2:121385/1‑133 (MQ=255) | TGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCAT‑CCC‑T‑GGGGGG‑CTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGA > NC_000913/1555573‑1555824 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |