New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3820587 =	36 (0.670)	4 (0.080)	4/244	NT	10.7%	coding (678/966 nt)	endA	endonuclease I
	?	NC_002947	3820646 =	33 (0.650)					coding (737/966 nt)	endA	endonuclease I

TGCTGATGGCCTGGGACAAGCAACACCCCGTGTCGGCCTGGGAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATCACCTGTCAGTCGGCTGCCCGGGCTACAACCAGGT  >  NC_002947/3820505‑3820793                                                                                                                                              |                                                                                                                                                    tgctgATGGCCTGGGACAAGCAACACCCCGTGTCGGCCTGGGAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTg                                                                                                                                              >  8:104677/1‑149 (MQ=255)   ctgATGGCCTGGGACAAGCAACACCCCGTGTCGGCCTGGGAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTgg                                                                                                                                             <  1:311318/148‑1 (MQ=255)   ctgATGGCCTGGGACAAGCAACACCCCGTGTCGGCCTGGGAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTg                                                                                                                                              <  7:104677/147‑1 (MQ=255)      aTGGCCTGGGACAAGCAACTCCCCGTGTCGGCCTGGGAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTaga                                                                                                                                         <  4:248268/149‑1 (MQ=255)          ccggcgACAAGCAACACCCCGTGTCGGCCTGGGAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCg                                                                                                                                     <  3:188596/144‑1 (MQ=255)              ggACAAGCAACACCCCGTGTGGGCCGGGCAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAAGCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGt                                                                                                                                   <  6:231197/147‑1 (MQ=255)                    gCAACACCCCGTGTCGGCCTGGGAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTa                                                                                                                            >  5:293370/1‑148 (MQ=255)                     ggccaggcttggtttggcctctggcttggctgcgtttaccggcaccgcctctaccacaccgcTGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAgcg                                                                                                                                                   <  4:96856‑M1/62‑1 (MQ=255)                     ggccaggcttggtttggcctctggcttggctgcgtttaccggcaccgcctctaccacaccgcTGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAgcg                                                                                                                                                   >  3:96856‑M1/63‑124 (MQ=255)                      aaCACCCCGTGTCGGCCTGGGAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAgcg                                                                                                                                                   <  5:226435/123‑1 (MQ=255)                      aaCACCCCGTGTCGGCCTGGGAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAgcg                                                                                                                                                   >  6:226435/1‑123 (MQ=255)                        cacCCCGTGTCGGCCTGGGAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGc                                                                                                                       <  7:129049/149‑1 (MQ=255)                           cccGTGTCGGCCTGGGAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCagc                                                                                                                    >  7:196170/1‑149 (MQ=255)                              gtgtCGGCCTGGGAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCagcc                                                                                                                   >  4:72233/1‑147 (MQ=255)                                         ggAGAAAGAACGCGACCGGCGCATTGCCGCGATCAGGGGCCATGCCAACCCGTGCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCa                                                                                                      <  8:254437/149‑1 (MQ=255)                                                  aCGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGc                                                                                              <  2:176595/148‑1 (MQ=255)                                                    gcgACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCagccaa                                                                                                                 >  7:41376/1‑127 (MQ=255)                                                    gcgACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCagccaa                                                                                                                 <  8:41376/127‑1 (MQ=255)                                                    gcgACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTg                                                                                           <  6:293370/149‑1 (MQ=255)                                                                   ccGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGATGCCAAACCAAGCCTGGCCAGCGCCGGTtcg                                                                            >  4:159683/1‑149 (MQ=255)                                                                   ccGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTtcg                                                                            >  5:273129/1‑149 (MQ=255)                                                                   ccGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTtc                                                                             <  4:323561/148‑1 (MQ=255)                                                                     gcgATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGt                                                                               >  6:333753/1‑144 (MQ=255)                                                                     gcgATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGt                                                                               <  5:333753/144‑1 (MQ=255)                                                                         tCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGATGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGgc                                                                      <  3:159683/149‑1 (MQ=255)                                                                          cATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGgcg                                                                     >  4:7039/1‑149 (MQ=255)                                                                                cgctgccaaccggcttgtagccggcagtccacttgcgttcgccgcttacgaacgggttggcaGCGGTGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGcc                                                                                        <  4:96856‑M2/62‑1 (MQ=255)                                                                                cgctgccaaccggcttgtagccggcagtccacttgcgttcgccgcttacgaacgggttggcaGCGGTGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGcc                                                                                        >  3:96856‑M2/63‑124 (MQ=255)                                                                                 cATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCt                                                              <  3:72233/149‑1 (MQ=255)                                                                                   tgccaaccggcttgtagccggcagtccacttgcgttcgccgcttacgaacgggttggcaGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTg                                                            >  8:13671‑M2/60‑149 (MQ=255)                                                                                   TGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTg                                                            >  6:273129‑M1/1‑149 (MQ=255)                                                                                   TGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTg                                                            >  5:237555‑M1/1‑149 (MQ=255)                                                                                     ccaaccggcttgtagccggcagtccacttgcgttcgccgcttacgaacgggttggcaGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTg                                                            <  7:13671‑M2/90‑1 (MQ=255)                                                                                     CCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTgg                                                           >  3:7039‑M1/1‑148 (MQ=255)                                                                                      CAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGc                                                          <  8:295810‑M1/148‑1 (MQ=255)                                                                                          CCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTtcggt                                                                          <  6:188217‑M1/128‑1 (MQ=255)                                                                                          CCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTtcggt                                                                          >  5:188217‑M1/1‑128 (MQ=255)                                                                                                         gAACGCAAGTGGACTGCCGGCTACAAGCCGGATGGCAGCGGCGTGGGAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAACCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTa                                      >  7:228005/1‑149 (MQ=255)                                                                                                           aCGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATc                                    <  8:196170/149‑1 (MQ=255)                                                                                                           aCGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATc                                    <  7:114789/149‑1 (MQ=255)                                                                                                                 gTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGACGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATCACCTGt                              <  2:125751/149‑1 (MQ=255)                                                                                                                       tGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATCACCTGTCAGTCg                        >  8:48836/1‑149 (MQ=255)                                                                                                                       tGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATCACCTGTCAGTCg                        >  5:117713/1‑149 (MQ=255)                                                                                                                         ccGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATCACCTGTCAGTCGGc                      >  1:203764/1‑149 (MQ=255)                                                                                                                           ggCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATCACCTGTCAGTCGGCTg                    >  8:64794/1‑149 (MQ=255)                                                                                                                               aCAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATCACCTGTCAGTCGGCTGCCCg                >  1:4807/1‑149 (MQ=255)                                                                                                                                cAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATCACCTGTCAGTCGGCTGCCCgg               <  2:203764/149‑1 (MQ=255)                                                                                                                                 aaGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATCACCTGTCAGTCGGCTGCCCggg              >  4:114380/1‑149 (MQ=255)                                                                                                                                   gCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCAc                                           >  7:308142/1‑118 (MQ=255)                                                                                                                                   gCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCAc                                           <  8:308142/118‑1 (MQ=255)                                                                                                                                        ttggcaGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATCACCTGTCAGTCGGCTGCCCGGGCTACAAc       >  4:217089‑M2/7‑149 (MQ=255)                                                                                                                                            caGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATCACCTGTCAGTCGGCTGCCCGGGCTACAACCAgg   >  2:225887‑M2/3‑149 (MQ=255)                                                                                                                                             aGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATCACCTGTCAGTCGGCTGCCCGGGCTACAACCAGGt  <  6:117713‑M2/148‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    TGCTGATGGCCTGGGACAAGCAACACCCCGTGTCGGCCTGGGAGAAAGAACGCGACCGGCGCATTGCCGCGATCATGGGCCATGCCAACCCGTTCGTAAGCGGCGAACGCAAGTGGACTGCCGGCTACAAGCCGGTTGGCAGCGGCGTGGTAGAGGCGGTGCCGGTAAACGCAGCCAAGCCAGAGGCCAAACCAAGCCTGGCCAGCGCCGGTTCGGTCGGCGCCGTGCTTGGCAACCGCAACAGCCACGTCTATCACCTGTCAGTCGGCTGCCCGGGCTACAACCAGGT  >  NC_002947/3820505‑3820793

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