Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | NC_000913 | 1,396,311 | 0 | A→C | 100.0% | 1.5 | 4 | Q79P (CAG→CCG) | insH1 | IS5 transposase and trans‑activator |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (2/0): ref base (0/1): total (3/1) |
GGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTACAACCTGAGCGATGGCG > NC_000913/1396018‑1396337 | gggAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCCTGCTACGCCTTCCCTGCATGCCCCATTgg > 2:107116/1‑301 (MQ=255) aaggaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGATATCATGTTTTTCATCGGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACGTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTa < 2:152033/300‑1 (MQ=25) aCATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCACGGCTGGTAATGGCCGGCGACCTTATCCCATGGAAACCATGCTACCCATTCACTGCAGCCGGCATTTGTCCACccgaagagag > 1:230676/1‑293 (MQ=11) aGCTAAGGAAGGTGCGAACAAGTCACTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCCACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGCAAAGAGAGTTTCTTGTCGCGCCTGGAGCGGATTCTTCCATGGAAAAACATGGTGGGCGTCATCGAGCCGTGTTACCCCAAGGCTGGTAATGGCGGCCGACCTTATCCGCTGGCACCCCTGTCACGCCTTCACTCCATGCCGCGTTGGTCAAGACAGAGCGATGccg > 2:360760/1‑300 (MQ=11) | GGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTACAACCTGAGCGATGGCG > NC_000913/1396018‑1396337 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
Reads not counted as support for junction |
read_name Not counted due to insufficient overlap past the breakpoint. |
read_name Not counted due to not crossing MOB target site duplication. |