New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 296265940 (1.030)5 (0.130) 5/126 2.9 11.7% coding (106/324 nt) ppdC putative prepilin peptidase‑dependent protein
?NC_000913 = 2962936 37 (0.980)coding (221/408 nt) ygdB DUF2509 family protein
Rejected: Frequency below/above cutoff threshold.

GTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGC  >  NC_000913/2962592‑2962726
                                                                   |                                                                   
gTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGT                                                                     <  2:1453057‑M1/68‑1 (MQ=255)
     ctgcCAGCCATGCCGCCAGAGCTGTTGGCACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCt                                                                <  2:600485/68‑1 (MQ=255)
           gccaTGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAAc                                                          <  1:672166/68‑1 (MQ=255)
                   gccAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAAc                                                  >  2:494394/1‑68 (MQ=255)
                      agagCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTacgctgaaaccgatgcccaggt                                               >  2:432601‑M1/1‑46 (MQ=255)
                        agCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGt                                             <  1:1143331/68‑1 (MQ=255)
                        agCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGt                                             <  1:1326460/68‑1 (MQ=255)
                         gCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTg                                            >  1:1114611/1‑68 (MQ=255)
                          cTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACg                                                 <  1:1349605/62‑1 (MQ=255)
                          cTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACg                                                 >  2:1349600/1‑62 (MQ=255)
                              tGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTacgctgaaaccgatgcccaggtttgtttgc                                       >  2:1419289‑M1/1‑38 (MQ=255)
                                gTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAatca                                     >  2:1030530/1‑68 (MQ=255)
                                gTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAatca                                     <  1:734123/68‑1 (MQ=255)
                                gTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAatca                                     <  1:732282/68‑1 (MQ=255)
                                   cTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAAtcatca                                  >  2:859130/1‑68 (MQ=255)
                                    ttttttCTGCTGGCAAGACAGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGCGACAATCATCAc                                 <  2:1462638/65‑1 (MQ=255)
                                    tGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTTGATACCCGATAACGCCGTGACAATCATCAc                                 <  2:366773/68‑1 (MQ=255)
                                     ggTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCAcc                                <  1:257806/68‑1 (MQ=255)
                                      gTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCa                               <  1:562825/68‑1 (MQ=255)
                                       tttCTGCTGGCAAGACTGTTCATTAATGTacgctgaaaccgatgcccaggtttgtttgcgtttactgg                              >  1:1433144‑M1/1‑29 (MQ=255)
                                       tttCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGt                                             <  2:1327591/53‑1 (MQ=255)
                                       tttCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGt                                             >  1:1327595/1‑53 (MQ=255)
                                       tttCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCAt                              >  2:792035/1‑68 (MQ=255)
                                       tttCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCAt                              >  1:233329/1‑68 (MQ=255)
                                       tttCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCAt                              >  1:98754/1‑68 (MQ=255)
                                       tttCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCAt                              >  2:787132/1‑68 (MQ=255)
                                         tCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCAACATCa                            >  1:1372030/1‑68 (MQ=255)
                                            gctgGCAAGACTGTTCATTAATGTTCGCTGGAAAccc                                                        >  2:1227185/1‑37 (MQ=255)
                                            gctgGCAAGACTGTTCATTAATGTTCGCTGGAAAccc                                                        <  1:1227189/37‑1 (MQ=255)
                                            gctgGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAaca                         >  2:315345/1‑68 (MQ=255)
                                            gctgGCAAGACTGGGCATAAATGTTCGCAGGAAACCCGATAACGCAGTGACAATAATCACCATCAaca                         >  2:204344/1‑68 (MQ=255)
                                             ctgGCAAGACTGTTCATTAATGTacgctgaaaccgatgcccaggtttgtttgcgtttactggcagata                        <  2:1433139‑M1/68‑46 (MQ=255)
                                              tgGCAAGACTGTTCATTAATGTacgctgaaaccgatgcccaggtttgtttgcgtttactggcagataa                       <  1:1419294‑M1/68‑47 (MQ=255)
                                              tgGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACAcc                       <  2:398405/68‑1 (MQ=255)
                                                gCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCAt                     <  1:605720/68‑1 (MQ=255)
                                                   aGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGc                  <  2:10947/68‑1 (MQ=255)
                                                        gctCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACa             <  1:1178308/66‑1 (MQ=255)
                                                          tCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACAtt           >  1:596072/1‑68 (MQ=255)
                                                           cATTAATGTTCGCTGGAGACCCGATAACGCAGTGACAATCATCGCCATCAACACCATCGCCAACATTa          >  2:185276/1‑68 (MQ=255)
                                                               aaTGTTCGCTGGAAGCCCGATAACGCAGTGCCAATCATCACCATCAACACCATCGACAACATTACCta      >  1:721852/1‑67 (MQ=255)
                                                               aaTGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTc      >  2:1408595/1‑68 (MQ=255)
                                                               aaTGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTc      >  1:362709/1‑68 (MQ=255)
                                                               aaTGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTc      >  2:614188/1‑68 (MQ=255)
                                                                 tGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACGCCATCGCCAACATTACCTCCg    >  1:919725/1‑68 (MQ=255)
                                                                 tGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCg    >  2:744236/1‑68 (MQ=255)
                                                                 tGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACAGTCGCCAACATTACATCCg    >  1:1285256/1‑68 (MQ=255)
                                                                   tacgctgaaaccgatgcccaggtttgtttgcgtttactggcagataatgaagccttattgattgccgg  <  1:781604‑M1/68‑68 (MQ=255)
                                                                   |                                                                   
GTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGC  >  NC_000913/2962592‑2962726

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.