Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | AM406671 | 674,848 | Δ1 bp | intergenic (+56/‑28) | llmg_0683 → / → llmg_pseudo_19 | putative transposase helper protein for IS712B/pseudogene |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AM406671 | 674,847 | 0 | A | . | 80.0% | 18.4 / 12.4 | 20 | intergenic (+55/‑29) | llmg_0683/llmg_pseudo_19 | putative transposase helper protein for IS712B/pseudogene |
Reads supporting (aligned to +/- strand): ref base A (0/4); new base . (11/5); total (11/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.60e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTATAGTGAGAAGAAAACCAAAACCCCAAAAAGTAAATAAAGTGATGGGCCCAAAACCAACAATTTGGGTGGCCCATTTCCTACATTTTCAATGAACCCAAAACCAAC‑ATTTTTAGATGGCCCTTGACACTCTACTCTAAATCTAATGAAAGCTCCGTCTCTAAAGCTGATTTTCTAGAATTAAAAAGACAACTGGCTAAAGTTAAGGAAGAACGAGACATCTTAAAAAAAGTATTGACCA > AM406671/674753‑674993 | cTATAGTGAGAAGAAAACCAAAACCCAAAAAAGTAAATAAAGTGATGGGCCCAAAACCAACAAATTGGGTGGACCATTTCCTAAATTATAAATG‑ACACAAAACCAAC‑ATTATCAGATGGCTATAGACACGATACTCTAAATATAATGaa > 2:2549902/1‑149 (MQ=11) atataGTGAGAAGAAAACCAAAACCCCAAAAGGTAGATGAAGTCATGGGCCCAAAACCAACGATGGGGGTGGCGCGTTTCATACGTTTTCACTG‑ACCATAAAGCAAT‑ATTTTGAGATCGCCCTTGCCACTCTACTTTAAATCTATTGaa < 2:5387414/148‑1 (MQ=2) agaagaGAGCCAAAACCCCAAAAAGTAAATAAAGGGATGGGCCCAAAACCCACAACTTGCGTGGACAATTACCTACATTTTCAATG‑ACCCAAAACCAAC‑ACTTTTAGAAGGCCCTCGACACCCTACTCTACATCTACTGAGAGCTCCgc > 2:7015608/1‑148 (MQ=11) agaaAACGCAAAAACCAAGAAGGAAATATAGTGATGGGCCCAAAACCAACAATCTGCGTGGCCCATTTCCTACATTTTCAATG‑ACCCAAAACCAAC‑ATCTTTAGAGGGCCCTTGGAACACTACTATAAATCTAATGAAAGCTCCGtctc < 2:9195869/149‑1 (MQ=14) aaaccCCAAAAAGTAAATAAAGTGATAGGCCAAAAACCAACAAGTGGGGTGGCACATTTCCTACATTTTCAATG‑ACCCAAAACCAAA‑ATGTTTAGATGGCCCTCGAGACTCTACTCTAAAGCTAATGAAAGCTCCAGCTCTAAAGCTGa > 2:503668/1‑149 (MQ=17) ccccAAAAAGTAAATAAAGTGATGTGCCCAAAACCAACAATTTGGTAGGCCCAGAAACTACATTTTAAATG‑ACCCTAAAAAAAC‑ATTATTTGATGGCCCTTGAAAATCTACTCTAAATATAATGAAATCTACGACTATAAAGCTGAtct > 2:8266879/1‑147 (MQ=9) aGTAAATAAAGTGATGGGCACAAAACCAACAATTTGGGTGGCGCATTTCGTACATTTTCAATG‑ACCCAAAACGAAC‑ATTTTTAGATGGCCCTTGACACTCTACTATAAATCTAATGAAAGCTCCGACTCTAAAGCGGATTTTCTAGAAt > 2:6279105/1‑149 (MQ=18) aGTGATGGGCCCAAAACCAACAATTTGGGTGGCCCATTTCCTACATTTTCAATG‑ACCCAAAACCAACAATTTATAAATGG‑‑ATAGAC‑CTCTACTCTAAATCTAATGAAAACTCCGTCTCTAAAGCTGATTTTCTAGAATTAAAAAGACaa > 2:2470571/1‑149 (MQ=14) tGATAGGCGCAGAAGCAGCAATTGGAGTGGACGAGTTCCTACGTTTACAAT‑AACCCAAAACCGAC‑ATTTTTAGATGGCGCTTGACACTCTACGCTAAATCTAATGAAAGCTCCGTCTCTAACGCGGATTTTCTAGAATTAAAAAGACaa < 2:5104254/149‑1 (MQ=18) aTGGGCCCAAAACCAACAATTTGGTTGGCCCATTTCCTACATTTTCAATG‑ACCCAAAACCAAC‑ATTTTTAGATGGCCCTTGACACTCTACTCTAAATCTAATGACAACTCCGTCTCTAAAGCTGATTTTCTAGAATAAAAAAGACAAcc > 2:8922134/1‑148 (MQ=21) aTGGGCCCAAAACCAACAATTTGGGTGGCCCATTTCCTACATTTTCAATG‑ACCCAAAACCAAC‑ATTTTTAGATGGCCCTTGACACTCTACTCTAAATCTAATGAAAACTCCGTCTCTAAAGCTGATTTTCTAGAATTAAAAAGACAACt > 2:8922171/1‑149 (MQ=255) aCCAACAATTTGGGGGGCCCATATCCTACATTTTCAATG‑ACCCAAAACCAAC‑ATGTTTAGATGGCCCTTGACACTCCACTCTAAATCTAATGAAAGCTCCGTCTCTAAAACCGCTTTTCTAGAATTAAAAAGACAACAGGCGAAAGTTa > 2:5097550/1‑149 (MQ=21) aaTTTGGATTGCCCATTAACGACTTTTTCAATG‑ACCCAAAACCAAC‑AGTTTTATATAGCCTTAAACACTCTACTATAAAACTAAAGAAAACTACGTCTCTAAAGCAGATTTTCTAGAATTAAAAAGTCAACTGGCAAAAGTAAAAgaaa < 2:2391245/149‑2 (MQ=12) ctctttACTACGTTTGCAATG‑ATCCAACACCAAC‑ATGTTTAGATGACCTTAGACTCCCTACTCTAAATCTAATGAAAGCTCCGTCCCTAAAGCTGCGTTACTGGAAGTAACAAGACAGCTGGCTAAAGTGAAGGGAGCACGAGACATCt < 2:9082386/146‑1 (MQ=9) cATTTGCTACATTTTCAATG‑ACCCACACCCAAC‑ATTTTTTGATGGCCCTTGACACACTACTCTATATCTGATGAAGGCTCCGTCTCTAATGCTGATTTACTAGACTTAAAAAGACAACTTGCTAAAGTTAAGGAAGAACGAGACATCtt > 2:2693813/1‑149 (MQ=18) gtgtCCTACATTTGCAAAG‑ACCCAAACCAACC‑ATTGTTAGATGGCCCTCGACCCTCTACGCTAAATCTAATGAGATATCCGTCTCTAAAGGTCTTTCTGTTGAATTAAAAAGACAACTGGCTAAAGTTAAGGAAGACCGAGTCATCTTa < 2:5484729/146‑1 (MQ=9) ttCAATGAA‑ACAAAACCAAC‑ATTTTTAGACGGCCCTAGACAATATAATATAAATATAATGAAAGCTGAGTCTCTAAAGCTGATTTTCTAGAATTACAGAGACAACTGGCTAACGTTAATGAAGTACGAAAAAAAAAAAAAAAAGTATTg < 2:6694670/149‑1 (MQ=12) ttCAATG‑ACCCAAAACCAAC‑ATTTTTTTTTTTTTTTTTACACTCTACTCTAAATCTAATGAAAGCTCCGTGTCTAAAGCTGATTTTCTAGAATTAACAAGACTACTTTCTAAAGTTAAGAAAAAACGAGACATTTTAATAATAGTAttt > 2:158209/1‑148 (MQ=2) aaaTGAA‑CCAAAACCAAT‑AATTTTAGATAGCCCTTGAAACACTACTCTAAATCTAATGAAAACTCCGTCTCTAAAGATAATTTTCTAGAATTAAAAAGACAACTTGATAAAGTTAATGAAGTACGAAACATATTAAAAAAAGTATTGAc < 2:10805377/148‑1 (MQ=14) aTGAA‑CCAAAACCAAC‑ATTTTTAGATAGCACTAGACACAATACACTAAATCTAATGAAAGCTCCGTCACTAAAGCTGATTTTCTAGAATTAAAAAGAAAACTGGTTAAAGTTAAGGAAGAACGAAACAAAAAAAAAAAAGTATTGACCa < 2:10151976/149‑1 (MQ=14) | CTATAGTGAGAAGAAAACCAAAACCCCAAAAAGTAAATAAAGTGATGGGCCCAAAACCAACAATTTGGGTGGCCCATTTCCTACATTTTCAATGAACCCAAAACCAAC‑ATTTTTAGATGGCCCTTGACACTCTACTCTAAATCTAATGAAAGCTCCGTCTCTAAAGCTGATTTTCTAGAATTAAAAAGACAACTGGCTAAAGTTAAGGAAGAACGAGACATCTTAAAAAAAGTATTGACCA > AM406671/674753‑674993 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |