New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4159113 =31 (1.000)7 (0.250) 7/260 2.5 20.0% intergenic (+17/‑192) PP_3661/PP_3662 membrane protein/hypothetical protein
?NC_002947 4159201 = 28 (0.990)intergenic (+105/‑104) PP_3661/PP_3662 membrane protein/hypothetical protein
Rejected: Frequency below/above cutoff threshold.

GCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTT  >  NC_002947/4158966‑4159251
                                                                                                                                                   |                                                                                                                                          
gCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTTGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTCGTTGGCGGCTGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGgctg                                                                                                                                          <  2:171899/150‑1 (MQ=255)
ctttggggctgacggatgcattgatgacatggcgttgaggcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcctg                                                                                                                                          <  2:50627‑M1/3‑1 (MQ=255)
  ttggggctgacggatgcattgatgacatggcgttgaggcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCt                                                                                                                                        >  1:148055‑M1/146‑150 (MQ=255)
     gggctgacggatgcattgatgacatggcgttgaggcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTcg                                                                                                                                     <  1:371048‑M1/8‑1 (MQ=255)
      ggctgacggatgcattgatgacatggcgttgaggcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTcgc                                                                                                                                    <  1:76041‑M1/9‑1 (MQ=255)
      ggctgacggatgcattgatgacatggcgttgaggcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTcgc                                                                                                                                    >  2:493669‑M1/142‑150 (MQ=255)
       gaCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTcgcg                                                                                                                                   >  2:25841/1‑150 (MQ=255)
        ctgacggatgcattgatgacatggcgttgaggcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCgg                                                                                                                                  >  2:85170‑M1/140‑150 (MQ=255)
          gacggatgcattgatgacatggcgttgaggcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGcgggg                                                                                                                                <  2:701344‑M1/13‑2 (MQ=255)
          ccGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGt                                                                                                                                >  2:252368/1‑150 (MQ=255)
                           gCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAgg                                                                                                               >  2:544359/1‑150 (MQ=255)
                            ccAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggt                                                                                                              <  1:533886/150‑1 (MQ=255)
                                  tGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGtgc                                                                                                        >  2:104318/1‑150 (MQ=255)
                                   gggCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGtgct                                                                                                       >  1:71774/1‑150 (MQ=255)
                                        aCGGGCGCTGGCCACTGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCt                                                                                                  <  1:217267/150‑1 (MQ=255)
                                           ctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGaa                                                                                               <  2:704573‑M1/46‑1 (MQ=255)
                                            gcgcTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAg                                                                                              >  1:223685/1‑150 (MQ=255)
                                               cttgCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACg                                                                                           <  2:45580/147‑1 (MQ=255)
                                                 ggCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGc                                                                                         <  1:290893/150‑1 (MQ=255)
                                                            gggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGt                                                                              <  1:647320‑M1/63‑1 (MQ=255)
                                                              gtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAg                                                                            >  2:76041‑M1/86‑150 (MQ=255)
                                                                ccTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                          <  2:645087/150‑1 (MQ=255)
                                                                 cTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCCCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAg                                                                         >  1:614859/1‑150 (MQ=255)
                                                                    aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCgg                                                                      <  1:717381/150‑1 (MQ=255)
                                                                     atccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                     >  2:647320‑M1/79‑150 (MQ=255)
                                                                        tgtgGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTc                                                                  >  1:701344/1‑150 (MQ=255)
                                                                         ttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagtcggCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGGTCa                                                                 <  2:512120‑M1/76‑1 (MQ=255)
                                                                           ggCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAcc                                                               <  1:423870/150‑1 (MQ=255)
                                                                                    gCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgca                                                      <  1:104318/150‑1 (MQ=255)
                                                                                       gggggCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgcaggc                                                   <  1:85664/150‑1 (MQ=255)
                                                                                            cgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTc                                              >  1:704573‑M1/56‑150 (MQ=255)
                                                                                            cTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTc                                              <  1:476096/150‑1 (MQ=255)
                                                                                                 gATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGtct                                         >  1:99604/1‑150 (MQ=255)
                                                                                                 gATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGtct                                         >  1:99596/1‑150 (MQ=255)
                                                                                                 gATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGGTCAAACGATAATCAGGCAGCTCGGtcc                                         >  1:36840/1‑149 (MQ=38)
                                                                                                    tACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCaa                                      <  2:71774/150‑1 (MQ=255)
                                                                                                          ggcaacgaggggggctgagaaaaaccttgataccgagtcggCTGCTTCGCGGGTAAACCCGCTCCTAGCTGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGtttt                                >  1:512120‑M1/42‑150 (MQ=255)
                                                                                                             ggTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTca                             <  1:25841/150‑1 (MQ=255)
                                                                                                                      gCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAGGCAGGCAGCTCGGTCTCaa                                      <  2:460691/132‑1 (MQ=255)
                                                                                                                      gCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAGGCAGGCAGCTCGGTCTCaa                                      >  1:460691/1‑132 (MQ=255)
                                                                                                                        tGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTc                  <  1:252368/150‑1 (MQ=255)
                                                                                                                                    tGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCcac      >  2:371048/1‑150 (MQ=255)
                                                                                                                                        gaCCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGccccc                     <  1:205591/131‑1 (MQ=255)
                                                                                                                                        gaCCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGccccc                     >  2:205591/1‑131 (MQ=255)
                                                                                                                                        gaCCGAGTCGGCTGCTTCGCGGGTAAAAATGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACAttt  >  1:18964/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
GCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTT  >  NC_002947/4158966‑4159251

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