Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | NC_000913 | 274,861 | 0 | T→G | 100.0% | 0.7 | 4 | D86A (GAT→GCT) | insH1 | IS5 transposase and trans‑activator |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (0/2): ref base (1/0): total (1/3) |
ATCGCGCCATCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTT > NC_000913/274852‑275158 | ctCGCGCCAGCGCTCAGGGTGTAACAATGCTGCATGCAGGTAATGGTTAGCATGGGTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCGTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGa < 2:479725/300‑1 (MQ=2) cgcgCCCGCGCTTAGGGTGTAACAATGCGGCATGCAGTGAATGCGTAGCATGGTCTCCAGCGGATAAGGTCGCCGGACATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTGGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAgg < 2:410515/301‑1 (MQ=255) cgcgCCACCGCGCAGGTTGTACGAAGGCTGCATGAAGTGAATGCGTAGCATGGTGTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTGCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAgg < 2:331335/301‑1 (MQ=9) ccATCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTTGTCTGACGGCGCTTCCTGCTGAATTAACTGTCGGCGAAGGTAAGTTGATGACTCCTGATGAACCCTGTTCTATGGCTCAAGATGACAAAAATAATATTATATACGGGACTTTTTCGCAACTTCCGGAGGCTtt > 1:588206/1‑301 (MQ=11) | ATCGCGCCATCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTT > NC_000913/274852‑275158 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 27 ≤ 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. |