New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 361460417 (0.850)5 (0.280) 3/118 1.8 23.7% coding (938/1575 nt) nikA nickel/heme ABC transporter periplasmic binding protein
?NC_000913 = 3615035 17 (0.950)coding (1369/1575 nt) nikA nickel/heme ABC transporter periplasmic binding protein

CGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCG  >  NC_000913/3614973‑3615102
                                                              |                                                                   
cGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGAtcggc                                                                >  2:526104‑M2/1‑63 (MQ=255)
      acGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTg                                                          >  2:680011/1‑68 (MQ=255)
        gCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCt                                                        >  2:445467/1‑68 (MQ=255)
                    gCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGa                                            >  1:285782/1‑68 (MQ=255)
                    gCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGa                                            >  1:648956/1‑68 (MQ=255)
                           aaGGATTAGCCGACAAACCGCTGATTGATAAAGAGAtcggcgacctgctgggtgccatacaacgcgtt                                     >  1:378879‑M2/1‑36 (MQ=255)
                              gATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGAc                                           >  1:472391/1‑59 (MQ=255)
                              gATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGAc                                           <  2:472392/59‑1 (MQ=255)
                                            ccGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCAAGACGAAACGCAACGTCAGGCGCTGTa                    <  2:622837/68‑1 (MQ=255)
                                              gCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGt                     <  1:506795/65‑1 (MQ=255)
                                              gCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGt                     >  2:506796/1‑65 (MQ=255)
                                                 gattgatAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATcgcg               >  2:354186/1‑68 (MQ=255)
                                                   ttgatAAAGAGAtcggcgacctgctgggtgccatacaacgcgttatcaatcagcgattttttgtttac             <  1:457339‑M2/68‑57 (MQ=255)
                                                   ttgatAAAGAGAtcggcgacctgctgggtgccatacaacgcgttatcaatcagcgattt                      >  2:358903‑M2/1‑12 (MQ=255)
                                                   ttgatAAAGAGAtcggcgacctgctgggtgccatacaacgcgttatcaatcagcgattt                      <  1:358902‑M2/59‑48 (MQ=255)
                                                   ttgatAAAGAGAtcggcgacctgctgggtgccatacaacgcgttatcaatcagagattttttgtttac             <  2:286623‑M2/68‑57 (MQ=255)
                                                   ttgatAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAGCGCAACGTCAGGCGCTGGATCGGGAc             >  1:308575/1‑68 (MQ=255)
                                                     gatAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACAt           <  1:220288/68‑1 (MQ=255)
                                                     gatAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACAt           >  2:244339/1‑68 (MQ=255)
                                                         aaGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTg       >  1:750798/1‑68 (MQ=255)
                                                          agagATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTgc      >  2:605814/1‑67 (MQ=255)
                                                             gaTCGGCGAAGTGCTGGCGCCCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGAccc   >  1:47094/1‑68 (MQ=255)
                                                              aTCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCg  <  2:259402/68‑1 (MQ=255)
                                                              |                                                                   
CGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCG  >  NC_000913/3614973‑3615102

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