New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 35306616 (0.620)5 (0.210) 4/272 3.0 19.2% intergenic (+33/+85) hisF/cbcV imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein
?NC_002947 = 353120 27 (1.120)intergenic (+87/+31) hisF/cbcV imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein
Rejected: Frequency below/above cutoff threshold.

TTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACC  >  NC_002947/352921‑353211
                                                                                                                                                 |                                                                                                                                                 
ttCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCgg                                                                                                                                                 >  7:64570‑M1/1‑146 (MQ=255)
ttCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCgg                                                                                                                                                 <  5:98995‑M1/148‑3 (MQ=255)
        ctggccgccagcaCCGCGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCgggcacgccc                                                                                                                                         >  6:3566‑M1/7‑138 (MQ=255)
                                             aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCgg                                                                                                     <  8:55010/147‑1 (MQ=255)
                                                ttGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                                 >  3:89183/1‑148 (MQ=255)
                                                       gTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCgc                                                                                           <  7:105706‑M2/147‑3 (MQ=255)
                                                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCCCCCGCGCCCACCCGGTTGGTGAGCTCTGTGTGGGCGCGGGCGTGGCCGCGCAGGGCCGCGCACCGGcc                                                                        >  1:97756/1‑148 (MQ=255)
                                                                             gCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCCCAGGGTTGCTGCGATGGCTGGGGGAGCGGGCGGGGCCGCGCAGGGCCGCGAAGCGGccccca                                                                    >  1:79906/1‑146 (MQ=37)
                                                                               cTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCACAGGTTTGCGGCGATCTCTGTGGGAGCCGGCGTGCCCGCGACGGGCCCCGAAGCGGCCCCaaat                                                                  >  1:94986/1‑146 (MQ=38)
                                                                                   aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCt                                                               >  8:95340/1‑147 (MQ=255)
                                                                                       cTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCgggcacgcccgctcccacagagatctcagcaaacctgtgggag                                                                                                        >  4:80122‑M1/1‑59 (MQ=255)
                                                                                       cTTCGCGCGGGAGCGTCGATCGCTGATCCATTGCCCGGGGGGGCTTTGCAGCCCTGCGCGGGCACCCCCGCGCCCACAGGTTTCCGGACATCTCTGTCGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAAc                                                          <  3:59900/148‑1 (MQ=37)
                                                                                        ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCgggcacgcccgctcccacagagatctcagcaaacctgtgggagc                                                                                                       <  3:80122‑M1/102‑45 (MQ=255)
                                                                                        ttCGCGCGGGAGCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAAcc                                                         <  7:89371/148‑1 (MQ=255)
                                                                                              cgGGATCGGCGTGCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGtt                                                   <  1:87769/148‑1 (MQ=255)
                                                                                                                 cccATTGCCCGGGGCGGCTTGCCAGCCCTTCGCGGGCCCCCCCGCCCCCACAGGTTTCCTGAGAGGTCTGGGGGCGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTgg                                <  5:105930/147‑1 (MQ=37)
                                                                                                                  ccATTGACCGGGGCTGCTTTGCAGCCCTGCGCGGGCACGCCCGGTCCCACAGGTTTGCTGAGATCTCTGTGGGTGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTgg                                <  2:97756/147‑1 (MQ=255)
                                                                                                                    aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGGTGGGATGTCTGTGGGCGCGGGCGTGCCCGCGCAGGGCCGCGCAGCGGCCCCCCTTGCTTCAACCGGGGGTCTTGCCGCGCAGCGCGGGGt                               >  7:80176/1‑146 (MQ=25)
                                                                                                                     ttGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                            <  7:95340/148‑1 (MQ=255)
                                                                                                                          cGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                       <  4:89183/148‑1 (MQ=255)
                                                                                                                            gggCTGCTTTGCAGCCCTGCGCGGGCACGCCCGCGCCCACCGGTTTGGTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGt                     <  5:78824/148‑1 (MQ=38)
                                                                                                                                      gCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa             <  2:15892/146‑1 (MQ=255)
                                                                                                                                                gcgggcacgcccgctcccacagagatctcagcaaacctgtgggagcgggcgtgcccgcgaagggctgcaaagcagccccggtcaatggatcagcgaacgacgatcccgcgcgaagccatgtaggccttggcctccggcactgtgtac  <  7:105706‑M1/147‑146 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACC  >  NC_002947/352921‑353211

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