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

ATTACGCGGTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGACACCCTGTTTGCCCCTTCTGTGCCCTACGCCAACCTCGGCCTGAAACCGAGCCAGTACGATCCGCAA  >  NC_000913/3614541‑3614671
                                                               |                                                                   
aTTACGCGGTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGACAcc                                                                 >  2:509709/1‑68 (MQ=255)
aTTACGCGGTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGACAcc                                                                 >  1:47419/1‑68 (MQ=255)
  tACGCGGTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGACACCCt                                                               >  1:373948/1‑68 (MQ=255)
  tACGCGGTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGACACCCt                                                               >  1:433568/1‑68 (MQ=255)
       ggTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGACACCCTGTTTg                                                          <  2:471150/68‑1 (MQ=255)
        gTAAACAAAAAATCTCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGAtctctttatcaa                                                         <  2:286623‑M1/68‑13 (MQ=255)
        gTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGAtctctttatcaa                                                         <  1:457339‑M1/68‑13 (MQ=255)
               aaaaaTCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGACACCCTGTTTGCCCCTTCt                                                  >  1:341466/1‑68 (MQ=255)
                 aaaTCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGAtctctttatcaa                                                         <  1:358902‑M1/59‑13 (MQ=255)
                 aaaTCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGAtctctttatcaa                                                         >  2:358903‑M1/1‑47 (MQ=255)
                 aaaTCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGACACCCTGTTTGCCCCTTCtgt                                                <  2:351207/68‑1 (MQ=255)
                   aTCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGACACACTGTTTGCCCCTTCTGTGc                                              >  2:703149/1‑68 (MQ=255)
                                aaCGCGTTGTATGGCACCCAGCAGGTCGCCGAtctctttatcaatcagcggtttgtcggctaatcctt                                 >  1:378879‑M1/1‑32 (MQ=255)
                                    cgTTGTATGGCACCCAGCAGGTCGCCGACACCCTGTTTGCCCCTTCTGTGCCCTACGCCAACCTCGGc                             >  1:526103/1‑68 (MQ=255)
                                          aTGGCACCCAGCAGGTCGCCGACACCCTGTTTGCCCCTTCTGTGCCCTACGCCAACCTCGGCCTGaaa                       >  1:687/1‑68 (MQ=255)
                                             gCACCCAGCAGGTCGCCGACACCCTGTTTGCCCCTTCTGTGCCCTACGCCAACCTCGGCCTGAAACCg                    >  1:280659/1‑68 (MQ=255)
                                               aCCCAGCAGGTCGCCGACACCCTGTTTGCCCCTTCTGTGCCCTGCGCCAACCTCGGCCTGAAGCCgcg                  >  1:702237/1‑66 (MQ=255)
                                               aCCCAGCAGGTCGCCGACACCCTGTTTGCCCCTTCTGTGCCCTACGCCAACCTCGGCCTGAAACCGAg                  >  2:296364/1‑68 (MQ=255)
                                               aCCCAGCAGGTCGCCGACACCCTGTTTGCCCCTTCTGTGCCCTACGCCAACCTCGGCCTGAAACCGAg                  <  2:544936/68‑1 (MQ=255)
                                                       ggTCGCCGACACCCTGTTTGCCCCTTCTGTGCCCTACGCCAACCTCGGCCTGAAACCGAGCCAGTACg          >  1:95130/1‑68 (MQ=255)
                                                           gCCGAtctctttatcaatcagcggtttgtcggctaatccttgctgtgcctggaagtcagcgtgtgacg      >  2:526104‑M1/1‑5 (MQ=255)
                                                              gACACACTGTTTGCCCCTTCTGTGCCCTACGCCAACCTCGGCCTGAGACCGAGCCAGTACGATCCGCa   >  1:747480/1‑68 (MQ=255)
                                                               acacacTGTTTGAAACTTCTGTGCCCTACGCCAACCAAGGCCTGAAACCGAGCCAGTGCGATCCGCaa  >  1:749101/1‑68 (MQ=255)
                                                               |                                                                   
ATTACGCGGTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGACACCCTGTTTGCCCCTTCTGTGCCCTACGCCAACCTCGGCCTGAAACCGAGCCAGTACGATCCGCAA  >  NC_000913/3614541‑3614671

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