New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	255386 =	34 (1.140)	7 (0.250)	5/280	2.0	16.5%	coding (331/1458 nt)	pepD	aminoacyl‑histidine dipeptidase (peptidase D)
	?	NC_000913	255482 =	39 (1.390)					coding (235/1458 nt)	pepD	aminoacyl‑histidine dipeptidase (peptidase D)
Rejected: Frequency below cutoff threshold.

CCAGAACCGCCAGCGCAGAGGCCATACCAATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTT  >  NC_000913/255335‑255631                                                                                                                                                    |                                                                                                                                                      ccAGAACCGCCAGCGCAGAGGCCATACCAATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCatat                                                                                                                                                     <  1:389928/150‑1 (MQ=255)     aaCCGCCAGCGCAGAGGCCATACCAATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATGAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAgg                                                                                                                                                <  1:293578/151‑1 (MQ=255)     aaCCGCCAGCGCAGAGGCCATACCAATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAgg                                                                                                                                                <  1:153085/151‑1 (MQ=255)     aaCCGCCAGCGCAGAGGCCATACCAATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAgg                                                                                                                                                <  2:112695/151‑1 (MQ=255)       gtatggaaaatcgtaaaccggtcgtcttacaggcccacctcgatatggtgccgcagaaaaataacgacaccgtgcatgacttcacgaaagatcctatccagccttatattgatggcgaatgggttaaagcgcgcggcaccaTATCGAgg                                                                                                                                                <  2:426881‑M2/8‑1 (MQ=255)       gtatggaaaatcgtaaaccggtcgtcttacaggcccacctcgatatggtgccgcagaaaaataacgacaccgtgcatgacttcacgaaagatcctatccagccttatattgatggcgaatgggttaaagcgcgcggcaccaTATCGAgg                                                                                                                                                <  1:162279‑M2/8‑1 (MQ=255)       gtatggaaaatcgtaaaccggtcgtcttacaggcccacctcgatatggtgccgcagaaaaataacgacaccgtgcatgacttcacaaaagatcctatccagccttatattgatggcgaatgggttaaagcgcgcggcaccaTATCGAgg                                                                                                                                                <  2:168696‑M2/8‑1 (MQ=255)             gcgcAGAGGCCATACCAATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGc                                                                                                                                           >  1:234970/1‑148 (MQ=255)                   agGCCATACCAATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTaa                                                                                                                                     <  1:288016/148‑1 (MQ=255)                    ggCCATACCAATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAgacg                                                                                                                                 >  1:48166/1‑151 (MQ=255)                      ccATACCAATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAgacgac                                                                                                                               <  2:383146/151‑1 (MQ=255)                                 ccGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGAt                                                                                                                    >  1:96770/1‑151 (MQ=255)                                 ccGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGAt                                                                                                                    <  1:309442/151‑1 (MQ=255)                                      aTCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTcc                                                                                                               <  2:383738/151‑1 (MQ=255)                                       tCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCa                                                                                                              >  1:315978/1‑151 (MQ=255)                                          gCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATa                                                                                                            >  2:192130/1‑150 (MQ=255)                                            cctcgataTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTTGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATAc                                                                                                           <  2:168696‑M1/141‑1 (MQ=255)                                            cctcgataTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATAc                                                                                                           <  2:426881‑M1/141‑1 (MQ=255)                                            cctcgataTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATAc                                                                                                           <  1:162279‑M1/141‑1 (MQ=255)                                                 gCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGcc                                                                                                                                          >  2:199113/1‑113 (MQ=255)                                                 gCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGcc                                                                                                                                          <  1:199110/113‑1 (MQ=255)                                                               gcTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGaa                                                                                      >  2:337324/1‑151 (MQ=255)                                                               aaaaataacgacaccgtgcatgacttcacgaaagatcctatccagccttatattgatggcgaatgggttaaagcgcgcggcaccaTATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGa                                                                                       <  2:414316‑M2/65‑1 (MQ=255)                                                                cTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAAt                                                                                     >  1:57584/1‑151 (MQ=255)                                                                cTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAAt                                                                                     >  2:82249/1‑151 (MQ=255)                                                                           tCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATtac                                                                          <  2:181981/151‑1 (MQ=255)                                                                           accgtgcatgacttcacgaaagatcctatccagccttatattgatggcgaatgggttaaagcgcgcggcaccaTATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATtac                                                                          <  2:368253‑M2/78‑1 (MQ=255)                                                                            ccgtgcatgacttcacgaaagatcctatccagccttatattgatggcgaatgggttaaagcgcgcggcaccaTATCGAGGTGGGCCTATAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATtacc                                                                         <  2:153150‑M2/79‑1 (MQ=255)                                                                              ccATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATtaccta                                                                       >  1:412249/1‑151 (MQ=255)                                                                                aTCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATtacctacc                                                                     <  2:81648/151‑1 (MQ=255)                                                                                  cAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTg                                                                   <  2:171972/151‑1 (MQ=255)                                                                                    atatAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGAt                                                                 >  2:391863/1‑151 (MQ=255)                                                                                    atatAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGAt                                                                 >  1:174373/1‑151 (MQ=255)                                                                                    atatAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGAt                                                                 >  1:169842/1‑151 (MQ=255)                                                                                           gCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGtt                                                           <  1:299954/150‑1 (MQ=255)                                                                                              ggATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGa                                                        >  2:98941/1‑150 (MQ=255)                                                                                              ggATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGa                                                        >  2:290667/1‑150 (MQ=255)                                                                                              ggATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTCTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGa                                                        >  1:450503/1‑150 (MQ=255)                                                                                                aTAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACa                                                      <  2:318850/150‑1 (MQ=255)                                                                                                                 gTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATGTTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCtttt                                    >  2:310318/1‑151 (MQ=255)                                                                                                                 gTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCtttt                                    >  2:418794/1‑151 (MQ=255)                                                                                                                 gTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCtttt                                    >  2:111357/1‑151 (MQ=255)                                                                                                                 gTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCtttt                                    >  2:128854/1‑151 (MQ=255)                                                                                                                 gTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCtttt                                    >  2:349458/1‑151 (MQ=255)                                                                                                                    tgatggcgaatgggttaaagcgcgcggcgccaTATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGcc                                 <  1:403125‑M2/119‑1 (MQ=255)                                                                                                                     tGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGccc                                <  1:261616/151‑1 (MQ=255)                                                                                                                           gtgtCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTCCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCccaaccaa                          >  2:158229/1‑151 (MQ=255)                                                                                                                               cGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAAt                      >  1:451841/1‑151 (MQ=255)                                                                                                                                     ttttCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTc                <  1:328403/151‑1 (MQ=255)                                                                                                                                              gcaccaTATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTg       >  1:267485‑M2/7‑151 (MQ=255)                                                                                                                                                cccaTATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTCTACGAATCATGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTg       >  2:157931‑M2/5‑149 (MQ=255)                                                                                                                                                   aTATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCtt  >  2:12144‑M2/2‑151 (MQ=255)                                                                                                                                                   aTATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCtt  >  2:56469‑M2/2‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CCAGAACCGCCAGCGCAGAGGCCATACCAATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTT  >  NC_000913/255335‑255631

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