Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | NC_000913 | 331,807 | 0 | T→A | 100.0% | 7.8 | 2 | intergenic (+311/‑564) | betT/yahA | choline transporter of high affinity/c‑di‑GMP‑specific phosphodiesterase |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (1/1): ref base (0/0): total (1/1) |
GTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGA > NC_000913/331685‑331928 | gTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGAGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGa > 1:504854/1‑244 (MQ=255) gTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGAGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGa < 2:504868/244‑1 (MQ=255) | GTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGA > NC_000913/331685‑331928 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ 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. |