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	42 [0]	[0] 42	[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)	42 (0.650)	41/476	0.2	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

CCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCTGAATGTATAACTGAGATTGATGATTCTGGAAAATAATCACTAAGCCGCATAAATTATCCTTATTTTATTTTTAGGGTACCGTCCGGCGGTAAAAATAAAATACTCGTGCGGTTTAATGAAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGA  >  NC_000913/3079170‑3079635                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            aaacagacgctgtcagaaattacactcaccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAAcc                                                                                                                                                                                                                         >  1:315770‑M2/231‑251 (MQ=255)                gaaattacactcaccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAAc                                                                                                                                                                                                           >  2:33590‑M2/216‑250 (MQ=255)                       cactcaccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtggcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAatt                                                                                                                                                                                                   <  2:374582‑M2/43‑1 (MQ=255)                         ctcaccagtaattcattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAattat                                                                                                                                                                                                 >  1:143579‑M2/207‑251 (MQ=255)                                      cattatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTtatat                                                                                                                                                                                    <  1:358926‑M2/58‑1 (MQ=255)                                         tatcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAAc                                                                                                                                                                                 >  1:490797‑M2/191‑251 (MQ=255)                                          atcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAAc                                                                                                                                                                                 <  1:473583‑M2/61‑1 (MQ=255)                                           tcggtaatttctggcagttcaaatacgcgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgcccatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGGGAATTATTAATTTCTTATATAACAt                                                                                                                                                                               >  2:401792‑M2/189‑251 (MQ=255)                                                                 agacgcgcagacggacatcccgcttgccaaaaatagcagcaactttggttttaattaatacctctgcaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttaataatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAAACGGTGTTTCTAAAACAGTGAATTATTAATTGCTTATATAACATTGAGTTATAGATATAacaaca                                                                                                                                                          >  1:459811‑M2/167‑250 (MQ=255)                                                                      cgcagacggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTc                                                                                                                                                    <  1:263259‑M2/90‑1 (MQ=255)                                                                           acggacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttaggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAAAAAATGTCAATAc                                                                                                                                               >  1:394034‑M2/157‑251 (MQ=255)                                                                              gacatcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTTAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGc                                                                                                                                             >  2:47326‑M2/154‑250 (MQ=255)                                                                                 atcccgcttgccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCatat                                                                                                                                         <  1:270677‑M2/101‑1 (MQ=255)                                                                                           ccataaatagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAACACGCATATCGTGATGCat                                                                                                                               <  2:410866‑M2/111‑1 (MQ=255)                                                                                                 atagcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACa                                                                                                                         <  2:266314‑M2/117‑1 (MQ=255)                                                                                                    gcagcaactttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATa                                                                                                                       >  1:133147‑M2/132‑250 (MQ=255)                                                                                                           ctttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAAc                                                                                                                                                                                                           <  2:268421‑M2/35‑1 (MQ=255)                                                                                                           ctttggttttcatttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAAc                                                                                                                                                                                                           >  1:268421‑M2/125‑159 (MQ=255)                                                                                                                      atttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCatat                                                                                                    <  1:62366‑M2/138‑1 (MQ=255)                                                                                                                       tttttacctctgtaatcagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCatat                                                                                                    >  2:47607‑M2/113‑250 (MQ=255)                                                                                                                                       cagtgtttatgttcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTac                                                                                   >  1:393151‑M2/97‑251 (MQ=255)                                                                                                                                                 gtgcggctgttaattgattaaataaagtgtcgagttttggatcgctaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGATCAGGGAACCCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACAGTGAGTTATAGATATAACAACAAATGTCAATACTCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATATACAACACGCCtt                                                                         <  2:443527‑M2/165‑1 (MQ=255)                                                                                                                                                   tcggcggttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaataacagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTTCACATAAATATTTGTGATCTACAACACGCCTTa                                                                        <  1:283381‑M2/166‑1 (MQ=255)                                                                                                                                                         gttaattgattaaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTg                                                                 >  1:392794‑M2/79‑251 (MQ=255)                                                                                                                                                                    aaataaagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCa                                                       >  2:84273‑M2/68‑250 (MQ=255)                                                                                                                                                                         aagtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACgg                                                  >  1:216015‑M2/63‑250 (MQ=255)                                                                                                                                                                           gtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCat                                                                                                                                           >  2:212737‑M2/61‑159 (MQ=255)                                                                                                                                                                           gtgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCat                                                                                                                                           <  1:212737‑M2/99‑1 (MQ=255)                                                                                                                                                                            tgtcgagttttggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACat                                              >  1:427090‑M2/60‑251 (MQ=255)                                                                                                                                                                                       ggatcgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACCGTGAATTATTACTTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGtaat                                   >  2:423574‑M2/49‑251 (MQ=255)                                                                                                                                                                                           cgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATaaaa                               <  1:42720‑M2/207‑1 (MQ=255)                                                                                                                                                                                           cgccaatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATaaa                                <  2:216015‑M2/206‑1 (MQ=255)                                                                                                                                                                                              caatatagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATaaaaa                              >  1:97510‑M2/42‑249 (MQ=255)                                                                                                                                                                                                  atagttattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCg                        <  1:170707‑M2/214‑1 (MQ=255)                                                                                                                                                                                                       tattaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGa                    >  1:264499‑M2/33‑250 (MQ=255)                                                                                                                                                                                                        attaatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGa                    >  2:267847‑M2/32‑249 (MQ=255)                                                                                                                                                                                                          taatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGc                 <  2:404144‑M2/221‑1 (MQ=255)                                                                                                                                                                                                           aatcaatatcagtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCg                >  1:353176‑M2/29‑250 (MQ=255)                                                                                                                                                                                                              caatatcagttgtttatccgtcacaAAACAAGTGAGCAGGGGAACCGGTGTTTTTAAAACAGTGAATTATTAATTTCTTATATTACATTGAGTTATAGATATAACAACCAATGTCAATACGCATATCGTGATTGATATTACATAATTGATATTTATTGCATATAAATATTTGTTATCTACAACACGCCTTATCTATTGCTTGGCCGCAAACGGACATTTCAAGGTAATAAAAAAGGTCGCCGAAGCGACCt            <  1:405807‑M2/226‑1 (MQ=255)                                                                                                                                                                                                                     agtggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTAcc     >  1:189827‑M2/19‑251 (MQ=255)                                                                                                                                                                                                                       tggtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCg   >  1:239357‑M2/17‑251 (MQ=255)                                                                                                                                                                                                                         gtttatccgtcacaAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGa  >  1:393169‑M2/15‑250 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            CCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCTGAATGTATAACTGAGATTGATGATTCTGGAAAATAATCACTAAGCCGCATAAATTATCCTTATTTTATTTTTAGGGTACCGTCCGGCGGTAAAAATAAAATACTCGTGCGGTTTAATGAAAGCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGA  >  NC_000913/3079170‑3079635

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