New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3538165 =28 (0.910)4 (0.150) 4/118 2.3 16.7% coding (781/2322 nt) yhgF putative transcriptional accessory protein
?NC_000913 3538246 = 15 (0.550)coding (862/2322 nt) yhgF putative transcriptional accessory protein
Rejected: Frequency below/above cutoff threshold.

GCTGAATGCCGATCCACAGTTCGATGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTGCGCCTGAACAATGCCCCGGCGGATAGCTGGCGCAAAGGCGTGGTGAGCTGGACGTG  >  NC_000913/3538098‑3538232
                                                                   |                                                                   
cgttgattgcttcatcttccgcacgttcgcgcacggtgcccatcagttcggtttccagatgcagcagc                                                                     >  2:771443‑M1/68‑68 (MQ=255)
  tGAATGCCGATCCACAGTTCGATGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCAc                                                                   <  1:1141490/68‑1 (MQ=255)
  tGAATGCCGATCCACAGTTCGATGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCAc                                                                   <  2:516104/68‑1 (MQ=255)
  tGAATGCCGATCCACAGTTCGATGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCAc                                                                   <  1:390019/68‑1 (MQ=255)
             ccACAGTTCGATGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTgcg                                                        <  2:353449/68‑1 (MQ=255)
             ccACAGTTCGATGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTgcg                                                        <  2:756518/68‑1 (MQ=255)
             ccACAGTTCGATGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTgcg                                                        <  2:336201/68‑1 (MQ=255)
             ccACAGTTCGATGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTgcg                                                        <  1:595960/68‑1 (MQ=255)
             ccACAGTTCGATGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTgcg                                                        <  2:580716/68‑1 (MQ=255)
             ccACAGTTCGATGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTgcg                                                        <  2:640050/68‑1 (MQ=255)
               acaGTTCGATGAGCCGCCCAAAGAGAGCTAATGCGAGCAAATCATCATGGATCACCTTGGCCTGAGcc                                                      >  2:219127/1‑68 (MQ=255)
                 aGTTCGATGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTGCGCCTg                                                    >  1:353646/1‑68 (MQ=255)
                    tCGATGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTGCGCCTGAAc                                                 >  1:1059532/1‑68 (MQ=255)
                       aTGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTGCGCCTGAACAAt                                              <  2:623290/68‑1 (MQ=255)
                               cccAAAGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTGCGCCTGAACAATGCCCcggc                                      <  2:191395/68‑1 (MQ=255)
                                ccAAAGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTGCGCCTGAACAATGCCCcggcg                                     <  1:799838/68‑1 (MQ=255)
                                  aaaGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTGCGCCTGAACAATGCCCCGGCGGa                                   <  1:405475/68‑1 (MQ=255)
                                      agagATAATGAGAGCAAATCATCAAGGAACACCATGGCCTGCGCCTGACAACTGACCCGGCGGATAGc                               >  1:199332/1‑68 (MQ=255)
                                        catcagttcggtttccagatgcatcagCACCTTGGCCTGCGCCTGAACAATGCCCCGGCGGATAGCTg                             <  2:953370‑M1/41‑1 (MQ=255)
                                          tcagttcggtttccagatgcatcagCACCTTGGCCTGCGCCTGa                                                   >  1:506962‑M1/26‑44 (MQ=255)
                                          tcagttcggtttccagatgcatcagCACCTTGGCCTGCGCCTGa                                                   <  2:506962‑M1/19‑1 (MQ=255)
                                          cTATTGCGAGCAAATCATCATGGATCACCTTGGCCTGCGCCTGAACAATGCCCCGGCGGATAGCTGgc                           <  2:453893/68‑1 (MQ=255)
                                              tGCGAGCAAATCATCATGGATCACCTTGGCCTGCGCCTGAACAATGCCCCGGCGGATAGCTGGCGCaa                       >  1:753234/1‑68 (MQ=255)
                                                 gAGCAAATCATCATGGATCACCTTGGCCTGCGCCTGAACAATGCCCCGGCGGATAGCTGGCGCAAAgg                    <  1:671294/68‑1 (MQ=255)
                                                 gAGCAAATCATCATGGATCACCTTGGCCTGCGCCTGAACAATGCCACGGCGGATAGCTGGCGCAAAgg                    >  1:737617/1‑68 (MQ=255)
                                                   gCAAATCATCATGGATCACCTTGGCCTGCGCCTGAAAAATGCCCCGGCGGATAGCTGGCGCAAAGGCg                  >  2:980174/1‑68 (MQ=255)
                                                        tcatcaTGGATCGCCTTGGCCTGCGCCTGAACAATGCCCCGGCGGATAGCTGGCGCAGAGGCgtgggg             >  2:733457/1‑66 (MQ=255)
                                                        tcatcaTGGATCACCTTGGCCTGCGCCTGAACAATGCCCCGGCGGATAGCTGGCGCAAAGGCgtggtg             <  1:135022/68‑1 (MQ=255)
                                                           tcatGGATCACCTTGGCCTGCGCCTGAACAATGCCCCGGCGGATAGCTGGCGCAAAGGCGTGGTGAGc          >  1:216887/1‑68 (MQ=255)
                                                           tcatGGATCACCTTGGCCTGCGCCTGAACAATGCCCCGGCGGATAGCTGGCGCAAAGGCGTGGTGAGc          >  1:148324/1‑68 (MQ=255)
                                                             catcagCACCTTGGCCTGCGCCTGAACAATGCCCCGGCGGATAGCTGGCGCAAAGGCGTGGTGAGCTg        >  1:953373‑M1/7‑68 (MQ=255)
                                                               ggATCACCTTGGCCTGCGCCTGAACAATGCCCCGGCGGATAGCTGGCGCAAAGGCGTGGTGAGCTGGa      >  1:136695/1‑68 (MQ=255)
                                                                  tCACCTTGGCCTGCGCCTGAACAATGCCCCGGCGGATAGCTGGCGCAAAGGCGTGGTGAGCTGGACGt   >  2:1044600/1‑68 (MQ=255)
                                                                   CACCTTGGCCTGCGCCTGAACAATTCCACGGCGGATTGCTGGCGCTAAGTCGTGGTGAGCTGGACgtt  >  2:268833‑M1/1‑67 (MQ=255)
                                                                   |                                                                   
GCTGAATGCCGATCCACAGTTCGATGAGCCGCCCAAAGAGAGCTATTGCGAGCAAATCATCATGGATCACCTTGGCCTGCGCCTGAACAATGCCCCGGCGGATAGCTGGCGCAAAGGCGTGGTGAGCTGGACGTG  >  NC_000913/3538098‑3538232

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