Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,077,568	Δ1,832 bp		[cmtA]–cmtB	[cmtA], cmtB

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3077568	3079399	1832	22 [0]	[0] 22	[cmtA]–cmtB	[cmtA],cmtB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3077567	0 (0.000)	22 (0.460)	21/494	0.9	100%	coding (1293/1389 nt)	cmtA	mannitol‑specific PTS enzyme IICB component CmtA
	?	NC_000913	3079400 =	0 (0.000)					intergenic (‑70/+244)	cmtB/tktA	mannitol‑specific PTS enzyme IIA component CmtB/transketolase 1

AAACAGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACACGTTTCACGCGCCCTTCCAGACTGGCATGAGTAACGACGATATCCGCATCCGCAGGCACGTTTTCTATGGCGTAATGTTTTACTTCAATTGCCAGCCCCGCTTTTTCCAGGCGTTTACGGAAGGTGGTCGCGCCCATCGCACTGGAGCCCATCCCCGCATCGCAAACAAAGGCGATACGCTTAACCCTGCTTAACGAGAATGCACCCTCTTGTTTCATCGCCTTAACCGCATTGGCTGACTGAGCA  >  NC_000913/3077338‑3077813                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       aaaCAGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaacc                                                                                                                                                                                                                                   >  1:149905‑M1/1‑230 (MQ=255)     aGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGAGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtg                                                                                                                                                                                                                               >  1:46255‑M1/1‑226 (MQ=255)      gACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgt                                                                                                                                                                                                                              >  1:98878‑M1/1‑225 (MQ=255)            gTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttga                                                                                                                                                                                                                         <  2:182609‑M1/250‑32 (MQ=255)                       cacTCACCAGTAATTCATTATCTGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCAAATATAGTTATTAATCAATATCAGTTGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaatt                                                                                                                                                                                                             <  2:213747‑M1/251‑44 (MQ=255)                                gcaattCATTATCGGTAATGTCTGGCAGATCACATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCCGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGCGTTTTGGATCGCCAATATAGTTATTAATCCCTATCAGTGGTTTATCCGTCACCaagcaattgagcagggaaaccggagctttgcaaacagggaattattaatttc                                                                                                                                                                                                    >  2:191285‑M1/3‑199 (MQ=255)                                      cATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatat                                                                                                                                                                                              >  1:247381‑M1/1‑193 (MQ=255)                                                                                      gCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtaaattattaatttcttatataacatggagttatagatatagcaacaaatgtcaatacgcatatcgtga                                                                                                                                              <  2:337247‑M1/251‑107 (MQ=255)                                                                                                                          ttACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaat                                                                                                          >  2:137739‑M1/1‑109 (MQ=255)                                                                                                                                  gTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtg                                                                                                   <  2:173371‑M1/250‑150 (MQ=255)                                                                                                                                         gtgtTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaa                                                                                           <  1:218241‑M1/251‑158 (MQ=255)                                                                                                                                                                gATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtc                                                                     <  1:152206‑M1/250‑180 (MQ=255)                                                                                                                                                                gATTAAATAAAGTGTCGAGGTTTGGATCGCCAATATAGTTATTAAGCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatttcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtcc                                                                    <  2:33190‑M1/251‑181 (MQ=255)                                                                                                                                                                      ataaaGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaac                                                              <  1:52250‑M1/251‑187 (MQ=255)                                                                                                                                                                                           cGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaa                                          >  2:291571‑M1/1‑44 (MQ=255)                                                                                                                                                                                             ccAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaa                                        >  2:184295‑M1/1‑42 (MQ=255)                                                                                                                                                                                                           aaTCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagcga                         <  2:385584‑M1/251‑224 (MQ=255)                                                                                                                                                                                                                      gTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagcgacctttttaccc              >  2:134970‑M1/1‑17 (MQ=255)                                                                                                                                                                                                                         gTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagagaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgcagcgacctttttacccga            >  1:86400‑M1/1‑14 (MQ=255)                                                                                                                                                                                                                                cGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagcgacctttttacccgaaatgctaa    <  1:184295‑M1/251‑245 (MQ=255)                                                                                                                                                                                                                                 gTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagcgacctttttacccgaaatgctaat   <  1:53650‑M1/251‑246 (MQ=255)                                                                                                                                                                                                                                   cacaaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagcgacctttttacccgaaatgctaatt  <  2:292813‑M1/250‑247 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       AAACAGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACACGTTTCACGCGCCCTTCCAGACTGGCATGAGTAACGACGATATCCGCATCCGCAGGCACGTTTTCTATGGCGTAATGTTTTACTTCAATTGCCAGCCCCGCTTTTTCCAGGCGTTTACGGAAGGTGGTCGCGCCCATCGCACTGGAGCCCATCCCCGCATCGCAAACAAAGGCGATACGCTTAACCCTGCTTAACGAGAATGCACCCTCTTGTTTCATCGCCTTAACCGCATTGGCTGACTGAGCA  >  NC_000913/3077338‑3077813

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