New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 24300917 (0.660)4 (0.170) 4/118 1.9 16.6% coding (295/2445 nt) fadE acyl coenzyme A dehydrogenase
?NC_000913 = 243192 25 (1.080)coding (112/2445 nt) fadE acyl coenzyme A dehydrogenase
Rejected: Frequency below/above cutoff threshold.

TCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACC  >  NC_000913/242943‑243070
                                                                  |                                                             
tCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTcc                                                              <  1:883696/68‑1 (MQ=255)
gcAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTAGCCCTGGAACAAGTCGCCCTcc                                                              <  1:356004/67‑1 (MQ=255)
 cAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTccc                                                             <  2:953966/68‑1 (MQ=255)
   ggCGCGTCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGTCTTGCCCTGGAACAAGTCGCCCTCccac                                                           <  2:444056/68‑1 (MQ=255)
      gcgGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAg                                                        >  1:319294/1‑68 (MQ=255)
           tGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAggtggt                                                   <  2:816379/68‑1 (MQ=255)
           cgCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAggtggt                                                   <  2:899062/67‑1 (MQ=255)
              ggATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAGCGAGTCGCCCTCCCACCAGGTGGTGcc                                                >  2:356540/1‑68 (MQ=255)
              ggATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCggcgttgctggtctg                                                >  2:145110‑M1/1‑53 (MQ=255)
                     tATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCa                                         >  2:413468/1‑68 (MQ=255)
                      aTGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCggcgttgctggtctgtggtcggc                                        >  2:632309‑M1/1‑45 (MQ=255)
                      aTGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCaa                                        <  2:448989/68‑1 (MQ=255)
                          agatttttCCAGTCCGGCTTGCCCTGGAACAAGACGCCCTCCCACCAGGTGGTGCCAGCATCAATCGc                                    >  1:652039/4‑68 (MQ=255)
                            cTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGctt                                  <  1:607302/68‑1 (MQ=255)
                            cTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGctt                                  <  2:765186/68‑1 (MQ=255)
                                     gTCCGGCTTGCCCTGGAACAAGTCGCCCTCggcgttgctggtctgtggtcggcgtgggtactggtgcc                         >  1:726771‑M1/1‑30 (MQ=255)
                                        cGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTgc                      <  1:302435/68‑1 (MQ=255)
                                            ttGCCCTGGAACAAGTCGCCCTCggcgttgctggtctgtggtcggcgtgggtactggtgcctctggcc                  >  2:631327‑M1/1‑23 (MQ=255)
                                            ttGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGa                  <  2:771272/68‑1 (MQ=255)
                                            ttGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCAGCATCAATCGCTTCTTTCTCAGTGCGCGa                  >  2:457997/1‑68 (MQ=255)
                                                      aCAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCATCTTTCTCAGTGCGCGACATCggcggc        >  2:874088/1‑68 (MQ=255)
                                                            cGCCCTCggcgttgctggtctgtggtcggcgtgggtactggtgcctctggccattatcctcgtgccat  >  1:953841‑M1/1‑7 (MQ=255)
                                                            cGCCCTCggcgttgctggtcggtggtcggcgtgggtactggtgcctctggccattatcctcgtgccat  <  2:902164‑M1/68‑62 (MQ=255)
                                                                  |                                                             
TCAGGCGCGGCTGCGGATAGTTATGCAGCTTTTTCCAGTCCGGCTTGCCCTGGAACAAGTCGCCCTCCCACCAGGTGGTGCCCGCATCAATCGCTTCTTTCTCAGTGCGCGACATCGGCGGCATCACC  >  NC_000913/242943‑243070

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