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

CCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCTGAATGTATAACTGAGATTGATGATTCTGGAAAATAATCACTAAGCCGCATAAATTATCCTTATTTTATTTTTAGGGTACCGTCCGGCGGTAAAAATAAAATACTCGTGCGGTTTAATGAAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATT  >  NC_000913/3079170‑3079645                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      aaacagacgctgtcagaaattacactcaccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAAcc                                                                                                                                                                                                                                   >  1:149905‑M2/231‑251 (MQ=255)     agacgctgtcagaaattacactcaccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagagtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGgtg                                                                                                                                                                                                                               >  1:46255‑M2/227‑251 (MQ=255)      gacgctgtcagaaattacactcaccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGgtgt                                                                                                                                                                                                                              >  1:98878‑M2/226‑251 (MQ=255)            gtcagaaattacactcaccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGa                                                                                                                                                                                                                         <  2:182609‑M2/31‑1 (MQ=255)                       cactcaccagtaattcattatctgtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgcaaatatagttattaatcaatatcagttgtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAatt                                                                                                                                                                                                             <  2:213747‑M2/43‑1 (MQ=255)                                gcaattcattatcggtaatgtctggcagatcacatacgcgcagacggacatcccgcttgccataaatagccgcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgcgttttggatcgccaatatagttattaatccctatcagtggtttatccgtcaccAAGCAATTGAGCAGGGAAACCGGAGCTTTGCAAACAGGGAATTATTAATTTc                                                                                                                                                                                                    >  2:191285‑M2/200‑251 (MQ=255)                                      cattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTtatat                                                                                                                                                                                              >  1:247381‑M2/194‑251 (MQ=255)                                                                                      gcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTAAATTATTAATTTCTTATATAACATGGAGTTATAGATATAGCAACAAATGTCAATACGCATATCGTGa                                                                                                                                              <  2:337247‑M2/106‑1 (MQ=255)                                                                                                                          ttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAat                                                                                                          >  2:137739‑M2/110‑251 (MQ=255)                                                                                                                                  gtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtgtg                                                                                                   <  2:173371‑M2/149‑1 (MQ=255)                                                                                                                                         gtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTacaa                                                                                           <  1:218241‑M2/157‑1 (MQ=255)                                                                                                                                                                gattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTc                                                                     <  1:152206‑M2/179‑1 (MQ=255)                                                                                                                                                                gattaaataaagtgtcgaggtttggatcgccaatatagttattaagcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATTTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTcc                                                                    <  2:33190‑M2/180‑1 (MQ=255)                                                                                                                                                                      ataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAAc                                                              <  1:52250‑M2/186‑1 (MQ=255)                                                                                                                                                                                           cgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATaaa                                          >  2:291571‑M2/45‑250 (MQ=255)                                                                                                                                                                                             ccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATaaaaa                                        >  2:184295‑M2/43‑250 (MQ=255)                                                                                                                                                                                                           aatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGa                         <  2:385584‑M2/223‑1 (MQ=255)                                                                                                                                                                                                                      gtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTAccc              >  2:134970‑M2/18‑251 (MQ=255)                                                                                                                                                                                                                         gtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGAGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGCAGCGACCTTTTTACCCGa            >  1:86400‑M2/15‑250 (MQ=255)                                                                                                                                                                                                                                cgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTaa    <  1:184295‑M2/244‑1 (MQ=255)                                                                                                                                                                                                                                 gtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAAt   <  1:53650‑M2/245‑1 (MQ=255)                                                                                                                                                                                                                                   cacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAAtt  <  2:292813‑M2/246‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      CCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCTGAATGTATAACTGAGATTGATGATTCTGGAAAATAATCACTAAGCCGCATAAATTATCCTTATTTTATTTTTAGGGTACCGTCCGGCGGTAAAAATAAAATACTCGTGCGGTTTAATGAAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATT  >  NC_000913/3079170‑3079645

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