New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 665394	45 (0.550)	4 (0.050)	4/256	NT	8.3%	coding (920/1113 nt)	mrdB	SEDS family protein MrdB
	?	NC_000913	= 665423	46 (0.610)					coding (891/1113 nt)	mrdB	SEDS family protein MrdB

TAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGC  >  NC_000913/665276‑665571                                                                                                                                                    |                                                                                                                                                     tAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTct                                                                                                                                                    >  1:375025‑M2/1‑148 (MQ=255)  aGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTctg                                                                                                                                                   <  1:364768‑M2/150‑4 (MQ=255)    cgcCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTctggc                                                                                                                                                 <  2:694764‑M2/150‑6 (MQ=255)     gcCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTctggcg                                                                                                                                                <  2:940529‑M2/150‑7 (MQ=255)       cGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAATGTGGTTTGCGCTctggcggc                                                                                                                                              >  2:695329‑M2/1‑142 (MQ=255)       cGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTctggcggc                                                                                                                                              >  1:1092570‑M2/1‑142 (MQ=255)       cGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATAGCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTctggcggc                                                                                                                                              >  1:666401‑M2/1‑142 (MQ=255)               cATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCAT‑CCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAg                                                                                                                                     <  2:667810/150‑1 (MQ=255)                 tAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGc                                                                                                                                    >  2:792409/1‑150 (MQ=255)                 tAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGc                                                                                                                                    >  1:114708/1‑150 (MQ=255)                 aaaCTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGc                                                                                                                                    <  2:1417665/149‑1 (MQ=255)                       aCCAGTTGGAGCGGAACCCCTACAACCGCCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAATCCGCCAGCCATGACTCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCgcgc                                                                                                                              <  1:707349/150‑1 (MQ=255)                        ccAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCAACCAAAGGTGGTTTGCGCTCTGGCGGCTATCGAGAGCCCGCGCa                                                                                                                             <  2:1156683/150‑1 (MQ=255)                          aGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATg                                                                                                                           >  2:425849/1‑150 (MQ=255)                            tGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTTCCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGCCGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGAt                                                                                                                         >  2:1251082/1‑150 (MQ=255)                               gAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATcag                                                                                                                      >  2:72386/1‑150 (MQ=255)                               gAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATcag                                                                                                                      <  1:1117199/150‑1 (MQ=255)                                aGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATcagc                                                                                                                     <  2:876894/150‑1 (MQ=255)                                      aCCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATg                                                                                                               <  2:1122712/150‑1 (MQ=255)                                       ccccTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGt                                                                                                              <  2:393287/150‑1 (MQ=255)                                        cccTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTa                                                                                                             >  2:832548/1‑150 (MQ=255)                                            aCAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTagaga                                                                                                         <  1:940529/150‑1 (MQ=255)                                                  ggCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAgagcgagc                                                                                                   <  1:304689/150‑1 (MQ=255)                                                             gCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTGTCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGcagaatcagaat                                                                                        >  2:1393690/1‑150 (MQ=255)                                                             gCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGcagaatcagaat                                                                                        >  2:355704/1‑150 (MQ=255)                                                                taccataccaATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGACATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGcc                                                                                     >  1:153336/1‑150 (MQ=255)                                                                  ccataccaATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCa                                                                                   >  1:657599/1‑150 (MQ=255)                                                                      accaATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTaa                                                                               >  1:923996/1‑150 (MQ=255)                                                                        caATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATc                                                                             <  2:114708/150‑1 (MQ=255)                                                                            atTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAg                                                                         <  1:872671/150‑1 (MQ=255)                                                                               tACGAAGACATAAACGAATAATATAAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGctc                                                                      <  1:686031/150‑1 (MQ=255)                                                                                 cGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCtt                                                                    >  1:347229/1‑150 (MQ=255)                                                                                      aCATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGcc                                                               <  1:487165/150‑1 (MQ=255)                                                                                       cATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCa                                                              >  1:265276/1‑150 (MQ=255)                                                                                           aaCGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTAc                                                          >  2:314987/1‑150 (MQ=255)                                                                                             cGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACcg                                                        >  1:1312168/1‑150 (MQ=255)                                                                                                      tCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATaa                                               >  2:935982/1‑150 (MQ=255)                                                                                                      tCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATaa                                               <  1:792409/150‑1 (MQ=255)                                                                                                        aGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAg                                             >  2:786738/1‑150 (MQ=255)                                                                                                        aGCATTAACCCGCCAGCCATCACGCGACCAAAGGTGGATTGCGCTCTGGCGGCTAACCACAGCCCGCGCATCATCAGCAGCATGTAGAGAGAGAGCAGAATCAGAATGCCCACAAATCCCAGCTCTTCCGCCAGTACAGCGAAGATAAAg                                             >  2:756070/1‑150 (MQ=255)                                                                                                         gCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGt                                            >  2:37813/1‑150 (MQ=255)                                                                                                            ttAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTctggcggcttaatgctgatattattcgtttatgtcttcgtaaatattggtatggtaagcggtattctgctggttgtaggggttccgctcccagtggtcagttatggagg                                         <  1:243359‑M2/150‑110 (MQ=255)                                                                                                            ttAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAg                                         <  1:832548/150‑1 (MQ=255)                                                                                                              aaGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTa                                       >  1:1097296/1‑150 (MQ=255)                                                                                                               aGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTctggcggcttaatgctgatattattcgtttatgtcttcgtaaatattggtatggtaagcggtattccgccggttgtaggggttccgctcccactggtcagttatggaggatc                                      <  1:354652‑M2/150‑113 (MQ=255)                                                                                                               aGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTAt                                      >  2:35680/1‑150 (MQ=255)                                                                                                                gccgccAGCCATGACGCGACCAAAGGTGGTTTGCGCTctggcggcttaatgctgatattattcgtttatgtcttcgtaaatattggtatggtaagcggtattctgccggttgtaggggttccgctcccactggtcagttatggaggatcg                                     <  1:1214473‑M2/150‑114 (MQ=255)                                                                                                                       gccaTGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCg                              >  2:1178587/1‑150 (MQ=255)                                                                                                                          aTGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCgggg                           <  1:297604/150‑1 (MQ=255)                                                                                                                                cgACCAAAGGTGGTTTGCGCTCTGGCGGCTGTCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAAt                     <  1:1393690/150‑1 (MQ=255)                                                                                                                                  aCCAAAGGTGGTTTGCGCTctggcggcttaatgctgatattattcgtttatgtcttcgtaaatattggtatggtaagcggtattctgccggttgtaggggttccgctcccactggtcagttatggaggatcggcgctaattgtgctgatg                   >  2:1214473‑M2/1‑19 (MQ=255)                                                                                                                                  aCCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTc                   <  2:364768/150‑1 (MQ=255)                                                                                                                                   ccAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCa                  <  1:694764/150‑1 (MQ=255)                                                                                                                                    cAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCaa                 <  1:72386/150‑1 (MQ=255)                                                                                                                                     aaaGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAg                <  2:347229/150‑1 (MQ=255)                                                                                                                                       aGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCt              <  1:247073/150‑1 (MQ=255)                                                                                                                                        ggtggtTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGAAAAAGTCAGAATGGCGTTCGGGGAGAAATTCAAGCtg             >  2:207150/1‑150 (MQ=255)                                                                                                                                                   ctctGGCGGCTATCCACAGCCCGCGAATGATAAGAAGAATGTAGAGAGCGAGCAGAATCAGAATTCGCACAAATCCCAGCTATTCCGCCAGTACAGCGAAGATAAAGTCAGTATGGCGTTCCGGGAGAAATTCAAGCTGTGACTGAGGGc  >  2:915301/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCGCTTACCATACCAATATTTACGAAGACATAAACGAATAATATCAGCATTAAGCCGCCAGCCATGACGCGACCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGC  >  NC_000913/665276‑665571

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