Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,003,617 | C→G | 17.4% | Q452E (CAA→GAA) | PP_4409 → | site‑specific recombinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,003,617 | 0 | C | G | 17.4% | 47.4 / 3.3 | 22 | Q452E (CAA→GAA) | PP_4409 | site‑specific recombinase |
Reads supporting (aligned to +/- strand): ref base C (7/11); new base G (2/2); total (9/13) | |||||||||||
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 |
ATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAATCAACAGAGGCATTAATCATGATAAACATGACCAAGGAGCTTATCGTCCTTGCCAGCAGCAGTCAATTATACGTACTGGTGCGTCAA > NC_002947/5003471‑5003765 | aTAGCGAATTGCTATTGATTACTGGCCACGCTAAAATCAGAAGTCTTCAGAAATACGCGAAGATCGACCGAGAATTATTGTTGGAAGAGATCGCCTACGGAAACACCATAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCaaa > 2:32016/1‑150 (MQ=255) ttGATTACAGGCCATGCTGATATCAGAAGTCTTCAGAAAAACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGTCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTgaa < 2:132980/150‑1 (MQ=255) aTTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTgaagaa < 2:34607/150‑1 (MQ=255) cTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTGAAATTCAGtgcttg < 2:537936/140‑6 (MQ=255) cTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTGAAATTCAGtgcttg > 1:537936/1‑135 (MQ=255) cTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGcc > 2:134103/1‑150 (MQ=255) cAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAg < 2:107369/150‑1 (MQ=255) aGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAATc < 1:326130/150‑1 (MQ=255) gagaATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAATCAACAGAGGc < 1:190868/150‑1 (MQ=255) agaATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAATCAACAGAGGCa < 2:124992/150‑1 (MQ=255) ttgGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGc > 1:399764/1‑90 (MQ=255) ttgGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGc < 2:399764/90‑1 (MQ=255) aaaaTGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAAt > 1:400191/1‑95 (MQ=255) aaaaTGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAAt < 2:400191/95‑1 (MQ=255) atTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAATCAACAGAGGCATTAATCATGATAAACATGACCAAGGAGCTTATCGTCCTTGc < 1:75551/136‑1 (MQ=255) atTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAATCAACAGAGGCATTAATCATGATAAACATGACCAAGGAGCTTATCGTCCTTGc > 2:75551/1‑136 (MQ=255) ttAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAATCAACAGAGGCATTAATCATGa > 1:556583/1‑104 (MQ=255) ttAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAATCAACAGAGGCATTAATCATGa < 2:556583/104‑1 (MQ=255) ggtaCTGAATTTGAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAATCAACAGAGGCATTAATCATGATAAACATGACCAAGGAGCTTATCGTCCTTGCca < 2:267616/126‑1 (MQ=255) ggtaCTGAATTTGAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAATCAACAGAGGCATTAATCATGATAAACATGACCAAGGAGCTTATCGTCCTTGCca > 1:267616/4‑129 (MQ=255) gAATTTCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAATCAACAGAGGCATTAATCATGATAAACATGACCAAGGAGCTTATCGTCCTTGCCAGCAGCAGTCAATTATACGTACTGGTGc < 1:134103/150‑1 (MQ=255) tCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAATCAACAGAGGCATTAATCATGATAAACATGACCAAGGAGCTTATCGTCCTTGCCAGCAGCAGTCAATTATACGTACTGGTGCGTCaa > 1:701516/1‑150 (MQ=255) | ATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGGCAATGTTCAAAGAAAAACTAAAGCAACAAAACAATCAACAGAGGCATTAATCATGATAAACATGACCAAGGAGCTTATCGTCCTTGCCAGCAGCAGTCAATTATACGTACTGGTGCGTCAA > NC_002947/5003471‑5003765 |
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: ‑ |