Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	257,908	Δ776 bp		insB9–[crl]	insB9, insA9, [crl]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	257908–258674	258683	10–776	40 [0]	[1] 41	insB9–[crl]	insB9,insA9,[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	0 (0.000)	40 (1.020)	38/480	0.2	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 family transposase B
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

CGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCAC  >  NC_000913/257675‑258141                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           cGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCCCAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgat                                                                                                                                                                                                                          <  1:17637‑M1/251‑19 (MQ=255)     aaaTTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctt                                                                                                                                                                                                                      <  1:131663‑M1/251‑23 (MQ=255)                             ccGGGGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTGGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTTCCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattttctggctgtatgcgt                                                                                                                                                                                              <  2:51832‑M1/251‑47 (MQ=255)                                   tATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgt                                                                                                                                                                                        >  1:396679‑M1/1‑199 (MQ=255)                                    aTCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGCTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtg                                                                                                                                                                                       >  1:367786‑M1/1‑198 (MQ=255)                                             ttGACAAAGGTCTTCAGTAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCTATGACGTTACCGAGTGGACACCCGAATAGCAGATTTATCAAAAAATTTACCGCGCTAGGCCGGTATATTCGTGATGGTAagtgcagagagaatcggttcttttgcgtttgtctggttgtatgcgtcatcgtgaaacctgc                                                                                                                                                                               <  2:425390‑M1/250‑62 (MQ=255)                                             ttGACAAAGGTCTTCAGGAAATCATCGACTGGAGCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataagcgattctttttcgatggtcttgctgtatgcttcaacgtgaaaccggc                                                                                                                                                                               <  2:295019‑M1/250‑62 (MQ=255)                                             ttGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGATTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAAATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggc                                                                                                                                                                               >  2:413264‑M1/1‑189 (MQ=255)                                                aCAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAAGTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCTCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgagtctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccg                                                                                                                                                                           <  2:175877‑M1/251‑66 (MQ=255)                                                                aaTCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattct                                                                                                                                                           >  1:312325‑M1/1‑170 (MQ=255)                                                                        aCTGGATCGAAAAACGCTTGTGTTAAAAATTTTCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACATTACCGAGAGGACACCCGAAGAGCAGATTGATCAAAAAATATACCACACTAGGCCCGTATATTCGAGAAGATAagtgcaaagataatagattctttttcgattgtctggctgtatgcgtcaacgtgaacccggcaccggaagtgcgtgaattctggggctgg                                                                                                                                                   >  2:84100‑M1/1‑162 (MQ=255)                                                                               cGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaacggaagtgcgtgaattctggggctggtggatgg                                                                                                                                            >  2:335444‑M1/1‑155 (MQ=255)                                                                                gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatgga                                                                                                                                           <  1:257960‑M1/251‑98 (MQ=255)                                                                                   aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagct                                                                                                                                        >  1:7358‑M1/1‑151 (MQ=255)                                                                                       gCTTGTGTTAACAATTTCCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAACAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggc                                                                                                                                     >  2:91286‑M1/1‑147 (MQ=255)                                                                                          tgtgTTAAAAATTTGCTAAATGTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATATCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAATAAGTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgtctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggatctggaagc                                                                                                                                  <  2:168205‑M1/250‑107 (MQ=255)                                                                                                aaaaaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCTCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgactccggcaccggaagtgcgtgaatgctggggctggtggatggagctggaagcgcagga                                                                                                                            <  1:220832‑M1/250‑113 (MQ=255)                                                                                                aaaaaTTTGCTAAATTTTGCCAATTTGGTAAAACACCTGCATCACCACCGGCGATAGCAATAACGTTACCGAGTGGACACCCGACAAGCAGATTGATCAAAAAATATACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgcttctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgcattctgggactggtggatgaagctggaagcgcagga                                                                                                                            >  2:451499‑M1/1‑138 (MQ=255)                                                                                                  aaaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCACCACAACAGCAGATCGCAATGACGTTGCCCAGCGGACACCCGAAGCCCAGATTGCTCAAAAAATTTACCGCACTAGGCCCGTCTATTCGCGAAGGAAagtgcaaagataatcgattctttttcgattgtctggatgtatgcgccaacgtgaagccggcaccggaagtgcgtggattctggggctgctggatggagctggaagcgcaggaatc                                                                                                                         >  1:30315‑M1/1‑136 (MQ=255)                                                                                                    aTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatccc                                                                                                                       <  1:315423‑M1/251‑118 (MQ=255)                                                                                                           aaaTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccgggagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtttt                                                                                                                  <  1:281634‑M1/249‑123 (MQ=255)                                                                                                           aaaTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttac                                                                                                                <  1:48454‑M1/251‑125 (MQ=255)                                                                                                              ttttGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacct                                                                                                              >  1:324282‑M1/1‑124 (MQ=255)                                                                                                               tttGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctac                                                                                                            >  1:82781‑M1/1‑123 (MQ=255)                                                                                                                                aCAGTTGCATCCCAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGAGACGGTAagtgcaaagataa                                                                                                                                                                                                                               >  2:290836‑M1/1‑106 (MQ=255)                                                                                                                                aCAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                                                                                                                                                                                                               <  1:290836‑M1/119‑14 (MQ=255)                                                                                                                                aCAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATGCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggttgatggagctggaagcgcagtaatcccgttttacctacagttaccagtttggtct                                                                                           <  2:312861‑M1/251‑146 (MQ=255)                                                                                                                                 cAGTTGCATCACAACAGGACATAGCAATGACGTTCCAGAGAGGCCACACGAAGCGCAGATTGATCAAAAACTTTACCGCACTAGGCCCGTATATTCGTGAAGGCAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgagaaaccggcaccggaagtgcgtgaattctggggctggtggctggagctggaagcgcaggaatcccgttttacctacagttacccgtttggtctg                                                                                          >  2:210855‑M1/1‑105 (MQ=255)                                                                                                                                           acaacaGGAGAGAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGAACAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcgtttgcgattgtctggctgtaggcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaagcccggtgtgcgtacagttaccagtttggtctgttcgataaag                                                                                <  1:411966‑M1/251‑157 (MQ=255)                                                                                                                                                  gagaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcg                                                                          >  1:222438‑M1/1‑88 (MQ=255)                                                                                                                                                    gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgact                                                                       >  1:258523‑M1/1‑86 (MQ=255)                                                                                                                                                          aaTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattatttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggacg                                                                  >  2:360882‑M1/1‑80 (MQ=255)                                                                                                                                                                            acacCCGATGAGCAGAGTGATCAAAAAATTGATCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataagcgattctttttcgattgtcttgctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtttgaccgacagttaccagtttggtcttgtcgataaagcaggctattggaagagtgtgccggtaaaagaca                                               <  2:294217‑M1/251‑190 (MQ=255)                                                                                                                                                                             cacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATGCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggccccggaagtgcgtgaattctggggctggtggatggagctggaatcgcaggaatcccgttgtacctacagttagcagtttgggcggttcgatagagcatgcgactggaagagtgttccggtaaaagacac                                              <  1:149098‑M1/251‑191 (MQ=255)                                                                                                                                                                                  gAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactga                                           >  1:107985‑M1/1‑56 (MQ=255)                                                                                                                                                                                                                  tAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtg           >  1:448428‑M1/1‑24 (MQ=255)                                                                                                                                                                                                                  tAGGCCCGTATATTAGTGCAGGTAtgtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttaactacagttaccagtttggtctgttcgctaaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaccgactggagcacaccctgcgtg           >  2:263019‑M1/1‑24 (MQ=255)                                                                                                                                                                                                                   aGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagt        >  1:3824‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                       ccGTATATTCGTGAAGGTAagtgcaacgataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgaaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgcacgactggagcacaccctgcgtgagtttc     >  2:409392‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                           atatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttcacg  >  1:354978‑M1/1‑15 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           CGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCAC  >  NC_000913/257675‑258141

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