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)	69 (1.760)	57/260	0.0	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)	44 (1.130)	41/258	0.1	100%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	= 4620788	0 (0.000)					coding (1186/1224 nt)	deoB	phosphopentomutase

CATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCAT  >  NC_000913/3584058‑3584343                                                                                                                                          |                                                                                                                                                     cATAACCTGAATCTGAGGCAGCACCTGGCACGTCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAg                                                                                                                                           <  1:192318/149‑1 (MQ=25)  aTAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAg                                                                                                                                           >  2:46928/1‑148 (MQ=35)       cTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtt                                                                                                                                     >  2:26488/1‑149 (MQ=40)           aTCTGAGGCAGCACCTGGCACGGCTGGGAAGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATATGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                                                                                                                                 <  8:68054/149‑1 (MQ=255)                         cTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa                                                                                                                   >  1:146427/1‑149 (MQ=255)                              aCGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTg                                                                                                              >  6:155797/1‑149 (MQ=255)                                 gCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGca                                                                                                           <  1:46928/149‑1 (MQ=255)                                  cTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGcac                                                                                                          <  1:26488/149‑1 (MQ=255)                                   tGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGca                                                                                                           >  7:53802/1‑147 (MQ=255)                                     ggACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAg                                                                                                        >  2:1231/1‑148 (MQ=255)                                                    tCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTc                                                                                        >  4:65482/1‑149 (MQ=255)                                                      gTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGt                                                                                      >  5:40755/1‑149 (MQ=255)                                                       ttCTCAAAATGGGTGGAGCTGCATGACAAAGTCAGCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGTAGTCATTACCGACATTAAAAGGAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGtt                                                                                     <  8:151774/149‑1 (MQ=255)                                                            aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGt                                                                                 <  1:1231/148‑1 (MQ=255)                                                              aaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGAt                                                                              <  3:65482/149‑1 (MQ=255)                                                                    tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTctct                                                                        <  2:146427/149‑1 (MQ=255)                                                                           gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCt                                                                 >  7:162154/1‑149 (MQ=255)                                                                            cATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTc                                                                  <  5:155797/147‑1 (MQ=255)                                                                               gACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGcaca                                                                                                         <  8:18062/105‑1 (MQ=255)                                                                               gACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGcaca                                                                                                         >  7:18062/1‑105 (MQ=255)                                                                                aCAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTaaa                                                             <  8:53802/148‑1 (MQ=255)                                                                                 cAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAAc                                                            >  2:32328/1‑148 (MQ=255)                                                                                        aTCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCa                                                     <  6:40755/148‑1 (MQ=255)                                                                                         tCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTc                                                                      <  1:195699/130‑1 (MQ=255)                                                                                         tCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTc                                                                      >  2:195699/1‑130 (MQ=255)                                                                                                ttATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCtt                                                                       >  8:47297/1‑122 (MQ=255)                                                                                                 tATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAgttg                                           >  5:8768/1‑149 (MQ=255)                                                                                                        aCATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAat                                                                                                   <  5:18300/86‑1 (MQ=255)                                                                                                        aCATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAat                                                                                                   >  6:18300/1‑86 (MQ=255)                                                                                                           tAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATgg                                                  >  5:55202/1‑132 (MQ=255)                                                                                                           tAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATgg                                                  <  6:55202/132‑1 (MQ=255)                                                                                                            aaaaCACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATaataat                                >  8:1822/1‑149 (MQ=255)                                                                                                            aaaaCACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATaataat                                >  1:32577/1‑149 (MQ=255)                                                                                                                       aaTAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGaaa                      >  2:197070/1‑148 (MQ=255)                                                                                                                            gTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAg                                              <  4:147129/119‑1 (MQ=255)                                                                                                                            gTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAg                                              >  3:147129/1‑119 (MQ=255)                                                                                                                             ttGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGaaaagatcg                >  3:71623/1‑142 (MQ=255)                                                                                                                             ttGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGaa                       <  4:71623/141‑1 (MQ=255)                                                                                                                                gAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGCGAATAGTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCtt                                                                       >  4:41414/1‑90 (MQ=255)                                                                                                                                gAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCAt            >  8:151164/1‑149 (MQ=255)                                                                                                                                 aGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCAtt           <  7:151164/149‑1 (MQ=255)                                                                                                                                  gTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATAtt                                                                                              >  2:47571/1‑65 (MQ=255)                                                                                                                                  gTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATAtt                                                                                              <  1:47571/65‑1 (MQ=255)                                                                                                                                   tCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTg          >  4:27017/1‑148 (MQ=255)                                                                                                                                    cATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAAc                                                            >  6:103501/1‑97 (MQ=255)                                                                                                                                    cATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAAc                                                            <  5:103501/97‑1 (MQ=255)                                                                                                                                          cGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCAt  >  1:50482/1‑149 (MQ=255)                                                                                                                                          |                                                                                                                                                     CATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCAT  >  NC_000913/3584058‑3584343

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