Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,131,341	Δ128 bp	IS5‑mediated	yghQ ←	putative transport protein YghQ

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3130367–3131345	3131468	124–1102	9 [8]	[1] 41	[insH9]–[yghQ]	[insH9],[yghQ]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2288919 =	NA (NA)	40 (0.840)	39/484	0.3	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	3131469 =	0 (0.000)					coding (940/1068 nt)	yghQ	putative transport protein YghQ

TAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAA  >  NC_000913/3131233‑3131705                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             aaaaacgataaccaacaggcgatgttattcacgctggccaacctggttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccgtaaatggcgaagaaacggtctaaataggctgattcaaggcagttacgggagaaaaaagcggctcaaacaggaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTc                                                                                                                                                                                                                                <  2:75740‑M2/15‑1 (MQ=255)          aaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGgtggtg                                                                                                                                                                                                                        <  2:104213‑M2/23‑1 (MQ=255)                 gggcgatgttattcacgcgggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaattgccccgcttattcgcaccttccCTAACAACCTCGGTCTGGTGGTGCGCGGATc                                                                                                                                                                                                                <  2:425503‑M2/31‑1 (MQ=255)                              cacgctggccaacctgattcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggcgcaaacatgaagaaatgaaatgactgagtcagcagagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCAcc                                                                                                                                                                                                   >  1:112221‑M2/208‑251 (MQ=255)                                                 cgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGAAAAAACTTTTAc                                                                                                                                                                                >  1:388574‑M2/189‑251 (MQ=255)                                                      ggaccaaattatacgtcagggggagagttctcactaaaaactggggataacgccttaaatggccaagaaacggtctaaataggctgattcaaggcatttacgggggaaaaaagcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGc                                                                                                                                                                           <  2:340161‑M2/68‑1 (MQ=255)                                                        accaaatgatacgtcagtgggagagatctcaataaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcacaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCa                                                                                                                                                                          >  1:201078‑M2/182‑250 (MQ=255)                                                         ccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCaaa                                                                                                                                                                        <  2:422799‑M2/71‑1 (MQ=255)                                                          caaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAAc                                                                                                                                                                       >  2:170197‑M2/180‑251 (MQ=255)                                                             atgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCg                                                                                                                                                                    <  1:150674‑M2/75‑1 (MQ=255)                                                               gatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagacatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATATAAACGCATCACCTCTAGGTAAAAACTTTTACCCAGcaaaccagcc                                                                                                                                                                  >  2:499629‑M2/175‑247 (MQ=255)                                                                  acgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGc                                                                                                                                                               <  1:204064‑M2/80‑1 (MQ=255)                                                                      cagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTc                                                                                                                                                           >  2:50995‑M2/168‑251 (MQ=255)                                                                       agtggtagagctctcactaaaaactggggataacgccttcaatggggaagaaacggtctaaataggctgattcgaggcatttacgggagaaaaaatcggcccacacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGTAGGTGCGCGTCTCTAAACGCATCACCACTGGGTAAAAACTTTTACCCAGCCACCCAGCTGGCGTTcc                                                                                                                                                          >  2:322018‑M2/167‑251 (MQ=255)                                                                          gggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGCTCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcggc                                                                                                                                                       >  1:423015‑M2/164‑251 (MQ=255)                                                                            gagagatctcactaaaaactggggataactccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcggcgg                                                                                                                                                     <  2:448764‑M2/90‑1 (MQ=255)                                                                               agatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGAGGTGCGCGGATCTAAACGCATCACCTCTGGGAAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGt                                                                                                                                                  >  1:1273‑M2/159‑251 (MQ=255)                                                                                gatctcactaaaacctggggataacgccttaaatgccgacgaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagacatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAAAGCATCACCTCTGGGCAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTc                                                                                                                                                 >  2:198209‑M2/158‑251 (MQ=255)                                                                                            aactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaagacctttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTTCGCGGACCTAAACGCATCACCTCTGGGGAAAAACTTTTACCCAGCAAACCCGCAGGCGTTCCGGCGGCGTCGCAGCATGTCAt                                                                                                                                     >  1:434461‑M2/146‑251 (MQ=255)                                                                                                 gggataacgccttaaatggcgaagaaacggtctcactaggctggttcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtacgccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCg                                                                                                                                 >  2:288047‑M2/141‑250 (MQ=255)                                                                                                   gataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAtt                                                                                                                              >  1:331116‑M2/139‑251 (MQ=255)                                                                                                               aatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCcg                                                                                                                  <  1:277759‑M2/125‑1 (MQ=255)                                                                                                                      aagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCggcggc                                                                                                             >  2:331470‑M2/120‑249 (MQ=255)                                                                                                                         aaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCgg                                                                                                        >  1:334925‑M2/117‑251 (MQ=255)                                                                                                                           acggactaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCCTGGCGATTTTAAATAACCCGGCGGCAGcgggg                                                                                                      >  2:403211‑M2/115‑250 (MQ=255)                                                                                                                               tctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATc                                                                                                                                                                                                      >  2:175536‑M2/111‑151 (MQ=255)                                                                                                                               tctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATc                                                                                                                                                                                                      <  1:175536‑M2/41‑1 (MQ=255)                                                                                                                                       ggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATc                                                                                          >  1:228298‑M2/103‑251 (MQ=255)                                                                                                                                        gctgattcaaggcatttacgggagaaaaaatcggcgcaaacatgaagaaacgaaatgactgagtcagccgagaagaagttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCGCCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATcc                                                                                         <  2:424843‑M2/150‑1 (MQ=255)                                                                                                                                                             gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACg                                                                    <  1:113037‑M2/171‑1 (MQ=255)                                                                                                                                                                   aaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAt                                                                                                                               >  2:337135‑M2/75‑186 (MQ=255)                                                                                                                                                                   aaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAt                                                                                                                               <  1:337135‑M2/112‑1 (MQ=255)                                                                                                                                                                   aaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGAATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCAACCAACACTGTGCTGCACGAGTTAc                                                              >  1:327006‑M2/75‑251 (MQ=255)                                                                                                                                                                                     agaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGt                                            <  1:253721‑M2/195‑1 (MQ=255)                                                                                                                                                                                             aatgactgggtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATg                                    <  2:228298‑M2/203‑1 (MQ=255)                                                                                                                                                                                                 actgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAAcc                                                                                                                    >  2:90973‑M2/45‑167 (MQ=255)                                                                                                                                                                                                 actgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAAcc                                                                                                                    <  1:90973‑M2/123‑1 (MQ=255)                                                                                                                                                                                                      gtcagccgagaagaatttccccgcttattcgcaccttccCGAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGAc                           <  2:157181‑M2/212‑1 (MQ=255)                                                                                                                                                                                                                       tccccgcttattcgcaccttcccTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCACACCCGCTGGCGTTCCGCCGACGTCGAAGACTGTCATGGCAATTTTCACTAACCCGGCGGCCGCGGGTCCTACCACGATCCCCACCACAACTGTCCTACACGAGTTACGAGCCGACCAGATGGGATGGGCAATGTTGGTTGCCCGGGCAAAACTcccccc          >  2:421904/22‑247 (MQ=255)                                                                                                                                                                                                                             cttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGCTTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCAACCAACACTGTGCTTCACGAGTTACGCGCCGACAAGATGGAGGGGGCAATATTGGTTGACCAGACAAAGCTCCACGCGCCttt    >  2:16041‑M2/17‑251 (MQ=255)                                                                                                                                                                                                                                attcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTaa  <  1:184075‑M2/237‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             TAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAA  >  NC_000913/3131233‑3131705

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