New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 = 516767043 (0.980)3 (0.070) 3/322 NT 12.6% coding (210/3906 nt) hrpA ATP‑dependent helicase HrpA
?AE015451 = 5168146 1 (0.020)intergenic (‑267/‑237) hrpA/PP_4547 ATP‑dependent helicase HrpA/putative Glutamine synthetase

TTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGGAAGCGGCCTTGTGTCGCGATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCTTACGCCGCATGTGAGTATCAAAGGCATGACTCCCGCCCAGCCCATCGCCCCACCGTCGTCCCGCCCTTCGCTACCGAAGGCCCGGAGCCGTGCCGGGGAAAATCCGCTGGTCCACATCATCCTCGCCTTCTGGGCCTTGTGGCACTGCCGTCACGCGCGCCCACCGGACACTTCGCTCACGCCATTGTGACCAAACCCGGCCAGCGTCTGGCCGTTTGACTCAACCGGGCCG  >  AE015451/5167926‑5168339
                                                                                                                                                                                                                            |                                                                                                                                                                                                 
ttCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGGAAGCGGCCTTGTGTCGCGATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCTTACGCCGCATGTGAGTATCAAAGGCATGACTCCCgccca                                                                                                                                                                                                   <  2:657621/221‑1 (MQ=255)
                                                                                                                  ggtgCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCTTACGCCGCATGTGAGTATCAAAGGCATGACTCCCGCCCAgcccatcgccgccaagcgtgacgaaatcaagaaggccctg                                                                                                                                                           <  1:739108‑M2/147‑41 (MQ=255)
                                                                                                                  ggtgCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCTTACGCCGCATGTGAGTATCAAAGGCATGACTCCCGCCCAgcccatcgccgccaagcgtgacgaaatcaagaaggccctg                                                                                                                                                           >  2:739106‑M2/1‑107 (MQ=255)
                                                                                                                                                                                           cgcATGTGAGTATCAAAGGCATGACTCCCGCCCAgcccatcgccgccaagcgtgacgaaatcaagaaggccctggccgaaaaccaggtattggtgatcgccggcgaaaccggctctggcaagacgacccagttgccgaaaatctgcctggagcttggccgtggcagccatgggttgatcgcccatacccagccacgcaggatcgcagcgcgtagtgtggccgcgcga  >  2:309116‑M2/1‑34 (MQ=255)
                                                                                                                                                                                                                    tCCCGCCCAGCCCATCGCCCCACCGTCGTCCCGCCCTTCGCTACCGAAGGCCCGGAGCCGTGCCGGGGAAAATCCGCTGGTCCACATCATCCTCGCCTTCTGGGCCTTGTGGCACTGCCGTCACGCGCGCCCAcc                                                                     <  2:209195/135‑1 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                 
TTCGGAAAGTCGGGCGATTGTCGCAGATTTGCGGATTCTCTGCTGCCTGCACGCTCCCTGTGGAAGCGGCCTTGTGTCGCGATAGGGCCGCAAGGCCGCTTCCACTTTGACCGTGGTGCTTTCAACACGAATGTTTAGCCAAGTGCGATGTAAAGGTTGCCATTCACTGCTTTAATCAACCTTACGCCGCATGTGAGTATCAAAGGCATGACTCCCGCCCAGCCCATCGCCCCACCGTCGTCCCGCCCTTCGCTACCGAAGGCCCGGAGCCGTGCCGGGGAAAATCCGCTGGTCCACATCATCCTCGCCTTCTGGGCCTTGTGGCACTGCCGTCACGCGCGCCCACCGGACACTTCGCTCACGCCATTGTGACCAAACCCGGCCAGCGTCTGGCCGTTTGACTCAACCGGGCCG  >  AE015451/5167926‑5168339

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ 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.