| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913 | = 1978502 | 0 (0.000) | 25 (0.260) | 20/288 | 4.0 | 100% | intergenic (‑305/+16) | flhD/insB1 | flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B |
| ? | NC_000913 | 1979279 = | 0 (0.000) | intergenic (‑64/‑474) | insA/uspC | IS1 repressor TnpA/universal stress protein | |||||
| Rejected: Position hash score below cutoff. | |||||||||||
AGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAA > NC_000913/1978306‑1978516 | aGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTGAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGAAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATgtat < 2:186497/200‑4 (MQ=255) tGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATgtataag < 2:785895/200‑7 (MQ=255) aGAAAT‑GGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTCCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATgtat < 1:624517/150‑4 (MQ=255) aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATgtataag < 1:684116/150‑7 (MQ=255) tGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATgtataagt < 1:658617/150‑8 (MQ=255) cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATgtataagtc < 3:287079/141‑9 (MQ=255) tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATgtataagtcatactt < 3:265136/141‑15 (MQ=255) cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATgtataagt < 1:88732/88‑8 (MQ=255) cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATgtataagt > 2:88739/1‑81 (MQ=255) aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATgtat > 3:705638/1‑57 (MQ=255) aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATgtat < 4:705608/60‑4 (MQ=255) | AGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAA > NC_000913/1978306‑1978516 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ 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. |