| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NZ_CP009273 | = 2000783 | 9 (0.620) | 3 (0.230) | 3/154 | 1.3 | 24.8% | coding (1147/1488 nt) | amyA | alpha‑amylase |
| ? | NZ_CP009273 | = 2000836 | 10 (0.760) | coding (1200/1488 nt) | amyA | alpha‑amylase | |||||
AATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACAGACGTT > NZ_CP009273/2000698‑2000871 | aaTTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAAcc < 1:42838/91‑1 (MQ=255) ttttGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAACCTa < 2:276319/91‑1 (MQ=255) cGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGGAAACCTATCCGATAGATATGCCaata > 1:261268/1‑91 (MQ=255) gACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCgagaattaactcatcaagctgttcgattattggcatatctatcggata > 2:306580‑M1/1‑43 (MQ=255) gcgcATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCgatagata > 1:113543/1‑49 (MQ=255) gcgcATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCgatagata < 2:113543/49‑1 (MQ=255) gcgcATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCgcccgtc < 2:55344‑M2/91‑8 (MQ=255) gcATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCgcccgtcac > 1:159967‑M2/1‑82 (MQ=255) cATTACGAAGATGTCGGTGGTGACGGGCgagaattaactcatcaagctgttcgattattggcatatctatcggataggtttgcccgtcacc > 2:173522‑M1/1‑28 (MQ=255) tGTCGGTGGTGACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCAc < 1:55344/91‑1 (MQ=255) ggtgACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCgcccgtcaccaccgacatcttcgtaatg > 2:173522‑M2/1‑63 (MQ=255) gtgACGGGCgagaattaactcatcaagctgttcgattattggcatatctatcggataggtttgcccgtcaccaccgacatcttcgtaatgc > 1:159967‑M1/1‑9 (MQ=255) gtgACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACa < 1:129356/91‑1 (MQ=255) gACGGGCgagaattaactcatcaagctgttcgattattggcatatctatcggataggtttgcccgtcaccaccgacatcttcgtaatgcgc < 2:55344‑M1/91‑85 (MQ=255) ggCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACAGACGtt < 1:182601/91‑1 (MQ=255) ggCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACAGACGtt < 2:334806/91‑1 (MQ=255) | AATTTTGTTGCGGGAAAATGGCGTTCCTTCGGTATTCTATCCGGACCTCTACGGTGCGCATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACAGACGTT > NZ_CP009273/2000698‑2000871 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ 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. |