Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,546,574	IS2 (+) +5 bp	coding (16‑20/1317 nt)	gntT →	high‑affinity gluconate transporter

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2996361 =	NA (NA)	47 (1.810)	42/492	0.0	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	3546574 =	0 (0.000)					coding (16/1317 nt)	gntT	high‑affinity gluconate transporter
*	?	NC_000913	3186096 =	NA (NA)	31 (1.200)	29/488	0.1	100%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 3546578	0 (0.000)					coding (20/1317 nt)	gntT	high‑affinity gluconate transporter

TAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGC  >  NC_000913/2996182‑2996606                                                                                                                                                                                    |                                                                                                                                                                                                                                                      tAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAACGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTcc                                                                                                                                                                                >  1:26851/1‑251 (MQ=255)                          gTGCCTTTTCATTTTTTTCCAGCGGCATAGCTGATAAATAATCCCTTCGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCACGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCAc                                                                                                                                                       >  2:179320/1‑251 (MQ=255)                             ccTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAt                                                                                                                                                    <  2:52871/251‑1 (MQ=255)                               ttttCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGc                                                                                                                                                  >  1:263038/1‑251 (MQ=255)                                        atttCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAt                                                                                                                                         <  2:263038/251‑1 (MQ=255)                                          ttCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTACAAAAAAACCCGAGTTCTGCAAACCTTTCTCGATACAAGTGACAACCACTGAGTCAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCAGCTGCCGCAGATATTCACGTGGCGAGCGATAACCCAGCGCACTCTGCGGATGCCATTc                                                                                                                                       >  2:263138/1‑251 (MQ=37)                                                     taGCTGATAAATCATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCt                                                                                                                            >  2:266991/1‑251 (MQ=255)                                                                      cTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATaa                                                                                                                                >  1:121086/1‑230 (MQ=255)                                                                      cTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATaa                                                                                                                                <  2:121086/230‑1 (MQ=255)                                                                                                                 aCCTTTCTCGATAAACGTGACAACCACTGAGTCAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTCCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACCAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATACTGCTCGACCGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTg                                                                 >  1:163139/1‑250 (MQ=255)                                                                                                                                   gACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCg                                              >  2:240939/1‑251 (MQ=255)                                                                                                                                         cACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCa                                        <  1:242462/251‑1 (MQ=255)                                                                                                                                                                   aaaTTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTActc              >  1:108706/1‑251 (MQ=35)                                                                                                                                                                    aaTTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTActct             <  2:90950/251‑1 (MQ=35)                                                                                                                                                                            ttGTTTATTGGATTTGCCCCAATATTTCCAGACATCTGTTATCTCTTAACCCATTACAAGCCCGCTGCCGCATATATTCCCGTGGAAAGCGATAAACCAGCGCACTATGCGGATGCAACTCGTTTTAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTATGGTTTGGGCATGATACTGATGTAGACACGCTTTAACGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTc     >  2:187654/1‑251 (MQ=11)                                                                                                                                                                               tttATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGc  >  1:140777/1‑251 (MQ=31)                                                                                                                                                                                    |                                                                                                                                                                                                                                                      TAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGC  >  NC_000913/2996182‑2996606

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