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.

ACAATGCCCCGGCGGATAGCTGGCGCAAAGGCGTGGTGAGCTGGACGTGGCGCATCAAGGTGCTGATGCATCTGGAAACCGAACTGATGGGCACCGTGCGCGAACGTGCGGAAGATGAAGCAATCAACGT  >  NC_000913/3538184‑3538313
                                                              |                                                                   
cagctcaccacgcctttgcgccagctatccgccggggcattgttcaggcgcaggccaaggtgCTGATg                                                                >  1:953373‑M2/63‑68 (MQ=255)
aCAATGCCCCGGCGGATAGCTGGCGCAAAGGCGTGGTGAGCTGGACGTGGCGCGTCAAGGTGCTGATg                                                                >  2:1132559/1‑68 (MQ=255)
aCAATGCCCCGGCGGATAGCTGGCGCAAAGGCGTGGTGAGCTGGACGTGGCGCATCAAGGTGCTGagg                                                                >  2:584344/1‑66 (MQ=255)
aCAATGCCCCGGCGGATAGCTGGCGCAAAGGCGTGGTGAGCTGGACGTGGCGCATCAAGGTGCTGATg                                                                >  2:877711/1‑68 (MQ=255)
    tGCCCCGGCGGATAGCTGGCGCAAAGGCGTGGTGAGCTGGACGTGGCGCATCAAGGTGCTGATGCATc                                                            <  1:388843/68‑1 (MQ=255)
       cccGGCGGATAGCTGGCGCAAAGGCGTGGTGAGCTGGACGTGGCGCATCAAGGTGCTGATGCATCTgg                                                         >  2:340343/1‑68 (MQ=255)
            cggATAGCTGGCGCAAAGGCGTGGTGAGCTGGACGTGGCGCATCAAGGTGCTGATGCATCTGGAAAcc                                                    <  2:91072/68‑1 (MQ=255)
              gATTGCTGGCGCAAAGGCGGGGTGAGCTGGACGTGGCGCATCAAGGTGCTGATGCATCTGGAAACCGa                                                  <  1:106021/68‑1 (MQ=255)
                   cTGGCGCAAAGGCGTGGTGAGCTGGACGTGGCGCATCAAGGTGCTGATGCATCTGGAAACCGAACTGa                                             <  2:667140/68‑1 (MQ=255)
                   cTGGCGCAAAGGCGGGGTGAGCTGGACGTGGCGCATCAAGGTGCTGATGCAGCTGGAAACCGAACTGa                                             <  2:1059527/68‑1 (MQ=255)
                     cagctatccgccggggcattgttcaggcgcaggccaaggtgCTGATGCATCTGGAAACCGAACTGATg                                           <  2:953370‑M2/27‑1 (MQ=255)
                             ggCGTGGTGAGCTGGACGTGGCGCATCAAGGTGCTGATGCATCTGGAAACCGAACTGATGGGCACCGt                                   <  1:739671/68‑1 (MQ=255)
                                  ggtgAGCTGGACGTGGCGCATCAAGGTGCTGATGCATCTGGAAACCGAACTGATGGGCACCGTGcgcg                              >  1:288360/1‑68 (MQ=255)
                                  ggtgAGCTGGACGTGGCGCATCAAGGTGCTGATGCATCTGGAAACCGAACTGATGGGCACCGTGcgcg                              >  1:1093744/1‑68 (MQ=255)
                                      aGCTGGACGTGGCGCATCAAGGTGCTGATGCATCTGGAAACCGAACTGATGGGCACCGTGCGCGAACg                          <  1:625813/68‑1 (MQ=255)
                                           tcaggcgcaggccaaggtgCTGATGCATCTGGAAACCGAACTGa                                             <  2:506962‑M2/25‑1 (MQ=255)
                                           tcaggcgcaggccaaggtgCTGATGCATCTGGAAACCGAACTGa                                             >  1:506962‑M2/20‑44 (MQ=255)
                                              gcggCGCATCAAGGTGCTGATGCATCTGGAAACCGAACTGATGGGCACCGTGCGCGAACGTGCGGAAg                  <  2:753232/66‑1 (MQ=255)
                                              gTGGCGCATCAAGGTGCTGATGCATCTGGAAACCGAACTGATGGGCACCGTGCGCGAACGTGCGGAAg                  >  2:804582/1‑68 (MQ=255)
                                                             gCTGCTGCATCTGGAAACCGAACTGATGGGCACCGTGCGCGAACGTGCGGAAGATGAAGCAATCAACg   >  2:771443‑M2/2‑68 (MQ=255)
                                                              cTGATGCATCTGTAAACCGAACTGATGGGCACCGTGCGCGAACGTGCTGAAGATGAAGCAATCTActt  <  1:219127/68‑3 (MQ=255)
                                                              |                                                                   
ACAATGCCCCGGCGGATAGCTGGCGCAAAGGCGTGGTGAGCTGGACGTGGCGCATCAAGGTGCTGATGCATCTGGAAACCGAACTGATGGGCACCGTGCGCGAACGTGCGGAAGATGAAGCAATCAACGT  >  NC_000913/3538184‑3538313

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.