New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 = 479189768 (0.950)3 (0.040) 3/334 NT 4.0% coding (1943/6525 nt) pvdI non‑ribosomal peptide synthetase (subunit of ferribactin synthase)
?AE015451 = 4791939 80 (1.170)coding (1901/6525 nt) pvdI non‑ribosomal peptide synthetase (subunit of ferribactin synthase)
Rejected: Frequency below/above cutoff threshold.

CGACAGGTAGGCGGCGAAACCGTCGGCATCGAGCACCGCCTCGCGCTGCAGCAGATGTACCGTCTTGCCGGCCCAGAGCGCGCCAAACAGCATGGTATGCCCAAGGTCGGCGGCAGGCGTGGAGACCATGGCGAAGTGTTCCAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/4791659‑4791897
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggtgcccgTGGCGGTTCACTCCGGCAACCTGGCTTATGTGATCTACACCTCGGGCACCACCGGTATGCCCAAGGGGGTCGCGGTGTCTCATGGGGCGCTGGCCAATTACCTTGAAGGCTTGTCCGAGCGGATTGCCCT  <  AE015451/4791939‑4791810
                                                                                                                                                                                                                                                                                                                                                                                 
CGACAGGTAGGCGGCGAAACCGTCGGCATCGAGCACCGCCTCGCGCTGCAGCAGATGTACCGTCTTGCCGGCCCAGAGCGCGCCAAACAGCATGGTATGCCCAAGGTCGGCGGCAGGCGTGGAGACCATGGCGAAGTGTTCCAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGC                                                                                                                                       <  1:39709/236‑1
CGACAGGGAGGCGGCGAAACCGTCGGCATCGAGCCCCGCCTCGCGCTGCAGCAGATGTACCGTCTTGCCGGCCCAGAGCGCGCCAAACAGCATGGTATGCCCAAGGTCGGCGGCAGGCGTGGAGACCATGGCGAAGTGTTCCAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGC                                                                                                                                       <  1:147219/236‑1
 GACAGGTAGGCGGCGAAACCGTCGGCATCGAGCACCGCCTCGCGCTGCAGCAGATGTACCGTCTTGCCGGCCCAGAGCGCGCCAAACAGCATGGTATGCCCAAGGTCGGCGGCAGGCGTGGAGACCATGGCGAAGTGTTCCAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTG                                                                                                                                        <  2:1200622/234‑1
 GACAGGTAGGCGGCGAAACCGTCGGCATCGAGCACCGCCTCGCGCTGCAGCAGATGTACCGTCTTGCCGGCCCAGAGCGCGCCAAACAGCATGGTATGCCCAAGGTCGGCGGCAGGCGTGGAGACCATGGCGAAGTGTTCCAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTG                                                                                                                                        <  2:586593/234‑1
   CAGGTAGGCGGCGAAACCGTCGGCATCGAGCACCGCCTCGCGCTGCAGCAGATGTACCGTCTTGCCGGCCCAGAGCGCGCCAAACAGCATGGTATGCCCAAGGTCGGCGGCAGGCGTGGAGACCATGGCGAAGTGTTCCAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCC                                                                                                                                     <  1:298913/235‑1
                                                                                                   CCCAAGGTCGGCGGCAGGCGTGGAGACCATGGCGAAGTGTTCCAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCGTGGCGGTTCACTCCGGCAACCTGGCTTATGTGATCTACACCTC                                                                                         <  2:1136051/183‑1
                                                                                                   CCCAAGGTCGGCGGCAGGCGTGGAGACCATGGCGAAGTGTTCCAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATGGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCC                                                                                                                                     <  1:704293/139‑1
                                                                                                                                                       AGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCGTGGCGGTTCACTCCGGCAACCTGGCTTATGTGATCTACACCTCGGGCA                                                                                    >  1:671481/1‑136
                                                                                                                                                       AGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCGTGGCGGTTCACTCCGGCAACCTGGCTTATGTGATCTACACCTCGGGCA                                                                                    <  2:671476/136‑1
                                                                                                                                                                                                                                          gcagGTGGCGGTTCACTCCGGCAACCTGGCTTATGTGATCTACACCTCGGGCACCACCGGTATGCCCAAGGGGGTCGCGGTGTCTCATGGGGCGCTGGCCAATTACCTTGAAGGCTTGTCCGAGCGGATTGCCCT  <  1:800850/131‑1
                                                                                                                                                                                                                                           CCCGTGGCGGTTCTCTCCGGCA                                                                                                                  >  2:221317/1‑22
                                                                                                                                                                                                                                                                                                                                                                                 
CGACAGGTAGGCGGCGAAACCGTCGGCATCGAGCACCGCCTCGCGCTGCAGCAGATGTACCGTCTTGCCGGCCCAGAGCGCGCCAAACAGCATGGTATGCCCAAGGTCGGCGGCAGGCGTGGAGACCATGGCGAAGTGTTCCAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/4791659‑4791897
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggtgcccgTGGCGGTTCACTCCGGCAACCTGGCTTATGTGATCTACACCTCGGGCACCACCGGTATGCCCAAGGGGGTCGCGGTGTCTCATGGGGCGCTGGCCAATTACCTTGAAGGCTTGTCCGAGCGGATTGCCCT  <  AE015451/4791939‑4791810

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ 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.