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	33 [0]	[1] 34	[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)	33 (0.560)	32/492	0.5	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

AATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCTGAATGTATAACTGAGATTGATGATTCTGGAAAATAATCACTAAGCCGCATAAATTATCCTTATTTTATTTTTAGGGTACCGTCCGGCGGTAAAAATAAAATACTCGTGCGGTTTAATGAAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAG  >  NC_000913/3079178‑3079649                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                          gctgtcagaaattacactcaccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGtttt                                                                                                                                                                                                                               <  1:326882‑M2/29‑1 (MQ=255)             tacactcaccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtccccAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGa                                                                                                                                                                                                                    >  2:240273‑M2/211‑250 (MQ=255)              acacgcaccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaacgttggttttcgtttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatagcagtgttttatccgtcacaAAGCAATTGAGCAGTGAAACCGGTGTTTTGAAAACAGTGaa                                                                                                                                                                                                                   <  2:172435‑M2/41‑1 (MQ=255)               cactcaccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAatt                                                                                                                                                                                                                 <  1:240273‑M2/43‑1 (MQ=255)                  tcaccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAattatt                                                                                                                                                                                                              >  1:444806‑M2/206‑251 (MQ=255)                    accagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTaa                                                                                                                                                                                                            <  1:349971‑M2/48‑1 (MQ=255)                     ccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcgtgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggattatccgtcacaAAGTAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTTAATTATTAAt                                                                                                                                                                                                           <  2:169081‑M2/49‑1 (MQ=255)                                    cggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACAtt                                                                                                                                                                                            <  1:65420‑M2/64‑1 (MQ=255)                                        aatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggagcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATACCATTGAg                                                                                                                                                                                         <  2:312092‑M2/67‑1 (MQ=255)                                         atttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGtt                                                                                                                                                                                       <  1:320253‑M2/69‑1 (MQ=255)                                                                                  gccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatagcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGc                                                                                                                                               <  2:189287‑M2/109‑1 (MQ=255)                                                                                   ccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCat                                                                                                                                             >  2:356444‑M2/141‑251 (MQ=255)                                                                                     ataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCatat                                                                                                                                           <  1:411061‑M2/113‑1 (MQ=255)                                                                                                   ctttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTg                                                                                                                                                                                                                     <  1:398220‑M2/39‑1 (MQ=255)                                                                                                   ctttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGatat                                                                                                                             >  1:330234‑M2/125‑251 (MQ=255)                                                                                                          tttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTg                                                                                                                       <  1:303208‑M2/133‑1 (MQ=255)                                                                                                             catttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATAGTTATTGCata                                                                                                                   <  1:71602‑M2/137‑1 (MQ=255)                                                                                                                  ttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaaa                                                                                                               <  2:242619‑M2/141‑1 (MQ=255)                                                                                                                        ctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAttt                                                                                                          <  1:140661‑M2/146‑1 (MQ=255)                                                                                                                          gtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtgtg                                                                                                       <  1:44106‑M2/149‑1 (MQ=255)                                                                                                                            aatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATc                                                                                                    <  2:237734‑M2/152‑1 (MQ=255)                                                                                                                            aatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATc                                                                                                    <  2:20677‑M2/152‑1 (MQ=255)                                                                                                                               cagtgtttatgttcggcggttaattgattaaataaagtgtcgagctttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATGTGTGATCTac                                                                                                 <  2:295570‑M2/155‑1 (MQ=255)                                                                                                                                       atgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGcc                                                                                         <  2:384006‑M2/163‑1 (MQ=255)                                                                                                                                                      ttgattacataaagtgtcgagttttggatcgccaatatagttattactcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAACCCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGt                                                                          >  2:182070‑M2/74‑251 (MQ=255)                                                                                                                                                        gattaaataaagtgtagagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTcc                                                                        >  2:172678‑M2/72‑251 (MQ=255)                                                                                                                                                                           ttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAgg                                                     >  2:18455‑M2/53‑251 (MQ=255)                                                                                                                                                                                       aatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAgg                                         >  1:91463‑M2/41‑251 (MQ=255)                                                                                                                                                                                          atagttattaatcaatatcagtggtttatccgacacaAAGCAATTGAGCAGGGAAACAGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCg                                      >  2:131447‑M2/38‑251 (MQ=255)                                                                                                                                                                                                   aatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCtt                                                                            <  1:32743‑M2/176‑1 (MQ=255)                                                                                                                                                                                                   aatcaatatcagcggtttatccgtcacaAAGCAATTGGGCAGGGAAAACGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCtt                                                                            >  2:32743‑M2/29‑204 (MQ=255)                                                                                                                                                                                                          agcattggtttatccgtctcaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCttttt                      <  1:46437‑M2/230‑1 (MQ=255)                                                                                                                                                                                                                       ccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCt          >  1:70358‑M2/9‑250 (MQ=255)                                                                                                                                                                                                                              aaaGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACAGTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcag  <  2:121800/251‑1 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                          AATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCTGAATGTATAACTGAGATTGATGATTCTGGAAAATAATCACTAAGCCGCATAAATTATCCTTATTTTATTTTTAGGGTACCGTCCGGCGGTAAAAATAAAATACTCGTGCGGTTTAATGAAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAG  >  NC_000913/3079178‑3079649

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: ‑