Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,966,174	IS1 (–) +9 bp	intergenic (+84/‑235)	trxA → / → rho	thioredoxin 1/transcription termination factor

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3583428 =	NA (NA)	30 (1.010)	23/298	0.2	100%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	3966174 =	0 (0.000)					intergenic (+84/‑243)	trxA/rho	thioredoxin 1/transcription termination factor
*	?	NC_000913	= 3584195	NA (NA)	49 (1.640)	34/298	0.0	100%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	= 3966182	0 (0.000)					intergenic (+92/‑235)	trxA/rho	thioredoxin 1/transcription termination factor

AATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACC  >  NC_000913/3584048‑3584341                                                                                                                                                    |                                                                                                                                                   aaTTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGAc                                                                                                                                                 <  1:64071/151‑1 (MQ=31)               aCCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtt                                                                                                                                   >  1:357074/1‑151 (MQ=39)                ccTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttt                                                                                                                                  >  2:40610/1‑151 (MQ=39)                ccTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttt                                                                                                                                  >  2:242875/1‑151 (MQ=39)                     aTCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGa                                                                                                                             >  2:143099/1‑151 (MQ=255)                                    tGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAatat                                                                                                              <  2:465828/151‑1 (MQ=255)                                    tGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAatat                                                                                                              <  2:228178/151‑1 (MQ=255)                                                                     cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGa                                                                             >  2:155408/1‑151 (MQ=255)                                                                              tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTc                                                                    <  2:258918/151‑1 (MQ=255)                                                                                  gCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCAtt                                                                 >  2:89715/1‑150 (MQ=255)                                                                                     gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTa                                                             >  2:296378/1‑151 (MQ=255)                                                                                     gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTa                                                             >  2:344649/1‑151 (MQ=255)                                                                                                    cGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGAc                                              <  1:287561/151‑1 (MQ=255)                                                                                                      ggCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAg                                            >  2:275360/1‑151 (MQ=255)                                                                                                      ggCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAg                                            >  2:415164/1‑151 (MQ=255)                                                                                                              cTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCat                                    <  2:222325/151‑1 (MQ=255)                                                                                                                 aaCATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCataat                                 <  2:384987/151‑1 (MQ=255)                                                                                                                 aaCATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCataat                                 <  2:392110/151‑1 (MQ=255)                                                                                                                   cATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCAtaataa                               >  1:37429/1‑151 (MQ=255)                                                                                                                      aaaaCACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCa                            <  1:66351/151‑1 (MQ=255)                                                                                                                            cTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGa                      >  2:341321/1‑151 (MQ=255)                                                                                                                              aTCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGaaa                    >  1:304327/1‑151 (MQ=255)                                                                                                                                               aTTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAAc   <  2:446209/151‑1 (MQ=255)                                                                                                                                                ttACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAAcc  <  1:308226/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   AATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACC  >  NC_000913/3584048‑3584341

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