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.

GTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAG  >  NC_002947/352976‑353267
                                                                                                                                                |                                                                                                                                                   
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)
                                cTTCGCGCGGGAGCGTCGATCGCTGATCCATTGCCCGGGGGGGCTTTGCAGCCCTGCGCGGGCACCCCCGCGCCCACAGGTTTCCGGACATCTCTGTCGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAAc                                                                                                                  <  3:59900/148‑1 (MQ=37)
                                 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)
                                                                                           gggCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCaccacc                                                       <  8:64570/147‑1 (MQ=255)
                                                                                                  cccGCTCCCACAGGTTTGCGGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGt                                                <  6:51007/148‑1 (MQ=255)
                                                                                                     gCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCgcgaagggctgcaaagcagccccggtcaatggatcagcgaacgacgatcccgcgcgaa                                                                                           <  3:80122‑M2/102‑59 (MQ=255)
                                                                                                      cTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCgcgaagggctgcaaagcagccccggtcaatggatcagcgaacgacgatcccgcgcgaag                                                                                          >  4:80122‑M2/1‑43 (MQ=255)
                                                                                                       tCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATcc                                           >  7:13460/1‑148 (MQ=255)
                                                                                                       tCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTGATcc                                           >  6:19083/1‑148 (MQ=255)
                                                                                                         ccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCACGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCt                                          >  1:100216/1‑147 (MQ=255)
                                                                                                          cacaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCt                                          >  6:131108/1‑146 (MQ=255)
                                                                                                                     tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCtt                             >  1:100674/1‑148 (MQ=255)
                                                                                                                     tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCtt                             <  3:22913/148‑1 (MQ=255)
                                                                                                                        gttgtctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttg                            <  2:94986/142‑1 (MQ=255)
                                                                                                                            tctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCg                      >  6:130714/1‑148 (MQ=255)
                                                                                                                               gtgGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGtt                   <  6:28155/148‑1 (MQ=255)
                                                                                                                                   gAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGa                <  5:3566/147‑1 (MQ=255)
                                                                                                                                       gggCGTGCCCgcgaagggctgcaaagcagccccggtcaatggatcagcgaacgacgatcccgcgcgaagccatgtaggccttggcctccggcactgtgtactcgccaaagtggaagatgctggccgccagcgcggtgctggcggccag           >  6:3566‑M2/1‑10 (MQ=255)
                                                                                                                                               ccgcgaagggctgcaaagcagccccggtcaatggatcagcgaacgacgatcccgcgcgaagccatgtaggccttggcctccggcactgtgtactcgccaaagtggaagatgctggccgccagcaccgcgctggcatgcccttccagaa   >  7:64570‑M2/1‑2 (MQ=255)
                                                                                                                                               ccgcgaagggctgcaaagcagccccggtcaatggatcagcgaacgacgatcccgcgcgaagccatgtaggccttggcctccggcactgtgtactcgccaaagtggaagatgctggccgccagcaccgcgctggcatgcccttccagaa   <  5:98995‑M2/148‑147 (MQ=255)
                                                                                                                                                cgcgAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAg  <  8:105706/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAG  >  NC_002947/352976‑353267

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.