Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,620,780	IS1 (–) +9 bp	coding (1178‑1186/1224 nt)	deoB →	phosphopentomutase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3583428 =	NA (NA)	30 (1.020)	26/260	0.3	100%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	4620780 =	0 (0.000)					coding (1178/1224 nt)	deoB	phosphopentomutase
*	?	NC_000913	= 3584195	NA (NA)	30 (1.020)	27/258	0.2	100%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	= 4620788	0 (0.000)					coding (1186/1224 nt)	deoB	phosphopentomutase

AGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTA  >  NC_000913/3583281‑3583571                                                                                                                                                    |                                                                                                                                                aGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATATAAATAATGATGATTTCATAAACCCTGATCTACAAGAACTGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAg                                                                                                                                                 <  6:128902/148‑1 (MQ=255)   ccTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTg                                                                                                                                              <  6:134302/149‑1 (MQ=255)    ctctCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGa                                                                                                                                             <  1:9211/149‑1 (MQ=255)    ctctCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGa                                                                                                                                             <  6:20579/149‑1 (MQ=255)     tctcGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGa                                                                                                                                             >  4:115641/1‑148 (MQ=255)     tctcGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGAt                                                                                                                                            >  1:106644/1‑149 (MQ=255)       tcGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGa                                                                                                                                             <  3:115641/146‑1 (MQ=255)       tcGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGAtgc                                                                                                                                          >  2:3352/1‑149 (MQ=255)         gATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGAtgc                                                                                                                                          <  2:106644/147‑1 (MQ=255)              cAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGa                                                                                                                                             <  8:118532/139‑1 (MQ=255)              cAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGa                                                                                                                                             >  7:118532/1‑139 (MQ=255)                   ccTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACt                                                                                                                              <  7:68961/149‑1 (MQ=255)                      tCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGt                                                                                                                                               >  1:137391/1‑129 (MQ=255)                         gTTAAACACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAtt                                                                                                                          >  4:57378/1‑147 (MQ=255)                                  cTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatga                                                                                                                <  5:145604/148‑1 (MQ=255)                                  cTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatga                                                                                                                >  6:145604/1‑148 (MQ=255)                                     tttGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatgatg                                                                                                              <  3:57378/147‑1 (MQ=255)                                        gCCGTAAATGAATTTTATAGAAATAACGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGgtgt                                                                                                          <  3:28368/148‑1 (MQ=255)                                               aTGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAgt                                                                                                                      <  7:5918/129‑1 (MQ=255)                                               aTGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAgt                                                                                                                      >  8:5918/1‑129 (MQ=255)                                                           gAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAtt                                                                                                                          >  4:92370/1‑113 (MQ=255)                                                           gAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAtt                                                                                                                          <  3:92370/113‑1 (MQ=255)                                                               taatGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCt                                                                                  <  1:3352/149‑1 (MQ=255)                                                                                           cAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAg                                                      >  5:30078/1‑149 (MQ=255)                                                                                            aaGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGtt                                                                                                         >  8:53002/1‑97 (MQ=255)                                                                                                      aGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATCGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACggg                                           >  6:112972/1‑149 (MQ=255)                                                                                                       gTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACggg                                           <  6:30078/148‑1 (MQ=255)                                                                                                                   tataGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGc                                    >  2:94275/1‑143 (MQ=255)                                                                                                                   tataGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGc                                    <  1:94275/143‑1 (MQ=255)                                                                                                                    ataGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGc                                                                                      >  2:78381/1‑92 (MQ=255)                                                                                                                                               ttAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTa  >  8:46574/1‑149 (MQ=32)                                                                                                                                                    |                                                                                                                                                AGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTA  >  NC_000913/3583281‑3583571

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.