Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,092,428 | G→A | 12.9% | intergenic (+31/+60) | deaD → / ← ygiD | ATP‑dependent DEAD‑box RNA helicase DeaD/4,5‑DOPA‑dioxygenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,092,428 | 0 | G | A | 12.9% | 72.2 / 3.1 | 31 | intergenic (+31/+60) | deaD/ygiD | ATP‑dependent DEAD‑box RNA helicase DeaD/4,5‑DOPA‑dioxygenase |
Reads supporting (aligned to +/- strand): ref base G (13/14); new base A (2/2); total (15/16) | |||||||||||
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.62e-01 |
TCAGCTTGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAA > NC_002947/2092288‑2092577 | tcagcttGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCGAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGAt > 1:425853/1‑150 (MQ=255) ttGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGAtttttt > 1:225002/1‑150 (MQ=255) tCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACAGTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCgg < 1:519646/137‑1 (MQ=255) tCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTAc < 1:504468/150‑1 (MQ=255) tCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAAACCCCGGCTTAGGTCgg > 2:519646/1‑137 (MQ=255) cGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGACCTAAGCCGGGGAtt > 1:84344/1‑133 (MQ=255) cGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGACCTAAGCCGGGGAtt < 2:84344/133‑1 (MQ=255) gACGGCCTCGAGAAGCTGCAGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GGCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCtt > 1:19150/1‑150 (MQ=255) gCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGTTCCTGGCGTAACGAAAAAACCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTcgc < 2:249264/150‑1 (MQ=255) aaGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGc < 2:225002/150‑1 (MQ=255) cAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCTCTTCGCGGGCAAGCCCGCGaaaggg > 1:396213/1‑146 (MQ=255) cAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAaga > 1:39105/1‑148 (MQ=255) cAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAaga < 2:39105/148‑1 (MQ=255) aaGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCTCTTCGCGGGCAAGCCCGCGaaagggc < 2:396213/150‑6 (MQ=255) aaGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCgg > 2:576080/1‑115 (MQ=255) aaGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCgg < 1:576080/115‑1 (MQ=255) ggACTCTCGCGTGGCTGGCAAGCAGCTGACGATGCAGGCG‑TACCGCGAGGATTGATCCTTGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGt < 2:299634/150‑1 (MQ=255) ggctggcAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCTTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGa < 1:187707/150‑1 (MQ=255) tggcAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCTCTTCGCGGGCAAGCCCGCGaaagg < 1:562710/125‑4 (MQ=38) tggcAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCTCTTCGCGGGCAAGCCCGCGaaagg > 2:562710/1‑122 (MQ=38) ggcAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTc < 1:521678/150‑1 (MQ=255) gcAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTcc > 1:176985/1‑150 (MQ=255) gCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTa < 2:425853/150‑1 (MQ=255) gAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTg > 1:313833/1‑150 (MQ=255) aTTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGc > 2:579350/1‑150 (MQ=255) tCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGc < 2:291540/87‑1 (MQ=255) tCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGc > 1:291540/1‑87 (MQ=255) cGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGaa < 1:579350/150‑1 (MQ=255) aaaaaaTCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCACTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCCGCGCGCCCATGGTAAACCCCTTATGCACCACGCCTAACTTGCTc > 2:76483/1‑150 (MQ=255) aaaaaTCCCCGACCTAAGCCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTcc > 1:209774/1‑150 (MQ=255) ccccGACCTAAGCCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGccc < 2:209774/150‑1 (MQ=255) ggTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCaaa < 1:51247/150‑1 (MQ=255) | TCAGCTTGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGC‑GTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAA > NC_002947/2092288‑2092577 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |