New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3951141	32 (0.610)	8 (0.200)	8/190	NT	21.7%	intergenic (‑61/+134)	PP_3486/PP_3487	cytochrome c/hypothetical protein
	?	NC_002947	= 3951201	33 (0.810)					intergenic (‑121/+74)	PP_3486/PP_3487	cytochrome c/hypothetical protein

GGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGG  >  NC_002947/3951053‑3951343                                                                                                                                                       |                                                                                                                                               gggCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTc                                                                                                                                               >  1:1409804‑M2/1‑149 (MQ=255)        ggCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgc                                                                                                                                        >  1:1194508‑M2/1‑142 (MQ=255)           tGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgac                                                                                                                                     >  1:1016247‑M2/1‑139 (MQ=255)               gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgacgcaa                                                                                                                                 <  1:504718‑M2/150‑16 (MQ=255)               gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgacgcaa                                                                                                                                 >  1:666047‑M2/1‑135 (MQ=255)               gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCCGCAGTGGCAGCTAAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgacgcaa                                                                                                                                 <  1:632091‑M2/150‑16 (MQ=255)                     gTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgacgcgaggccgc                                                                                                                           >  2:1086989‑M2/1‑129 (MQ=255)                           atgccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgacgcaaggccgctcccac                                                                                                                     <  1:731534‑M2/150‑28 (MQ=255)                                aaGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAg                                                                                                                >  2:646748/1‑150 (MQ=255)                                 aGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGc                                                                                                               <  1:646748/150‑1 (MQ=255)                                     cgatcttAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGAt                                                                                                           <  1:311469/144‑1 (MQ=255)                                      gCACGTAAGCCATATCGCAAGCCCTGCGGGAACTGCCTTGCGGGTGGAAAGAGCTGCGCAGCCACCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTCCGCG‑CCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATg                                                                                                          >  2:877116/1‑150 (MQ=37)                                            aaGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGACGACAGA‑GGGGAACGATGCCCG‑ACCCTTTCGCAACGCAACGCCGCTCCCACCAAAGCGGATgagatcg                                                                                                    >  2:311469/1‑144 (MQ=37)                                             aGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACaga                                                                                                   <  1:298156/150‑1 (MQ=255)                                               ccATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCCCAGAgg                                                                                                 >  2:891731/1‑150 (MQ=255)                                               ccATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAgg                                                                                                 >  1:481243/1‑150 (MQ=255)                                                   atCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCaa                                                                                             <  2:481243/150‑1 (MQ=255)                                                           cccTGGGGGAGCGGCCTTGGGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCa                                                                                      <  2:250342/150‑1 (MQ=255)                                                             cTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGGGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGGAGGATCAGAGACTTCCTGC‑CGACGCTAGGCCGCGCCCACAAAAGCGGCCGCCACCGAGGCCAATGCCCCAggc                                                                                   >  1:1236779/1‑150 (MQ=19)                                                             cTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGAt                                                                                                           <  1:343413/126‑1 (MQ=255)                                                             cTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGAt                                                                                                           >  1:657317/1‑126 (MQ=255)                                                             cTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGAt                                                                                                           >  2:343413/1‑126 (MQ=255)                                                             cTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGAt                                                                                                           <  2:657317/126‑1 (MQ=255)                                                                tgGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgac                                                                                                                                     >  2:1389874‑M2/1‑86 (MQ=255)                                                                tgGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgac                                                                                                                                     <  1:1389874‑M2/97‑12 (MQ=255)                                                                          cTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACt                                                                      >  2:1270720/1‑150 (MQ=255)                                                                          cTTGCGGCGCGAAATGGCTGCGGAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACt                                                                      >  2:1307257/1‑150 (MQ=255)                                                                           ttGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCCAGCATGAAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTc                                                                     >  1:969453/1‑150 (MQ=255)                                                                           ttGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGATGGCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTc                                                                     >  1:660555/1‑150 (MQ=255)                                                                               gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctca                                                                 <  1:1035590‑M2/150‑80 (MQ=255)                                                                                   cgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgaccg                                                             >  1:607421/1‑150 (MQ=255)                                                                                   cgAAAGGGCTGCGCAGCAGCCCCGGCAGTAGCAGCTCAAGCGCGGAAAGG‑TAGGAGCGCTGCGCG‑CTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccc                                                             <  2:1268714‑M2/150‑84 (MQ=255)                                                                                        gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacct                                                                   <  2:664373‑M2/139‑78 (MQ=255)                                                                                        gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacct                                                                   >  1:664373‑M2/1‑62 (MQ=255)                                                                                            tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgacgcaaggccgctcccacagggcttgcgatatg                                                                                                      >  2:1232103‑M2/1‑58 (MQ=255)                                                                                            tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgacgcaaggccgctcccacagggcttgcgatatg                                                                                                      <  1:1232103‑M2/100‑43 (MQ=255)                                                                                              cgcAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggccctggac                                                  >  2:1101686‑M2/1‑56 (MQ=255)                                                                                                 agcagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt                                               <  2:504718/150‑1 (MQ=255)                                                                                                    agcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCaa                                            <  2:731534/150‑1 (MQ=255)                                                                                                      ccccGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGTCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACa                                          >  2:213230/1‑150 (MQ=255)                                                                                                         cGGCAGTGCCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCt                                       <  1:1307257/150‑1 (MQ=255)                                                                                                             aGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCa                                                                             >  1:371847/1‑108 (MQ=255)                                                                                                             aGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCa                                                                             <  2:371847/108‑1 (MQ=255)                                                                                                                 gCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGc                                                       <  1:1291545/126‑1 (MQ=255)                                                                                                                 gCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGc                                                       >  2:1291545/1‑126 (MQ=255)                                                                                                                   aGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTg                             >  2:779247/1‑150 (MQ=255)                                                                                                                             cAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGTCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcgcc                   <  1:213230/150‑1 (MQ=255)                                                                                                                                aaaGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCt                >  2:566459/1‑150 (MQ=255)                                                                                                                                 aaGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTg               >  1:407931/1‑150 (MQ=255)                                                                                                                                    a‑gGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTgg                            <  2:997681/134‑1 (MQ=255)                                                                                                                                    a‑gGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTgg                            >  1:997681/1‑134 (MQ=255)                                                                                                                                             gcTGCGCG‑CTcctttcgcgacgcaaggccgctcccacagggcttgcgatatggcttacctgcgtgcttggcatggcaacttcgacctcagccctggccctggacctgaccgaccacgaacaagcaggcaaacgcctttaccgcgagggcg    <  1:228721‑M2/150‑141 (MQ=255)                                                                                                                                               tGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCCCAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCgg  <  1:891731/150‑1 (MQ=255)                                                                                                                                                    g‑cTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCac                                                                   >  2:1074456/1‑80 (MQ=255)                                                                                                                                                    g‑cTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCac                                                                   <  1:1074456/80‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                               GGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGA‑GGGGAGCGCTGCGCG‑CTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGG  >  NC_002947/3951053‑3951343

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