Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,910,754	Δ300 bp		[mazF]	[mazF]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2910754	2911053	300	83 [0]	[0] 83	[mazF]	[mazF]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2910753	0 (0.000)	83 (1.150)	68/488	0.0	100%	intergenic (‑68/+3)	mazG/mazF	nucleoside triphosphate pyrophosphohydrolase/endoribonuclease toxin MazF
	?	NC_000913	2911054 =	0 (0.000)					coding (38/336 nt)	mazF	endoribonuclease toxin MazF

GCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGCCGGGATCCGCACCGC  >  NC_000913/2910808‑2911300                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        ccagttcaccgcgaagatcgtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAg                                                                                                                                                                                                                                                     <  1:387718‑M2/5‑1 (MQ=255)    gttctccgcgaagatcgtcaaaaccttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggggtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagttagcaaacggtcgatttgattcattgaattgacctgaaaattgcgggactgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATc                                                                                                                                                                                                                                                  <  2:372826‑M2/8‑1 (MQ=255)       caccgcgaagatcgtcaaaatcttcacggggnatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGcc                                                                                                                                                                                                                                               <  2:261857‑M2/11‑1 (MQ=255)             gcagatagtcaaaatcttcacgggcgatggcgtccagcacttcgaaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatagcgcaggcgctgcataatagtgcgcaaacggtcgatttgattcattgcattgtcctgaaacttgcgggtctgtcaggtggaaacctgtacacagaatagaagtgagttattaaATCAGGTCGCCCACATc                                                                                                                                                                                                                                         >  2:398996‑M2/235‑251 (MQ=255)             gaagatcgtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATAAGATCGCCCATATc                                                                                                                                                                                                                                         >  1:347452‑M2/235‑251 (MQ=255)                    gtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACg                                                                                                                                                                                                                                  >  1:60674‑M2/228‑251 (MQ=255)                        aaatcttcacgggcgatggcgtcaagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATc                                                                                                                                                                                                                              >  1:511697‑M2/224‑251 (MQ=255)                                 cgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCttacc                                                                                                                                                                                                                      >  1:217219‑M2/215‑250 (MQ=255)                                    gcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATAGCGGGTACGTATCGGCttaccatta                                                                                                                                                                                                                  <  2:524142‑M2/40‑1 (MQ=255)                                     ggatggcttccagcacttcgtaggtttctgcaagggtgtaaggcgcaatggtggcaaatgtctgctctgtatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcagggggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCCTATCGGGTACGTATCGGCttaccattac                                                                                                                                                                                                                 <  2:471163‑M2/41‑1 (MQ=255)                                         ggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGa                                                                                                                                                                                                             <  1:52210‑M2/45‑1 (MQ=255)                                          gcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaatggtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGAc                                                                                                                                                                                                            <  1:339902‑M2/46‑1 (MQ=255)                                            gtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACtt                                                                                                                                                                                                          >  1:293296‑M2/204‑251 (MQ=255)                                            gtccagcaattcgtaggattcttcaagggggtaaggcgcaatggtgtcaaatttctgctctttatcccacgggcagccgttttcaggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgacccgactagaagtgagttagtaaATCAGTACGACCATATCGGGTACGTATCGGCGTACCATTCCCAGACtt                                                                                                                                                                                                          >  2:400383‑M2/204‑251 (MQ=255)                                               cagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCt                                                                                                                                                                                                       <  1:115844‑M2/51‑1 (MQ=255)                                                   acttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATc                                                                                                                                                                                                   >  1:86849‑M2/197‑251 (MQ=255)                                                   acttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATc                                                                                                                                                                                                   >  1:242581‑M2/197‑251 (MQ=255)                                                       cgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTc                                                                                                                                                                                               <  1:525080‑M2/59‑1 (MQ=255)                                                          aggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCgg                                                                                                                                                                                             <  1:481808‑M2/61‑1 (MQ=255)                                                           ggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGct                                                                                                                                                                                           <  1:85740‑M2/63‑1 (MQ=255)                                                            gtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgagggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGctc                                                                                                                                                                                          >  1:519666‑M2/188‑251 (MQ=255)                                                                        gtgtaaggcgcaatggggacaaatgtctgctctttttcccacgcgcagcagttttccggatcgcgcaggcgctgcatactagtgagcaaacggtcgatttgtttcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgaggtaggaaATCAGATCGCCCCTATCGGGTACGAATGGGCTTACCCTTACCAGAACTCCTTATCTTTCGGCTCTCCCCAGTCGa‑t                                                                                                                                                                             >  2:441583‑M2/176‑251 (MQ=255)                                                                          gtaatgctcaatggtggcaaattgcttctctntatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccataatagaagtgagttagtaaATCAGATCGCCCATATCGGTTACGTTTCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGa‑tat                                                                                                                                                                           <  2:272670‑M2/78‑1 (MQ=255)                                                                            aaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacgttcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATtc                                                                                                                                                                         <  2:382596‑M2/80‑1 (MQ=255)                                                                              ggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttacggatcgcgcaggcgctgcataaaagcgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtgtaaacctgagaccagaatagcagagagtcagtaaAACAGCTCGCCCATATCGGGAACGTAACGGCTTACCATTACCAGACTTCCTTACCTTTCGGTTATCCCCAGTCGA‑TATtctc                                                                                                                                                                       >  2:581900‑M2/170‑251 (MQ=255)                                                                               gcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCg                                                                                                                                                                      >  1:494757‑M2/169‑251 (MQ=255)                                                                                            aaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGCGGTTTTCCg                                                                                                                                                         >  2:311130‑M2/156‑251 (MQ=255)                                                                                            aaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCg                                                                                                                                                         <  1:298715‑M2/96‑1 (MQ=255)                                                                                            aaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCg                                                                                                                                                         <  1:429205‑M2/96‑1 (MQ=255)                                                                                                 tctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgacttgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagaaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGgcggg                                                                                                                                                    >  2:555543‑M2/151‑249 (MQ=255)                                                                                                    gctctttatcccacgggcagctgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATg                                                                                                                                                 >  2:218830‑M2/148‑251 (MQ=255)                                                                                                    gctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgctttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGt                                                                                                                                                                                                                                 >  1:360059‑M2/148‑172 (MQ=255)                                                                                                    gctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtgctgaaaattgcgggtctgtcaggtggaagcctgtgacccgaatagaagtgatttagtaaATCAGATCTCCCATATCGGGTACGt                                                                                                                                                                                                                                 <  2:360059‑M2/25‑1 (MQ=255)                                                                                                          tatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTg                                                                                                                                           >  1:500008‑M2/142‑251 (MQ=255)                                                                                                           atcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATc                                                                                                                                                                                                                                                       >  1:240017‑M2/141‑143 (MQ=255)                                                                                                           atcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATc                                                                                                                                                                                                                                                       <  2:240017‑M2/3‑1 (MQ=255)                                                                                                             cccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGAcc                                                                                                                                        <  1:541669‑M2/113‑1 (MQ=255)                                                                                                             cacacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagccaacggtcgatttgcttcattgaattgtcctgaaaattgcggatctgtcaggaggtaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCAGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGAGTA‑TCTCGTGGAGGTTTTCCGGCGGTATGTCGTTGAcc                                                                                                                                        >  2:408093‑M2/139‑251 (MQ=255)                                                                                                                  gggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtgtaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTCCGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTc                                                                                                                                   >  1:250318‑M2/134‑251 (MQ=255)                                                                                                                  gggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTc                                                                                                                                   >  1:547121‑M2/134‑251 (MQ=255)                                                                                                                       gccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCaa                                                                                                                              <  2:547121‑M2/123‑1 (MQ=255)                                                                                                                         cgtgttccggatcgcgcaggcgctgcataagagtgagcaaacggtcgattggagtcattgaatgggcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGc                                                                                                                            <  2:250318‑M2/125‑1 (MQ=255)                                                                                                                             ttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCTA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTa                                                                                                                         <  1:540572‑M2/128‑1 (MQ=255)                                                                                                                               ccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTCCCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAAt                                                                                                                      <  1:147261‑M2/131‑1 (MQ=255)                                                                                                                               ccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAAt                                                                                                                      <  2:174826‑M2/131‑1 (MQ=255)                                                                                                                                        gcaggagctgaataatagtgagaaaacggtcgatttgattcattgagttgtcctgcaaattgcgagtctgtaaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGATTACAATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGT‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGCTGGTCAGTTCAAAACGCGTAAATACGGActct                                                                                                             >  2:421922‑M2/112‑251 (MQ=255)                                                                                                                                            gcgctgtagaatagtgagcaaacggtcgatttgattaagtgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaatttggttagtaaATCAGATCGCCCATAGCGGGTACGTATCGGTTGACCATTACCAGACTTCCTTATCTTTCGGCTCTCTCCAGTCGA‑TATGCTCGTGGAGGTTTTCCGGCTTGATTTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTAc                                                                                                         <  2:486482‑M2/144‑1 (MQ=255)                                                                                                                                                tgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtgtaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCAc                                                                                                     >  1:164421‑M2/104‑251 (MQ=255)                                                                                                                                                 gcataatagtgagcaaacgggcgatttgagtcattgaattgtcctgaaaattgcgggtcggtcaggtggaaacctgtgaccagaatagaagtgagttagttaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACTCACt                                                                                                    <  2:86849‑M2/149‑1 (MQ=255)                                                                                                                                                   ataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCCTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGCCCAGTTCAGCAAGCGTAAATACGGGCTCTTTACgcacggg                                                                                                  >  2:410552‑M2/101‑248 (MQ=255)                                                                                                                                                      atagtgagcaagcggtcgatttgattcattgaattgtcctgaaaattgcgggtctggcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTc                                                                                               <  2:415721‑M2/154‑1 (MQ=255)                                                                                                                                                      atagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTc                                                                                               <  2:6232‑M2/154‑1 (MQ=255)                                                                                                                                                      atagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTc                                                                                               <  1:226347‑M2/154‑1 (MQ=255)                                                                                                                                                           gagcaaacggtcgatttgcttcattgcattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGCCCAGTTCAGCAAGCGTAAATACGGGCTCTTTCCGCACTGGCTCcata                                                                                           >  1:145939‑M2/93‑250 (MQ=255)                                                                                                                                                           gagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCaata                                                                                           <  1:202522‑M2/158‑1 (MQ=255)                                                                                                                                                           gagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCaata                                                                                           <  2:348790‑M2/158‑1 (MQ=255)                                                                                                                                                              caaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAAtaatt                                                                                        >  1:248404‑M2/90‑250 (MQ=255)                                                                                                                                                                      cgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGcc                                                                                >  2:223023‑M2/82‑250 (MQ=255)                                                                                                                                                                               tcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAg                                                                       <  2:97914‑M2/178‑1 (MQ=255)                                                                                                                                                                                    gaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAAt                                                                 <  1:434704‑M2/184‑1 (MQ=255)                                                                                                                                                                                       ttgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTc                                                                                                                                                                                 <  1:363227‑M2/73‑1 (MQ=255)                                                                                                                                                                                       ttgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTc                                                                                                                                                                                 >  2:363227‑M2/65‑137 (MQ=255)                                                                                                                                                                                            ctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGCCCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATcttcactt                                                         >  2:209080‑M2/60‑251 (MQ=255)                                                                                                                                                                                              gaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagtgagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTCCCAGACTGCCTTATCTTTCGGCTCTCCCCAGTCTA‑TATTCTCGTAGAGGTTGTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATcttcacttca                                                       <  2:176353‑M2/194‑1 (MQ=255)                                                                                                                                                                                                      ggggtctgtcatgtggaaacctgttaccagaatagaagtgagttagtaaATCAGATCGCCCAGATCGGGTTCGTATCGGCTTACCATTACCAGGCTTCCTTATCGTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAat                                               <  2:500008‑M2/202‑1 (MQ=255)                                                                                                                                                                                                             gtcaggtggaaaccggtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTGCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCa                                          <  2:93657‑M2/207‑1 (MQ=255)                                                                                                                                                                                                               caggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGAtt                                      >  2:192502‑M2/41‑251 (MQ=255)                                                                                                                                                                                                                 ggtggaaacctgttaccagaatagaactgagttagtaaATCATATCGCCCATATCGGGTAAGTATAGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCCGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTACTTTGCCATCCACCAGTTCCATCTTCACTTCCTCCTCAATATTCCGATTg                                     >  2:250356‑M2/39‑250 (MQ=255)                                                                                                                                                                                                                              gaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTa                        <  1:400383‑M2/225‑1 (MQ=255)                                                                                                                                                                                                                               accagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAAc                      <  1:218830‑M2/227‑1 (MQ=255)                                                                                                                                                                                                                               acaagaatagaagtgatttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAAc                      <  2:526323‑M2/227‑1 (MQ=255)                                                                                                                                                                                                                                  agaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGt                    >  1:12387‑M2/22‑250 (MQ=255)                                                                                                                                                                                                                                  agaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTACTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGGTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGt                    >  1:106795‑M2/22‑250 (MQ=255)                                                                                                                                                                                                                                   gaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAg                  <  1:577335‑M2/231‑1 (MQ=255)                                                                                                                                                                                                                                   gaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAg                  <  1:396681‑M2/231‑1 (MQ=255)                                                                                                                                                                                                                                       agaagtgagttagtaaATCAGATAGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGCCgg              >  1:512668‑M2/17‑251 (MQ=255)                                                                                                                                                                                                                                            tgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGCCGGGATcc         >  1:148368‑M2/12‑251 (MQ=255)                                                                                                                                                                                                                                             gagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGC‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGCCACTCTTCACTTCATCATCAATATTCAGANTGCGCGCCTGCATTAACGTAGCCGGGATCCg        >  2:254326‑M2/11‑251 (MQ=255)                                                                                                                                                                                                                                             gagttaggaaAACAGATCGCCCATATCGGGTACGTATCGGCTTTCCATTACCAGGCTTCCTTATCTGTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGTCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTGCATCATCAATATTCAGATTGAGCGCCTGCAGTAACGTAGCCGGGATCCg        <  2:494757‑M2/241‑1 (MQ=255)                                                                                                                                                                                                                                               gttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCCTCAATATTCAGATTGAGCGCCTGCATTAACGTACCCGTGATCCgcc      >  2:13352‑M2/9‑250 (MQ=255)                                                                                                                                                                                                                                                 tagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGCCGGGATCCGCAcc    >  1:185468‑M2/7‑251 (MQ=255)                                                                                                                                                                                                                                                   gtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TCTTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATCTTCAGATTGAGCGCCTGCATTAACGTACCCGGGATCCGCAccgc  >  2:409174‑M2/5‑251 (MQ=255)                                                                                                                                                                                                                                                   gtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGCCGGGATCCGCAccgc  >  1:497978‑M2/5‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        GCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGA‑TATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGCCGGGATCCGCACCGC  >  NC_000913/2910808‑2911300

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