New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 516635 =64 (1.100)3 (0.060) 3/110 2.8 4.8% coding (53/780 nt) fetB putative iron ABC exporter membrane subunit FetB
?NC_000913 = 516631 63 (1.260)coding (49/780 nt) fetB putative iron ABC exporter membrane subunit FetB
Rejected: Frequency below/above cutoff threshold.

CGCCCGACGCTCCAGAGAATATCTTTCTCCAGCGCCAGTTTTTCTTTATGGCTAATTAAGATTGCCACCACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/516707‑516635
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaccaccaGCATTAATGCCAGTGCTAATGATTCGTTAGTAATATTATGCGAGTTCATAGCGTGCTTCCTG  <  NC_000913/516631‑516570
                                                                                                                                       
CGCCCGACGCTCCAGAGAATATCTTTCTCCAGCGCCAGTTTTTCTTTATGGCTAATTAAGATTGC                                                                        <  1:3230297/65‑1
   CCGACGCTCCAGAGAATATCTTTCTCCAGCGCCAGTTTTTCTTTATGGCTAATTAAGATTGCCAC                                                                     <  1:2081056/65‑1
    CGACGCTCCAGAGAATATCTTTCTCCAGCGCCAGTTTTTCTTTATGGCTAATTAAGATTGCCACC                                                                    <  1:171865/65‑1
         CTCCAGAGAATATCTTTCTCCAGCGCCAGTTTTTCTTTATGGCTAATTAAGATTGCCACCACCAG                                                               >  1:985263/1‑65
                         TCTCCAGCGCCAGTTTTTCTTTATGGCTAATTAAGATTGCCACCACCAGCATTAATGCCAGTGCT                                               <  1:4269626/65‑1
                                                    TAATTAAGATTGCCACCACCAGCATTAATGCCAGTGCTAATGATTCGTTAGTAATATTATGCGAG                    <  1:1766342/65‑1
                                                                CCACCACCAGCATTAATGCCAGTGCTAATGATTCGTTAGTAATATTATGCGAGTTCATAGCGTGC        <  1:4319628/65‑1
                                                                 CACCACCAGCATTAATGCCAGTGCTAATGATTCGTTAGTAATATTATGCGAGTTCATAGCGTGCT       <  1:4066470/65‑1
                                                                  ACCACCAGCATTAATGCCAGTGCTAATGATTCGTTAGTAATATTATGCGAGTTCATAGCGTGCTT      <  1:2245544/65‑1
                                                                  ACCACCAGCATTAATGCCAGTGCTAATGATTCGTTAGTAATATTATGCGAGTTCATAGCGTGCTT      <  1:2583882/65‑1
                                                                  ACCACCAGCATTAATGCCAGTGCTAATGATTCGTTAGTAATATTATGCGAGTTCATAGCGTGCTT      >  1:70841/1‑65
                                                                     ACCAGCATTAATGCCAGTGCTAATGATTCGTTAGTAATATTATGCGAGTTCATAGCGTGCTTCCT   <  1:1307335/65‑1
                                                                     ACCAGCATTAATGCCAGTGCTAATGATTCGTTAGTAATATTATGCGAGTTCATAGCGTGCTTCCT   <  1:226457/65‑1
                                                                      CCAGCATTAATGCCAGTGCTAATGATTCGTTAGTAATATTATGCGAGTTCATAGCGTGCTTCCTG  >  1:3963705/1‑65
                                                                                                                                       
CGCCCGACGCTCCAGAGAATATCTTTCTCCAGCGCCAGTTTTTCTTTATGGCTAATTAAGATTGCCACCACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/516707‑516635
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaccaccaGCATTAATGCCAGTGCTAATGATTCGTTAGTAATATTATGCGAGTTCATAGCGTGCTTCCTG  <  NC_000913/516631‑516570

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