New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4164667 =34 (1.110)7 (0.250) 3/122 2.7 21.7% coding (1029/1845 nt) btuB vitamin B12/cobalamin outer membrane transporter
?NC_000913 4164716 = 19 (0.670)coding (1078/1845 nt) btuB vitamin B12/cobalamin outer membrane transporter

AAACAATGTCATCGTTGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCT  >  NC_000913/4164601‑4164732
                                                                  |                                                                 
aaaCTATGTCATCGTTGGGCACGGTAGTATTGGTGCGGGTGGCGACTGGCAGAAACAGACTACGACGc                                                                  <  2:377639/68‑1 (MQ=255)
  cgtcactgcgtgctgcgccttcaaaggtaaaatcgccgacttgttgcagcccggtcagatagatGCCg                                                                <  1:1127197‑M1/4‑1 (MQ=255)
  cgtcactgcgtgctgcgccttcaaaggtaaaatcgccgacttgttgcagcccggtcagatagatGCCg                                                                <  2:575338‑M1/4‑1 (MQ=255)
   gtcactgcgtgctgcgccttcaaaggtaaaatcgccgacttgttgcagcccggtcagatagatGCCgg                                                               >  1:575338‑M1/64‑68 (MQ=255)
   cAATGTCATCGTTGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCgg                                                               >  1:1110444/1‑68 (MQ=255)
   cAATGTCATCGTTGGGCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCgg                                                               <  1:727721/68‑1 (MQ=255)
     aTGTCATCGTTGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGt                                                             <  2:38803/68‑1 (MQ=255)
       gTCATCGTTGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTAc                                                           >  2:406658/1‑68 (MQ=255)
        tCATCGTTGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGACGGGTACa                                                          >  2:702202/1‑68 (MQ=255)
        gcATCGTTGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACa                                                          <  2:163138/67‑1 (MQ=255)
         cATCGTTGGTCACGTTATTATTGGTGAGGTTGGCGACTTGCAGAAACAGACTACTACTCCGGGTACat                                                         <  2:604581/68‑2 (MQ=255)
         cATCGTTGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAg                                                         <  1:249840/68‑1 (MQ=255)
         cATCGTTGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAg                                                         <  2:285565/68‑1 (MQ=255)
           tCGTTGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGt                                                       <  2:532173/68‑1 (MQ=255)
           tCGTTGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGt                                                       <  1:834230/68‑1 (MQ=255)
           tCGTTGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGt                                                       <  2:2070/68‑1 (MQ=255)
            cGTTGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGtt                                                      >  1:146790/1‑68 (MQ=255)
              ttGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGTTAt                                                    <  1:948157/68‑1 (MQ=255)
              ctgcgccttcaaaggtaaaatcgccgacttgttgcagcccggtcagatagctGCCGGGTACAGGTTAt                                                    <  1:93369‑M1/16‑1 (MQ=255)
              ctgcgccttcaaaggtaaaatcgccgacttgttgcagcccggtcagatagatGCCGGGTACAGGTTAt                                                    <  1:200678‑M1/16‑1 (MQ=255)
              ctgcgccttcaaaggtaaaatcgccgacttgttgcagcccggtcagatagatGCCGGGTACAGGTTAt                                                    <  1:94822‑M1/16‑1 (MQ=255)
                ggTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGTGTACAGGTTATGt                                                  <  1:541160/68‑1 (MQ=255)
                ggTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGTTATGt                                                  >  1:550036/1‑68 (MQ=255)
                 gTCACGGTAGTATTGGTGCGGGCGGCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGTTATGtt                                                 <  1:295047/68‑1 (MQ=255)
                          gtaTTTGTGCGGGTGTCGACTGTCAGAAACAGACTACGACGCCGTGTACAGGTTATGTTGAggatgga                                        <  1:215834/68‑1 (MQ=255)
                          gtaTTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGTTATGTTGAggatgga                                        <  2:1110439/68‑1 (MQ=255)
                          gtaTTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGTTATGTTGAggatgga                                        <  1:679853/68‑1 (MQ=255)
                          gtaTTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGTTATGTTGAggatgga                                        <  2:417184/68‑1 (MQ=255)
                            aTTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGTTATGTTGAGGATGGata                                      >  2:936239/1‑68 (MQ=255)
                              tGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGTTATGTTGAGGATGGATATg                                    <  2:283380/68‑1 (MQ=255)
                                 tGCGGGTGTCGACTGGCAGAAACAGACTACGCCGCCGGGTACAGGTTATGTTGAGGATGGATATGATc                                 >  1:582180/1‑68 (MQ=255)
                                 tGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGt                                                       >  2:236871/1‑46 (MQ=255)
                                 tGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGt                                                       <  1:236871/46‑1 (MQ=255)
                                 tGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGTTATGTTGAGGATGGATATGATc                                 >  2:849744/1‑68 (MQ=255)
                                    gggTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAAc                              <  1:284127/68‑1 (MQ=255)
                                    gggTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAAc                              <  2:403523/68‑1 (MQ=255)
                                    gccgacttgttgcagcccggtcagatagatGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAAc                              >  2:200678‑M1/31‑68 (MQ=255)
                                    gccgacttgttgcagcccggtcagatagatGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAAc                              >  2:850426‑M1/31‑68 (MQ=255)
                                     ccgacttgttgcagcccggtcagatagatGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACg                             >  2:94822‑M1/30‑68 (MQ=255)
                                     ccgacttgttgcagcccggtcagatagatGCCGGGTACAGGGTATGTTGAGGATGGATATGATCAACg                             >  2:93369‑M1/30‑68 (MQ=255)
                                            aCTGGCAGAAACAGACTACGACGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACGTAATAcc                      >  2:334644/1‑68 (MQ=255)
                                                           cTACGACGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGAcc       >  2:952209/1‑68 (MQ=255)
                                                              cgacgCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCggg    <  1:78870/68‑1 (MQ=255)
                                                                acgCCGGGTACAGGTTATGTTGAGGATGGATATGATCAGCGTAATACCGGCATCTATCTGACCGGGCt  >  2:1127192/1‑68 (MQ=255)
                                                                  GCCGGGTACAGGTTATGTTGAGGATGGATATGATCAAc                              <  1:621595‑M1/38‑1 (MQ=255)
                                                                  GCCGGGTACAGGTTATGTTGAGGATGGATATGATCAAc                              >  2:621595‑M1/1‑38 (MQ=255)
                                                                  |                                                                 
AAACAATGTCATCGTTGGTCACGGTAGTATTGGTGCGGGTGTCGACTGGCAGAAACAGACTACGACGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCT  >  NC_000913/4164601‑4164732

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.