Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,092,428 | G→A | 15.8% | 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 | 15.8% | 73.3 / 5.9 | 33 | intergenic (+31/+60) | deaD/ygiD | ATP‑dependent DEAD‑box RNA helicase DeaD/4,5‑DOPA‑dioxygenase |
Reads supporting (aligned to +/- strand): ref base G (14/13); new base A (3/2); total (18/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCTTGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGC > NC_002947/2092291‑2092565 | gcttGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGAtttt > 2:502658/1‑150 (MQ=255) cttGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGAttttt > 2:252795/1‑150 (MQ=255) tGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGAttttttt > 2:249730/1‑150 (MQ=255) tGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGAttttttc < 1:485757/150‑2 (MQ=255) gcatCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTcgc < 1:502658/147‑1 (MQ=255) ccTCGAGAAGCTGCTGTCCCAGCTCAAGGAAACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGCCTTAGGTCGGGGGTTTTTTTTACCTGTACCGGCCCTTTcccg > 2:154708/1‑150 (MQ=255) tCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCggg > 2:558459/1‑150 (MQ=255) tgTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGCTCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGaa < 1:610971/150‑1 (MQ=255) aaGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGACCTAAGCCGGGGAttttttc > 1:383770/1‑108 (MQ=38) aaGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGACCTAAGCCGGGGAttttttc < 2:383770/109‑2 (MQ=38) aGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGcc < 1:252795/150‑1 (MQ=255) aGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATAAACCACATTCGTCGGGGATTTTTTTTACCTGTACCGGCACTTTCGCGGGCGTGCCCGCGAAGAGGCCCTTCAg > 1:13404/1‑150 (MQ=255) acTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTAc > 1:598960/1‑100 (MQ=255) acTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTAc < 2:598960/100‑1 (MQ=255) tCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCa > 1:220465/1‑142 (MQ=255) tCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCa < 2:220465/142‑1 (MQ=255) aaGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCtt < 2:96018/127‑1 (MQ=255) aaGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCtt > 1:96018/1‑127 (MQ=255) aGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAgcgcg < 1:249730/150‑1 (MQ=255) tGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCa < 1:660894/142‑1 (MQ=255) tGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCa > 2:660894/1‑142 (MQ=255) tACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATc < 1:440916/130‑1 (MQ=255) tACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATc > 2:440916/1‑130 (MQ=255) gcgAGGATTGCTCCTGGCGTAACGAAAAAACCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCACTTTCGCGGGCTTGACCGCGAAGAGGCCCGTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAAcccc < 1:154708/150‑1 (MQ=255) tGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTTAGTCGAACCGATAGATGTCCAGCGCGAGTGCGCCCAGGGTAAAACCCTGATGCACCACGcc > 2:552360/1‑150 (MQ=255) cGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAAccc < 2:306330/132‑1 (MQ=255) cGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAAccc > 1:306330/1‑132 (MQ=255) aaaaaTCCCCGACCTAAGCCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCa < 2:108466/105‑1 (MQ=38) aaaaaTCCCCGACCTAAGCCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCa > 1:108466/1‑105 (MQ=38) aaaaTCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTccc > 1:562855/1‑150 (MQ=255) aTCCCCGACCTAAGCCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTACCTGCAAAGTTGCCCTTAAGAGGATGCTGGAGATTGCGATCCCCAGCGCGCGCAGGGTAAACCCCTGCTGCACCACGCCGAACTTGCTCCCgcc > 1:279522/1‑148 (MQ=255) cccGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAg > 1:86119/1‑108 (MQ=255) cccGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAg < 2:86119/108‑1 (MQ=255) | GCTTGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGC > NC_002947/2092291‑2092565 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |