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	28 [0]	[2] 30	[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)	28 (0.470)	27/492	0.8	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

CACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACACGTTTCACGCGCCCTTCCAGACTGGCATGAGTAACGACGATATCCGCATCCGCAGGCACGTTTTCTATGGCGTAATGTTTTACTTCAATTGCCAGCCCCGCTTTTTCCAGGCGTTTACGGAAGGTGGTCGCGCCCATCGCACTGGAGCCCATCCCCGCATCGCAAACAAAGGCGATACGCTTAACCCTGCTTAACGAGAATGCACCCTCTTGTTTCATCGCCTTAACCGCATTGGCTGACTGAGCA  >  NC_000913/3077364‑3077813                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       cACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattatt                                                                                                                                                                                                          <  2:419671‑M1/250‑47 (MQ=255) cACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattat                                                                                                                                                                                                           <  1:49621‑M1/249‑46 (MQ=255)    cAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatt                                                                                                                                                                                                      <  1:114284‑M1/251‑51 (MQ=255)        aaTTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttctt                                                                                                                                                                                                  <  1:330218‑M1/251‑55 (MQ=255)                aTCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataac                                                                                                                                                                                           >  1:225906‑M1/1‑189 (MQ=255)                           cTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttat                                                                                                                                                                                 >  2:221459‑M1/1‑178 (MQ=255)                                     aaaTACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaac                                                                                                                                                                     >  1:359862‑M1/1‑168 (MQ=255)                                             gcAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtca                                                                                                                                                             <  1:238962‑M1/251‑92 (MQ=255)                                                 aCGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatctaacaacaaatgtcaatac                                                                                                                                                         >  2:29051‑M1/1‑156 (MQ=255)                                                         cccGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcg                                                                                                                                                 >  1:35827‑M1/1‑148 (MQ=255)                                                               tGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgat                                                                                                                                             <  2:35827‑M1/249‑108 (MQ=255)                                                                            agcaACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataaat                                                                                                                              >  2:229031‑M1/1‑129 (MQ=255)                                                                                                 tACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaat                                                                                                          <  1:29051‑M1/250‑143 (MQ=255)                                                                                                                 gtTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaac                                                                                          >  2:402269‑M1/1‑92 (MQ=255)                                                                                                                               gTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattg                                                                           >  2:44135‑M1/1‑78 (MQ=255)                                                                                                                                     tGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtc                                                                     >  2:409197‑M1/1‑72 (MQ=255)                                                                                                                                           aataaaGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaa                                                               <  1:343458‑M1/251‑186 (MQ=255)                                                                                                                                                     cGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatc                                                     >  1:431240‑M1/1‑56 (MQ=255)                                                                                                                                                        gTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaag                                                  <  1:357478‑M1/251‑199 (MQ=255)                                                                                                                                                             ggATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatcttttgcttgtccgcaaacgcacatatcagggtaat                                             <  2:442826‑M1/251‑204 (MQ=255)                                                                                                                                                                    cAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgcaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaa                                        >  1:158966‑M1/1‑41 (MQ=255)                                                                                                                                                                      atatAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcaggggaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggt                                    <  2:259810‑M1/251‑213 (MQ=255)                                                                                                                                                                      ataaagttaTTAATCAATATCAGTGGTTTATCCGTCACAacgcaattgagaagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataactatttgtgatctacaacacactttatctattgcttgtcagcaaacggacatatcatggcaataaaaaaggt                                    >  2:279512‑M1/5‑39 (MQ=255)                                                                                                                                                                               ttaATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagc                           >  2:323282‑M1/1‑30 (MQ=255)                                                                                                                                                                                      atatCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagcgaccttt                    >  1:3434‑M1/1‑23 (MQ=255)                                                                                                                                                                                           aGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagcgacctttttacc               >  2:146829‑M1/1‑18 (MQ=255)                                                                                                                                                                                                   aTCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagcgacctttttacccgaaatgc       >  2:455903‑M1/1‑10 (MQ=255)                                                                                                                                                                                                         cacaaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagcgacctttttacccgaaatgctaatt  >  2:460453‑M1/1‑4 (MQ=255)                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       CACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACACGTTTCACGCGCCCTTCCAGACTGGCATGAGTAACGACGATATCCGCATCCGCAGGCACGTTTTCTATGGCGTAATGTTTTACTTCAATTGCCAGCCCCGCTTTTTCCAGGCGTTTACGGAAGGTGGTCGCGCCCATCGCACTGGAGCCCATCCCCGCATCGCAAACAAAGGCGATACGCTTAACCCTGCTTAACGAGAATGCACCCTCTTGTTTCATCGCCTTAACCGCATTGGCTGACTGAGCA  >  NC_000913/3077364‑3077813

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