Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	1	1,016,059	Δ849 bp		[unknown]–[unknown]	[unknown], [unknown]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	1	1016059	1016907	849	33 [0]	[0] 33	[unknown]–[unknown]	[unknown],[unknown]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	1	= 1016058	0 (0.000)	33 (0.400)	23/426	2.4	100%	coding (74/627 nt)	unknown	PTS system, galactitol‑specific IIC component
	?	1	1016908 =	0 (0.000)					coding (275/1044 nt)	unknown	Galactitol‑1‑phosphate 5‑dehydrogenase (EC 1.1.1.251)

AATGATGTATTAATTGAAGCGACGTCTATCCGTTAAAAGACATTTTGCAGTTAGCGGTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTAACACCCATTGCCAAGCAAGCACGTAAAAAATTACAGGCGCGTTTCGGCGGCGAAGAATTTCTTATCGGCTTGGATCCGGCGCTCTTGCTTGGACACACCTCAGTGGTGTCCGCCAGTTTGATTTTTATTCCGCTGACTATTTTGGTTGCCGTTATCTTACCGGGAAATCAGGTGCTGCCTTTCGGCGATTTGGCCACTATCGGTTTCTTTGTCGCCATGGCGGTAGCGGTGCATCAGGGCAATTT  >  1/1015941‑1016304                                                                                                                      |                                                                                                                                                                                                                                                       ggattattatcggtattttggcgggctatccgttaaaagacattttgcagtTAGCGGTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaa                                                                                                                   >  1:830556‑M1/26‑118 (MQ=255)  gattattatcggtattttggcgggctatccgttaaaagacattttgcagTTAGCGGTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcacc                                                                                                                            >  1:308594‑M1/25‑117 (MQ=255)              tattttggcgggctatccgttaaaaGACATTTTGCAGTTAGCGGTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatga                                                                                                       <  1:522105‑M1/238‑146 (MQ=255)                           tATCCGTTAAAAGACATTTTGCAGTTAGCGGTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaac                                                                                         <  1:772758‑M1/251‑160 (MQ=255)                            aTCCGTTAAAAGACATTTTGCAGTTAGCGGTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaac                                                                                         >  1:880018‑M1/1‑91 (MQ=255)                                gTTAAAAGACATTTTGCAGTTAGCAGTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattc                                                                                                                                                                                                                            <  1:283564‑M1/115‑29 (MQ=255)                                 ttAAAAGACATTTTGCAGTTAGCGGTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgaaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaacc                                                                                                                 <  1:334519‑M1/221‑136 (MQ=255)                                    aaaGACATTTTGCAGTTAGCGGTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacgg                                                                                <  1:144938‑M1/251‑169 (MQ=255)                                    aaaGACATTTTGCAGTTAGCGGTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacgg                                                                                <  1:733984‑M1/251‑169 (MQ=255)                                    aaaGACATTTTGCAGTTAGCGGTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacgg                                                                                <  1:655604‑M1/251‑169 (MQ=255)                                                        ggTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcac                                                                     <  1:431969‑M1/242‑180 (MQ=255)                                                         gTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacgccattaat                                                            >  1:263935‑M1/1‑62 (MQ=255)                                                          tAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacgccattaatc                                                           <  1:402464‑M1/250‑190 (MQ=255)                                                          tAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacgccattaatc                                                           <  1:652046‑M1/250‑190 (MQ=255)                                                                ccGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaa                                                                                                                                                                                                                >  1:177772‑M1/1‑55 (MQ=255)                                                                          gATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggt                                                                               >  1:16234‑M1/1‑45 (MQ=255)                                                                            tgttgtTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacgccattaatctggcgcaggggtgtagcga                                        <  1:805503‑M1/251‑209 (MQ=255)                                                                                gttAATGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacgg                                                                                >  1:598786‑M1/1‑39 (MQ=255)                                                                                   aaTGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacgccattaatctggcgcaggggtgtagcgataaaaat                                 <  1:191295‑M1/251‑216 (MQ=255)                                                                                   aaTGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacgccattaatctggcgcaggggtgtagcgataaaaat                                 <  1:813162‑M1/251‑216 (MQ=255)                                                                                    aTGCCACGCGTCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacgccattaatctggcgcaggggtgtagcgataaaaatg                                <  1:859439‑M1/251‑217 (MQ=255)                                                                                              tCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacgccattaatctggcgcaggggtgtagcgataaaaatgtgattgtggt                      >  1:375359‑M1/1‑25 (MQ=255)                                                                                              tCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtc                                                                          >  1:874760‑M1/1‑25 (MQ=255)                                                                                              tCATTAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattc                                                                                                                                                                                                                            <  1:171925‑M1/53‑29 (MQ=255)                                                                                                  tAAACCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacgccattaatctggcgcaggggtgtagcgataaaaatgtgattgtggtagga                  <  1:80011‑M1/251‑231 (MQ=255)                                                                                                   aaaCCAATCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacg                                                                    >  1:82354‑M1/1‑20 (MQ=255)                                                                                                         aTCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacgccattaatctggcgcaggggtgtagcgataaaaatg                                <  1:68160‑M1/230‑217 (MQ=255)                                                                                                         aTCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacgccattaatctggcg                                                      >  1:74236‑M1/1‑14 (MQ=255)                                                                                                         aTCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtc                                                                          <  1:584577‑M1/188‑175 (MQ=255)                                                                                                          tCATGGACGGTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacgccattaatctggcgcaggggtgtagcgataaaaatgtgattgtggtaggagccggtac          <  1:567799‑M1/251‑239 (MQ=255)                                                                                                                   gTTTgcgacgaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatat                                                                                                                                                                  >  1:408939‑M1/1‑4 (MQ=255)                                                                                                                   gTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctattacggtgggtctgcacgccattaatctggcgcaggggtgtagcgataaaaatgtgattgtggtaggagccggtactattggct  <  1:493545‑M1/250‑247 (MQ=255)                                                                                                                   gTTTgcgaagaatgtcacaaaggatattattcattatgtaaaaactacaatttcatcggttcccgtcgcttcggcgccaatgcggaatatatttgcgtcaacaagaaaaatttggtgaaagtatcaccagattccaaccctttgcatgatgcttttatcgaacctatta                                                                                   >  1:195043‑M1/1‑4 (MQ=255)                                                                                                                      |                                                                                                                                                                                                                                                       AATGATGTATTAATTGAAGCGACGTCTATCCGTTAAAAGACATTTTGCAGTTAGCGGTAAAAACCGCTGCAGTGATGTTGTTAATGCCACGCGTCATTAAACCAATCATGGACGGTTTAACACCCATTGCCAAGCAAGCACGTAAAAAATTACAGGCGCGTTTCGGCGGCGAAGAATTTCTTATCGGCTTGGATCCGGCGCTCTTGCTTGGACACACCTCAGTGGTGTCCGCCAGTTTGATTTTTATTCCGCTGACTATTTTGGTTGCCGTTATCTTACCGGGAAATCAGGTGCTGCCTTTCGGCGATTTGGCCACTATCGGTTTCTTTGTCGCCATGGCGGTAGCGGTGCATCAGGGCAATTT  >  1/1015941‑1016304

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑