Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,361,226	Δ9,811 bp		[gltF]–[nanE]	[gltF], yhcA, yhcD, yhcE, insH10, yhcE, yhcF, yhcG, nanQ, nanK, [nanE]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3361226	3371036	9811	38 [0]	[0] 40	[gltF]–[nanE]	[gltF],yhcA,yhcD,yhcE,insH10,yhcE,yhcF,yhcG,nanQ,nanK,[nanE]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3361225	0 (0.000)	38 (0.520)	35/492	0.6	100%	coding (50/765 nt)	gltF	periplasmic protein GltF
	?	NC_000913	3371037 =	0 (0.000)					intergenic (‑1/+47)	nanE/nanT	putative N‑acetylmannosamine‑6‑phosphate 2‑epimerase/N‑acetylneuraminate:H(+) symporter

CGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTACCC  >  NC_000913/3370833‑3371282                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      tcacatgatcgtaactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAgg                                                                                                                                                                                                          <  2:164276‑M2/46‑1 (MQ=255)     atgatcgtaactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTa                                                                                                                                                                                                      >  1:537212‑M2/201‑250 (MQ=255)         tcgtaactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACttt                                                                                                                                                                                                 >  1:511209‑M2/197‑251 (MQ=255)             aactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGt                                                                                                                                                                                             <  1:164358‑M2/59‑1 (MQ=255)                     aattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTa                                                                                                                                                                                     <  1:44124‑M2/67‑1 (MQ=255)                               ataaatacatcccagcaagttggctatttttnagctgtttatttatttacattcacaaacagatatatatacttcttcagctttttattttaaaaaataaatttgttattaaggttttgttcatatgtttttcaaaaagaacctcacaacatctgctatttgtgcagcactttcACAGCTGCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGCCTAAATCGttttt                                                                                                                                                                           <  2:277317‑M2/77‑1 (MQ=255)                                               aagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCTTTTTTGGCGCTGCCAAacggc                                                                                                                                                           <  2:236255‑M2/93‑1 (MQ=255)                                                 gttggcttattatgaacgaattatttaagtacattcgcaaacatatatagatactgcttcatcgatttattgtataaaataaatttggtattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctattggtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATGATCTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggcac                                                                                                                                                         <  2:226821‑M2/95‑1 (MQ=255)                                                   tggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggcacgg                                                                                                                                                       >  1:97573‑M2/155‑251 (MQ=255)                                                              gaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGaa                                                                                                                                             <  2:163730‑M2/107‑1 (MQ=255)                                                                                                  tatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacg                                                                                                                                                             >  2:365562‑M2/108‑198 (MQ=255)                                                                                                     actgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTcc                                                                                                     >  1:518844‑M2/105‑251 (MQ=255)                                                                                                         cttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGCTTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGgc                                                                                                 >  2:121057‑M2/101‑251 (MQ=255)                                                                                                            cagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGgcgc                                                                                               >  2:191368‑M2/98‑250 (MQ=255)                                                                                                                    tattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgttcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTTGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCACCCAACGc                                                                                      >  1:275055‑M2/90‑251 (MQ=255)                                                                                                                          aaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAc                                                                                 <  2:573595‑M2/167‑1 (MQ=255)                                                                                                                           aaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAcgc                                                                               >  1:101307‑M2/83‑251 (MQ=255)                                                                                                                            aaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCCCCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAcgcg                                                                              >  2:381980‑M2/82‑251 (MQ=255)                                                                                                                             aaataaatttgttattaccgatatgttcatctgattttcgaaaagaacatcacaacagctgctatttgtgcagcaatttcACATCTTCCCTTAGCGAAAGGCCCGATAACTAGCCCAGGCCACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCTCAACCAAAGCTGAAcgcg                                                                              >  2:475578‑M2/81‑250 (MQ=255)                                                                                                                                          attaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGCCTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCACCCAACGCTGAACGCGAGAAGGCATATcc                                                                 >  2:232154‑M2/68‑250 (MQ=255)                                                                                                                                                       tcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATcagc                                                     >  1:259562‑M2/55‑249 (MQ=255)                                                                                                                                                        catatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAgg                                                  >  1:333084‑M2/54‑251 (MQ=255)                                                                                                                                                             gtttttaaaaaagaacctcacaaccgctgctttttgtgcagctctttaCCATCTTCCCTTAGCGAAAGGCCCGGTGCATAGACCGGGCATCATGATTAACGTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACTCAACGCTTCCGTGCGCAACCAATGCTGAACGCGATAAGGCATCTCCAGCCCTATCCGCAGGATca                                              <  2:545920‑M2/202‑1 (MQ=255)                                                                                                                                                                   caaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCNTATCCAGCCCAATCAGCAGGATCACCACGa                                        >  2:270798‑M2/43‑250 (MQ=255)                                                                                                                                                                    aaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggcac                                                                                                                                                         <  2:182779‑M2/95‑1 (MQ=255)                                                                                                                                                                    aaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggcac                                                                                                                                                         >  1:182779‑M2/42‑136 (MQ=255)                                                                                                                                                                      aaagaacctcacaacagctgctatttgcgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAACTCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAAcgg                                    >  2:28742‑M2/40‑250 (MQ=255)                                                                                                                                                                        agaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTc                                   >  2:116970‑M2/38‑250 (MQ=255)                                                                                                                                                                                 caacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACg                         >  2:54638‑M2/29‑251 (MQ=255)                                                                                                                                                                                 caacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAAGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACg                         >  1:549556‑M2/29‑251 (MQ=255)                                                                                                                                                                                  aacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTa                                                                                                                                                                                     >  2:372437‑M2/28‑94 (MQ=255)                                                                                                                                                                                  aacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTa                                                                                                                                                                                     <  1:372437‑M2/67‑1 (MQ=255)                                                                                                                                                                                    cagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACgaga                      <  1:83531‑M2/226‑1 (MQ=255)                                                                                                                                                                                      gctgctatttgtgcagcactttcACAGCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCGGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGc                    <  1:269308‑M2/228‑1 (MQ=255)                                                                                                                                                                                       ctgctatttgtgcagcactttcACATCTTCCCTTAGCGCAAGGCCCGGTACATAGACCGGGCAACAGGCTTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCAGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGACCGTCAGACCGcacgagaccg                   >  2:215118‑M2/23‑248 (MQ=255)                                                                                                                                                                                         gctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGCCTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGANCACCACGAACGTCAGACTGAACGAGAGCGAt                 >  2:283616‑M2/21‑251 (MQ=255)                                                                                                                                                                                           tatttgagcagcactttcACATCTTCCCCTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTATGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAAAGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATAAGCACGCAACGCTTCCGGGCGCCACCAACGCTGACAGCGAGAAGGCATATCCAGCCCAATCAGCAGAATAACAACGAAAGTACGAATGAACGAGAGCGATg                >  2:384083‑M2/19‑250 (MQ=255)                                                                                                                                                                                                          ttcACATCTTCCCTTAGCGAAATGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCCTATCCAGCCCACTACGCAGGATCACCACGAACGTCAGACTACACGAGAGCGATGCCAGCGAAGGAccc  >  2:4239‑M2/4‑249 (MQ=255)                                                                                                                                                                                                             acatCTTCCCTTAGCGAAAGGCCCGGTACATAGCCCGGGCAACAGGAATAAATTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCAACACGAACGTAAGACTGAACGAGAGCGAAGCCAGCGCacaaccc  >  2:456710/1‑240 (MQ=255)                                                                                                                                                                                                             acatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTAccc  >  2:455641/1‑246 (MQ=255)                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      CGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTACCC  >  NC_000913/3370833‑3371282

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