New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1304426 =34 (1.350)3 (0.140) 3/116 2.1 9.6% coding (593/909 nt) oppC oligopeptide ABC transporter permease
?NC_000913 1304482 = 27 (1.220)coding (649/909 nt) oppC oligopeptide ABC transporter permease
Rejected: Frequency below/above cutoff threshold.

GCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGT  >  NC_000913/1304361‑1304493
                                                                 |                                                                   
gattcaaagaggatcatgctgggcaccagaagtgatgcgtagaccaccaccacaccgagtacattCgg                                                                   >  2:376608‑M1/66‑68 (MQ=255)
 ctggATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAggcggc                                                                  <  2:198997/68‑1 (MQ=255)
 ctggATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAggcggc                                                                  >  1:712896/1‑68 (MQ=255)
     aaagaggatcttgctgggcaccagaagtgacgcgtagaccaccaccacaccgagtacattCGGCACaa                                                              >  1:649922‑M1/61‑68 (MQ=255)
     aaagaggatcatgctgggcaccagaagtgatgcgtagaccaccaccacaccgagtacattCGGCACaa                                                              >  1:136689‑M1/61‑68 (MQ=255)
     aaagaggatcatgctgggcaccagaagtgatgcgtagaccaccaccacaccgagtacattCGGCACaa                                                              >  1:24746‑M1/61‑68 (MQ=255)
     aaagaggatcatgctgggcaccagaagtgatgcgtagaccaccaccacaccgagtacattCGGCACaa                                                              >  1:316228‑M1/61‑68 (MQ=255)
           cTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGc                                                        >  2:243014/1‑68 (MQ=255)
           cTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGc                                                        <  2:87175/68‑1 (MQ=255)
                aTTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGgtgt                                                   <  2:46027/68‑1 (MQ=255)
                aTTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGgtgt                                                   <  1:507214/68‑1 (MQ=255)
                 ttGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTa                                                  <  2:830735/68‑1 (MQ=255)
                     gCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGa                                              >  1:342880/1‑68 (MQ=255)
                      cGTGGGCAAACCCTGAGTCTGAAGCGCAAAGTGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGAc                                             <  1:908264/68‑1 (MQ=255)
                          ggCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCGCAAGTTGGCGGTGTATCGACATCg                                         >  2:910591/1‑68 (MQ=255)
                          ggCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCg                                         >  2:526694/1‑68 (MQ=255)
                           gCATACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCgg                                        <  2:213839/68‑1 (MQ=255)
                           gCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCgg                                        <  2:469822/68‑1 (MQ=255)
                             aagtgatgcgtagaccaccaccacaccgagtacattCTGCACAAGTTGGCGGTGTATCGACATCGGGc                                      >  2:764384‑M1/37‑68 (MQ=255)
                                tgatgcgtagactaccaccacaccttgtacattCGGCACAAGTTGGCTGTGTATCGACATCGGGCAtt                                   >  2:897729‑M1/34‑68 (MQ=255)
                                cccTGAGTCTGATGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGTCAtt                                   <  2:818506/68‑1 (MQ=255)
                                 ccTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGc                                                        >  1:367794/1‑46 (MQ=255)
                                 ccTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGc                                                        <  2:367796/46‑1 (MQ=255)
                                 ccTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTg                                  >  2:208838/1‑68 (MQ=255)
                                  cTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGt                                 <  1:777564/68‑1 (MQ=255)
                                   tGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGTGCATTGtt                                <  2:750096/68‑1 (MQ=255)
                                       tCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTc                            >  1:127531/1‑68 (MQ=255)
                                              cgcAAAGAGTTTATTGAGGCGGCACAAGCTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACAt                     <  1:241190/68‑1 (MQ=255)
                                                 aaaGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATtgt                  >  2:104305/1‑68 (MQ=255)
                                                 aaaGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATtgt                  <  2:712898/68‑1 (MQ=255)
                                                 aaaGAGTTTATTGAGGCGGCACAAGTTGGCGGGGTATCGACATCGGGCATTGTTATTCGCCACATtgt                  <  1:110945/68‑1 (MQ=255)
                                                 aaaGAGTTTATAGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATtgt                  <  1:919161/68‑1 (MQ=255)
                                                   agagTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGc                >  1:863572/1‑68 (MQ=255)
                                                      gTTTTTTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGa             <  1:870073/68‑1 (MQ=255)
                                                       tttATTGAGTCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTTTTATTCGCCACATTGTGCCGaa            <  2:488688/68‑1 (MQ=255)
                                                       tttATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGaa            <  2:349600/68‑1 (MQ=255)
                                                           ttGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTa        >  2:382743/1‑68 (MQ=255)
                                                           ttGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTa        >  2:103321/1‑68 (MQ=255)
                                                           ttGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTa        >  1:272153/1‑68 (MQ=255)
                                                           ttGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTa        >  1:109194/1‑68 (MQ=255)
                                                            tGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTAc       >  1:17470/1‑68 (MQ=255)
                                                            acattCGGCACAAGTTGGCGGTGTATAGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGtaa       <  2:914166‑M1/63‑2 (MQ=255)
                                                                 CGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGgt  >  2:863575‑M1/1‑68 (MQ=255)
                                                                 |                                                                   
GCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGT  >  NC_000913/1304361‑1304493

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.