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

AAACAGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACACGTTTCACGCGCCCTTCCAGACTGGCATGAGTAACGACGATATCCGCATCCGCAGGCACGTTTTCTATGGCGTAATGTTTTACTTCAATTGCCAGCCCCGCTTTTTCCAGGCGTTTACGGAAGGTGGTCGCGCCCATCGCACTGGAGCCCATCCCCGCATCGCAAACAAAGGCGATACGCTTAACCCTGCTTAACGAGAATGCACCCTCTTGTTTCATCGCCTTAACCGCATTGGC  >  NC_000913/3077338‑3077803                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             aaaCAGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaacc                                                                                                                                                                                                                         >  1:315770‑M1/1‑230 (MQ=255)                gAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaac                                                                                                                                                                                                           >  2:33590‑M1/1‑215 (MQ=255)                       cacTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGGCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaatt                                                                                                                                                                                                   <  2:374582‑M1/251‑44 (MQ=255)                         cTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattat                                                                                                                                                                                                 >  1:143579‑M1/1‑206 (MQ=255)                                      cATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatat                                                                                                                                                                                    <  1:358926‑M1/251‑59 (MQ=255)                                         tATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataac                                                                                                                                                                                 >  1:490797‑M1/1‑190 (MQ=255)                                          aTCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataac                                                                                                                                                                                 <  1:473583‑M1/250‑62 (MQ=255)                                           tCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCCATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagggaattattaatttcttatataacat                                                                                                                                                                               >  2:401792‑M1/1‑188 (MQ=255)                                                                 agacgCGCAGACGGACATCCCGCTTGCCAAAAATAGCAGCAACTTTGGTTTTAATTAATACCTCTGCAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTAATAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaacggtgtttctaaaacagtgaattattaattgcttatataacattgagttatagatataacaaca                                                                                                                                                          >  1:459811‑M1/3‑166 (MQ=255)                                                                      cgcAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtc                                                                                                                                                    <  1:263259‑M1/251‑91 (MQ=255)                                                                           aCGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTAGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaaaaaatgtcaatac                                                                                                                                               >  1:394034‑M1/1‑156 (MQ=255)                                                                              gACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaatttagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgc                                                                                                                                             >  2:47326‑M1/1‑153 (MQ=255)                                                                                 aTCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatat                                                                                                                                         <  1:270677‑M1/251‑102 (MQ=255)                                                                                           ccATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaacacgcatatcgtgatgcat                                                                                                                               <  2:410866‑M1/251‑112 (MQ=255)                                                                                                 aTAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattaca                                                                                                                         <  2:266314‑M1/251‑118 (MQ=255)                                                                                                    gcagcaACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacata                                                                                                                       >  1:133147‑M1/1‑131 (MQ=255)                                                                                                           cTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaac                                                                                                                                                                                                           <  2:268421‑M1/159‑36 (MQ=255)                                                                                                           cTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaac                                                                                                                                                                                                           >  1:268421‑M1/1‑124 (MQ=255)                                                                                                                      aTTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatat                                                                                                    <  1:62366‑M1/251‑139 (MQ=255)                                                                                                                       tttttACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatat                                                                                                    >  2:47607‑M1/1‑112 (MQ=255)                                                                                                                                       cAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctac                                                                                   >  1:393151‑M1/1‑96 (MQ=255)                                                                                                                                                 gtgcggctgttaATTGATTAAATAAAGTGTCGAGTTTTGGATCGCTAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgatcagggaacccggtgttttgaaaacagtgaattattaatttcttatataacagtgagttatagatataacaacaaatgtcaatactcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatatacaacacgcctt                                                                         <  2:443527‑M1/248‑166 (MQ=255)                                                                                                                                                   tCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATAACAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttatttcacataaatatttgtgatctacaacacgcctta                                                                        <  1:283381‑M1/250‑167 (MQ=255)                                                                                                                                                         gTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattg                                                                 >  1:392794‑M1/1‑78 (MQ=255)                                                                                                                                                                    aaataaaGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgca                                                       >  2:84273‑M1/1‑67 (MQ=255)                                                                                                                                                                         aaGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacgg                                                  >  1:216015‑M1/1‑62 (MQ=255)                                                                                                                                                                           gtgtCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcat                                                                                                                                           >  2:212737‑M1/1‑60 (MQ=255)                                                                                                                                                                           gtgtCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcat                                                                                                                                           <  1:212737‑M1/159‑100 (MQ=255)                                                                                                                                                                            tgtCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacat                                              >  1:427090‑M1/1‑59 (MQ=255)                                                                                                                                                                                       ggATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaaccgtgaattattactttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaat                                   >  2:423574‑M1/1‑48 (MQ=255)                                                                                                                                                                                           cGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaa                               <  1:42720‑M1/251‑208 (MQ=255)                                                                                                                                                                                           cGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaa                                <  2:216015‑M1/250‑207 (MQ=255)                                                                                                                                                                                              cAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaa                              >  1:97510‑M1/1‑41 (MQ=255)                                                                                                                                                                                                  ataGTTATTAATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcg                        <  1:170707‑M1/251‑215 (MQ=255)                                                                                                                                                                                                       tattaATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccga                    >  1:264499‑M1/1‑32 (MQ=255)                                                                                                                                                                                                        attaATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccga                    >  2:267847‑M1/1‑31 (MQ=255)                                                                                                                                                                                                          taATCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagc                 <  2:404144‑M1/250‑222 (MQ=255)                                                                                                                                                                                                           aaTCAATATCAGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagcg                >  1:353176‑M1/1‑28 (MQ=255)                                                                                                                                                                                                              cAATATCAGTTGTTTATCCGTCACAaaacaagtgagcaggggaaccggtgtttttaaaacagtgaattattaatttcttatattacattgagttatagatataacaaccaatgtcaatacgcatatcgtgattgatattacataattgatatttattgcatataaatatttgttatctacaacacgccttatctattgcttggccgcaaacggacatttcaaggtaataaaaaaggtcgccgaagcgacct            <  1:405807‑M1/251‑227 (MQ=255)                                                                                                                                                                                                                     aGTGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagcgacctttttacc     >  1:189827‑M1/1‑18 (MQ=255)                                                                                                                                                                                                                       tGGTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagcgacctttttacccg   >  1:239357‑M1/1‑16 (MQ=255)                                                                                                                                                                                                                         gTTTATCCGTCACAaagcaattgagcagggaaaccggtgttttgaaaacagtgaattattaatttcttatataacattgagttatagatataacaacaaatgtcaatacgcatatcgtgatgcatattacataattgatatttattgcatataaatatttgtgatctacaacacgccttatctattgcttgtccgcaaacggacatatcaaggtaataaaaaaggtcgccgaagcgacctttttacccga  >  1:393169‑M1/1‑14 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             AAACAGACGCTGTCAGAAATTACACTCACCAGTAATTCATTATCGGTAATTTCTGGCAGTTCAAATACGCGCAGACGGACATCCCGCTTGCCATAAATAGCAGCAACTTTGGTTTTCATTTTTACCTCTGTAATCAGTGTTTATGTTCGGCGGTTAATTGATTAAATAAAGTGTCGAGTTTTGGATCGCCAATATAGTTATTAATCAATATCAGTGGTTTATCCGTCACACGTTTCACGCGCCCTTCCAGACTGGCATGAGTAACGACGATATCCGCATCCGCAGGCACGTTTTCTATGGCGTAATGTTTTACTTCAATTGCCAGCCCCGCTTTTTCCAGGCGTTTACGGAAGGTGGTCGCGCCCATCGCACTGGAGCCCATCCCCGCATCGCAAACAAAGGCGATACGCTTAACCCTGCTTAACGAGAATGCACCCTCTTGTTTCATCGCCTTAACCGCATTGGC  >  NC_000913/3077338‑3077803

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