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 | . | 84.6% | 12.6 / 6.4 | 14 | intergenic (+55/‑29) | llmg_0683/llmg_pseudo_19 | putative transposase helper protein for IS712B/pseudogene |
Reads supporting (aligned to +/- strand): ref base A (0/1); major base . (7/4); minor base T (0/2); total (7/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.92e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGTGAGAAGAAAACCAAAACCCCAAAAAGTAAATAAAGTGATGGGCCCAAAACCAACAATTTGGGTGGCCCATTTCCTACATTTTCAATGAACCCAAAACCAACATTTTTAGATGGCCCTTGACACTCTACTCTAAATCTAATGAAAGCTCCGTCTCTAAAGCTGATTTTCTAGAATTAAAAAGACAACTGGCTAAAGTTAAGGAAGAACGAGACATCTTAAAAAAAGTATTGACCAT > AM406671/674757‑674994 | agtgttaagacAAACAATACCTCAAACGGTATATAAAGTGATGGGCCCAAAATGACAGATCTGATTAGTCACTTTCGTACATTTTCAGTG‑CCCCAAAATCAACCTTTTTAGATTGCACTTAACACTCTACTCTTAAACTTATGAAAGCt > 2:4715943/7‑149 (MQ=2) gaaagtcaaaaccCCAAAAAGTAAAGAAAGTGATGGGCCCAAAAGCAACAATTTGCGTGGCCCATTTCCTACTTTTTCAATG‑ACCCAAGCCTTACATTTTTAGATGGCCGTTGTCGCTCTTCTCTAAATCTTATGCAAGCTCCGTctct < 2:4252775/143‑1 (MQ=14) acaccACAACAAGTAAATAAAGTGATGGTCCCAAGACCAACAGTTGGGGTGGCCCATCTCCTACACTTTCAATG‑ACCCAAGTCCAACATTTTTAGATGGCCCTTGGCACTCTACTCTAAATCTATTGAAAGCTCCGTCTCTAAAGCTGa < 2:8087765/147‑1 (MQ=17) accCCAAAAGGAAGATAAAGTGATGGGGCCCAAACCCACAACGTGGGGGGCCCAACGCCTACATTTACAATG‑ACCCACAACCAACACGTTTAGAGGGACCTTGACACTCTACTCTAAAACTAATGAAAGCTCCGTCTCTAAAGATGAtt > 2:1121729/1‑149 (MQ=19) aatacaTTGATGTGCTCAAAACCAACAATTTGAGTGGCCCATTTCCTACATTTTTAAT‑TACACAAAACCAACATATTTAGATGGCCCTTGACACTATACTCTCAATCTAATGAAAGCTCCGACTCTAAACCTGATTATCTAGAATTaaa < 2:2142168/149‑1 (MQ=17) aataaaGTGATGGGCCCAAAACCAACAATTTGGGTGGCCCATTTCCTACATTTTCCATG‑ACCCAAAACCAACATTTTTATATGGCCCTTAACACTCTACTCTAAATCTAATGAAATCACCATCCCTAAAGCTGATTCTCTAGAATTaaa > 1:1065087/1‑149 (MQ=17) taaaGTGATGGGCCCAAAACCAACAATTTGGGTAGCCCATTTCCTACATTTTCAAT‑AACCCAAAACCAACATTTTTAGATTGCCCTTGACACTCTACTCTAAATCTAATGAAAACTCCGTCTCTAAAGCTGATTTTCTAGAATTaaaaa < 2:965882/149‑1 (MQ=255) aGTGATTGGAACAACACCAACAATTTGGGTGGCCCATTACCTACATTTTCAATG‑ACACAAAACCAACATTTTTAGATCGACCTTAACACTCCACTCTAAATCTAATGCAAGCTCCGTCCCAAAAGCGGCTTGTCTATAATTAAAAAGAc > 2:4703572/1‑149 (MQ=17) aTGGGCCCAAAACCAACAATTTGGGTGGCCCATTTCCTACATTTTCAATG‑ACCCAAAACCAACATTTTTAGATGGCCCTTGACACTCTACTCTAAATCTAATAAAAACTCCGTCTCTAAAGCTGATTTTCTAGAATTAAAAAGACAACt > 1:6589282/1‑149 (MQ=255) ggCACAAAACAAACAATTTGGGTGGCCGATTTCCTACATTTTCAATG‑ACCCAAAACCAACCGTTTTAGATGTCCCGGGACACACTACACTAAATCTAATGAAAGCTCCGGCTCTAAAGATGATTATCTACAATTAAAAAGACAACAGGc < 2:1649731/149‑1 (MQ=18) gCCCAAAACCAGCAATTTGGGTGGCTCATGGCCTACATTTTCAATG‑ACCCAAAACCACCATTTATAGGTGGCCCTTGACCCTCTACTCTAAATCTAATGAAAGCTCCGCCTCTAAAGCTGATTCTCTAGAAATACAAAGACAAcggtct > 2:1626407/1‑144 (MQ=18) aaacaaacATATTGAGTGGCCCATTACATACTGTTTCAA‑ATACCCAACACCAACCATATTAGATGGCCCATGACACTCTACTCTAAATCTAATAAAAGCTCCGTACATAAAGCTGATTTTCTAGAATTAAAAAGACAACTGGCTAAAGt < 2:6743084/149‑1 (MQ=12) aaTTTCATACATTATCAATG‑AGCCAAAACAAACATTCTAAGAGGGCCCTTAACACTCTACTCTAAATCTAAACAAAGCTCCGCCTCTAAAGCGAAGATACTACAAATAAAAAGACAACTGGCTAAAGTTAAGGAAGAACGAGACATCtt < 2:7633358/148‑1 (MQ=12) aTTTCCTACATTTTGAATG‑ACCCAAAACCAACATTTTTAGATGGCGCTTGACTCTCTACTCTAAATCTTATGAAGGCTCCGTCTCTAAAGCTGTTTTTCTAGAATTAAAAAGCCAACTGACTAAAGTTAAGGAGGAAAGAGACCTCTTa > 2:10064745/1‑149 (MQ=25) aaaCCCAAAACCAACATTTTTAGATGGCCCATGACACTCTACTCTAAATCTAATGAAAGCTCCGTCTCTAAAGCTGATTTTCTAGAATTAAAAAGACAACTGGCTAAAGTTAAGGAAGAACGAAAAAAAAAAAAAAAAGTATTGACCat < 2:9324704/148‑1 (MQ=17) | AGTGAGAAGAAAACCAAAACCCCAAAAAGTAAATAAAGTGATGGGCCCAAAACCAACAATTTGGGTGGCCCATTTCCTACATTTTCAATGAACCCAAAACCAACATTTTTAGATGGCCCTTGACACTCTACTCTAAATCTAATGAAAGCTCCGTCTCTAAAGCTGATTTTCTAGAATTAAAAAGACAACTGGCTAAAGTTAAGGAAGAACGAGACATCTTAAAAAAAGTATTGACCAT > AM406671/674757‑674994 |
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: ‑ |