New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	344650 =	38 (0.840)	6 (0.140)	6/242	NT	14.1%	coding (956/1068 nt)	mutY	A/G‑specific adenine glycosylase
	?	NC_002947	344698 =	36 (0.820)					coding (908/1068 nt)	mutY	A/G‑specific adenine glycosylase

GTTCTTCTTGGTATTTGCGGCACATCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGG  >  NC_002947/344502‑344799                                                                                                                                                      |                                                                                                                                                      ggcctgtggggtggtttgtggagcctcccggagctggacgacatcgcccagctcgacgacctggcctaccagcacggcctgcgcctggctggccgccacgccatggacggcctgacccataccttcagccacttccagctggcgatcgag                                                                                                                                                       >  1:494692‑M1/150‑150 (MQ=255) gTTCTTCTTGGTATTTGCGGCACATCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGtcggc                                                                                                                                                      >  1:1076834/1‑150 (MQ=255)  ttcttcTTGGTATTTGCGGCACATCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTGACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGtcggcc                                                                                                                                                     <  2:453189/150‑1 (MQ=255)  ttcttcTTGGTATTTGCGGCACATCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGtcggcc                                                                                                                                                     >  1:422249/1‑150 (MQ=255)  gcctgtggggtggtttgtggagcctcccggagctggacgacatcgcccagctcgacgacctggcctaccagcacggcctgcgcctggctggccgccacgccatggacggcctgacccataccttcagccacttccagctggcgatcgaGc                                                                                                                                                      >  1:652864‑M1/149‑150 (MQ=255)     ttcttGGTATTTGCGGCACATCACGGTGCGGGTCATCGTGCGTCTCCTGTGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTCCCCCTGGGCGGCGAGGCCCAGCCGCGGCGGGGTGGCGAG‑GTTATACCAGTGCCAGtcggcctcg                                                                                                                                                  <  2:194752/150‑1 (MQ=255)       ctnGGTATTTGCGGCACATCACGGTGCNGGTCATCGGGCGNCTCCTGNGNTCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGtcggcctcggc                                                                                                                                                >  1:650905/1‑150 (MQ=255)                   ggCACAGCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGcc                                                                                                                                    <  2:55912/150‑1 (MQ=255)                         tCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCgg                                                                                                                              >  1:1049978/1‑150 (MQ=255)                         tCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCgg                                                                                                                              <  2:748785/150‑1 (MQ=255)                                     gacgacatcgcccagctcgacgacctggcctaccagcacggcctgcgcctggctggccgccacgccatggacggcctgacccataccttcagccacttccagctggcgatcgaGCCATGGCT‑CGTGCTCGCCGACCGGGTCGACGTGCAc                                                                                                                  >  2:802182‑M1/114‑150 (MQ=255)                                                   cTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTGTACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGa                                                                                                    >  2:586159/1‑150 (MQ=255)                                                   cTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGa                                                                                                    >  2:123963/1‑150 (MQ=255)                                                      cgacgacctggcctaccagcacggcctgcgcctggctggccgccacgccatggacggcctgacccataccttcagccacttccagctggcgatcgaGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATc                                                                                                  >  2:502629‑M1/97‑150 (MQ=255)                                                             ttCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTc                                                                                                                                       >  1:1181940/1‑105 (MQ=255)                                                             ttCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTc                                                                                                                                       <  2:1181940/105‑1 (MQ=255)                                                              tggcctaccagcacggcctgcgcctggctggccgccacgccatggacggcctgacccataccttcagccacttccagctggcgatcgaGCCATGGCT‑CGTGCTCGCCGACCGGGTCGACGTGCACGAGc                                                                                                              <  1:705610‑M1/41‑1 (MQ=255)                                                              tggcctaccagcacggcctgcgcctggctggccgccacgccatggacggcctgacccataccttcagccacttccagctggcgatcgaGCCATGGCT‑CGTGCTCGCCGACCGGGTCGACGTGCACGAGc                                                                                                              >  2:705610‑M1/89‑129 (MQ=255)                                                                gtcgGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGaa                                                                                       <  1:708472/150‑1 (MQ=255)                                                                    accagcacggcctgcgcctggctggccgccacgccatggacggcctgacccataccttcagccacttccagctggcgatcgaGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGa                                                                                                    <  1:563227‑M1/52‑1 (MQ=255)                                                                    accagcacggcctgcgcctggctggccgccacgccatggacggcctgacccataccttcagccacttccagctggcgatcgaGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGa                                                                                                    >  2:563227‑M1/83‑134 (MQ=255)                                                                     ccGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGAAAAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGc                                                                                  <  1:616335/150‑1 (MQ=255)                                                                           gTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAgg                                                                            >  2:762325/1‑150 (MQ=255)                                                                               aaGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGCGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGCGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGaacatgat                                                                        <  2:100299/150‑7 (MQ=255)                                                                                agcagcTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATgg                                                                       <  2:1096209/150‑1 (MQ=255)                                                                                agcagcTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATgg                                                                       <  2:746405/150‑1 (MQ=255)                                                                                   agcTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTc                                                                    <  1:123963/150‑1 (MQ=255)                                                                                      tttttCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAg                                                                  <  1:1180920/149‑1 (MQ=255)                                                                                      tttttCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAg                                                                  >  2:1180920/1‑149 (MQ=255)                                                                                      tttttCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGACGG‑GAGTTTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAgg                                                                 <  1:616619/150‑1 (MQ=255)                                                                                                           cccAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGgcca                                            <  2:985008/150‑1 (MQ=255)                                                                                                           cccAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGgcca                                            <  1:463892/150‑1 (MQ=255)                                                                                                            ccAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCAc                                                                                                                                             <  1:925757/52‑1 (MQ=255)                                                                                                            ccAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCAc                                                                                                                                             >  2:925757/1‑52 (MQ=255)                                                                                                                       ggggggGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGcc                                                                                                             >  2:805840/1‑73 (MQ=255)                                                                                                                       ggggTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGcc                                                                                                             <  1:805840/73‑1 (MQ=255)                                                                                                                          gTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATgg                                                                                                         >  2:306441/1‑74 (MQ=255)                                                                                                                          gTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATgg                                                                                                         <  1:306441/74‑1 (MQ=255)                                                                                                                               gAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGTGGCCAGCCAGGCGCAGGCCGTGCTg                        >  1:631438/1‑150 (MQ=255)                                                                                                                               gAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTg                        <  1:762325/150‑1 (MQ=255)                                                                                                                               gAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTg                        >  2:438386/1‑150 (MQ=255)                                                                                                                                   gTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAgg                                >  1:168451/1‑139 (MQ=255)                                                                                                                                   gTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAgg                                <  2:168451/139‑1 (MQ=255)                                                                                                                                        aCCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCa                >  1:452695/1‑150 (MQ=255)                                                                                                                                        aCCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCa                >  1:908755/1‑150 (MQ=255)                                                                                                                                         ccAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCan               <  2:650905/150‑2 (MQ=255)                                                                                                                                          cAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCa                >  2:1203436/1‑148 (MQ=255)                                                                                                                                          cAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCa                <  1:1203436/148‑1 (MQ=255)                                                                                                                                          cAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAgg              >  1:26573/1‑150 (MQ=255)                                                                                                                                                      gcctcggccACGGGCTCGCCGACCGGGTCGACGTGCACGAGACCTGGGTCCCTCGCCACCGGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCCGCCAGGCCCGGCCCTCAAGCTgg  >  2:347706/1‑150 (MQ=255)                                                                                                                                                      GCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTgg  >  2:1049978‑M1/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GTTCTTCTTGGTATTTGCGGCACATCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAG‑GTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGG  >  NC_002947/344502‑344799

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.