Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,165,573	(C)6→7	intergenic (+52/‑95)	yegP → / → yegQ	DUF1508 domain‑containing protein YegP/putative peptidase YegQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,165,567	1	.	C	100.0%	189.2 / NA	51	intergenic (+46/‑101)	yegP/yegQ	DUF1508 domain‑containing protein YegP/putative peptidase YegQ
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (21/30);  total (21/30)

CAGCCAGATGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAG‑C‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAA  >  NC_000913/2165428‑2165701                                                                                                                                             |                                                                                                                                         cagccagATGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCagag                                                                                                                                           >  2:479498/1‑140 (MQ=255)        aTGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAG‑c‑cccc‑cc                                                                                                                                 >  1:850543/1‑139 (MQ=255)         tGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAG‑c‑cccc‑cca                                                                                                                                <  1:535950/140‑3 (MQ=255)          gTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAG‑c‑cccc‑ccac                                                                                                                               <  1:944912/140‑4 (MQ=255)             cGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAgc                                                                                                                             <  2:269597/140‑1 (MQ=255)                  ccGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTAc                                                                                                                        >  2:330829/1‑140 (MQ=255)                  ccGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTAc                                                                                                                        >  2:67664/1‑140 (MQ=255)                   cgcgCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACa                                                                                                                       >  1:1199505/1‑140 (MQ=255)                    gcgcAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACaa                                                                                                                      >  1:276805/1‑140 (MQ=255)                     cgcAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAAt                                                                                                                     >  2:46720/1‑140 (MQ=255)                        aaTCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATgcc                                                                                                                  <  1:399049/140‑1 (MQ=255)                        aaTCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATgcc                                                                                                                  <  1:247270/140‑1 (MQ=255)                        aaTCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATgcc                                                                                                                  <  1:257036/140‑1 (MQ=255)                          tctcGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATgcccg                                                                                                                <  2:1208630/140‑1 (MQ=255)                                     cGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTg                                                                                                     >  1:169128/1‑140 (MQ=255)                                      ggAATTGCCTCCGTTATAGCCAATGGCACAACCCAGACGGTGAAAGACAATACGTATTTACGATGCCGGGCGTGGGATTTTGACGTCCGGTGCTTTAACACAG‑CGCCCC‑CACAGCGCTACAATCCCCGCCCTTAAATTgg                                                                                                    >  2:1497513/1‑140 (MQ=255)                                        aaTTGCCTCCGTTAAAGCCAATGGCTCAAGCCAGACGGTGAAAGACAATACGGATTTACGGGGCGGGGCGTGGGTTGTTGGGGTCCGGCGCTTTACCAGAG‑C‑CCCCGCGCAGCGCTCCAATGCCGGCCCTTAGAGTgggg                                                                                                  <  2:23517/140‑1 (MQ=255)                                          ttGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGc                                                                                                <  2:565210/140‑1 (MQ=255)                                          ttGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGc                                                                                                <  2:676432/140‑1 (MQ=255)                                          ttGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGc                                                                                                <  1:1336773/140‑1 (MQ=255)                                          ttGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGc                                                                                                <  1:1311012/140‑1 (MQ=255)                                           tGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCa                                                                                               <  2:1351863/140‑1 (MQ=255)                                           tGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCa                                                                                               <  1:1526189/140‑1 (MQ=255)                                           tGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCa                                                                                               <  1:1420111/140‑1 (MQ=255)                                             ccTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACt                                                                                             >  1:1424936/1‑140 (MQ=255)                                                 cGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTcccc                                                                                         >  2:777519/1‑140 (MQ=255)                                                       aGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCg                                                                                   <  2:309111/140‑1 (MQ=255)                                                           aaTGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtc                                                                               >  2:577315/1‑140 (MQ=255)                                                            aTGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtca                                                                              >  2:443751/1‑140 (MQ=255)                                                                   aaGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTg                                                                       >  2:788330/1‑140 (MQ=255)                                                                     gCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGaa                                                                     >  2:1271418/1‑140 (MQ=255)                                                                      ccAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAg                                                                    <  1:1519562/140‑1 (MQ=255)                                                                      ccAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAg                                                                    <  2:824522/140‑1 (MQ=255)                                                                           cGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGAt                                                               <  1:1429605/140‑1 (MQ=255)                                                                                aaaGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGAc                                                          <  1:330829/140‑1 (MQ=255)                                                                                aaaGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGAc                                                          <  2:1178033/140‑1 (MQ=255)                                                                                aaaGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGAc                                                          <  2:1199505/140‑1 (MQ=255)                                                                                    aCAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGt                                                      >  2:99662/1‑140 (MQ=255)                                                                                       aTACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTCCAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtcat                                                   <  2:967278/140‑1 (MQ=255)                                                                                       aTACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtcat                                                   <  2:850543/140‑1 (MQ=255)                                                                                        tACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtcatc                                                  >  1:460255/1‑140 (MQ=255)                                                                                        tACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtcatc                                                  >  2:787668/1‑140 (MQ=255)                                                                                                 aCGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgaga                                         <  1:950615/140‑1 (MQ=255)                                                                                                 aCGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgaga                                         <  1:51568/140‑1 (MQ=255)                                                                                                  cGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgagag                                        >  2:933238/1‑140 (MQ=255)                                                                                                       ccGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTa                                   >  1:1202738/1‑140 (MQ=255)                                                                                                       ccGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAG‑c‑cccc‑c                                                                                                                                  <  2:1308475/44‑1 (MQ=255)                                                                                                       ccGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAG‑c‑cccc‑c                                                                                                                                  >  1:1308475/1‑44 (MQ=255)                                                                                                              tGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTaa                            >  2:296551/1‑140 (MQ=255)                                                                                                                ggTTGTTGGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGCGAGAATTATGTTTAAAc                          >  1:672528/1‑140 (MQ=255)                                                                                                                     ttgGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGaa                     <  1:1179911/140‑1 (MQ=255)                                                                                                                     ttgGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGaa                     >  2:1328194/1‑140 (MQ=255)                                                                                                                     ttgGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGaa                     <  2:1072612/140‑1 (MQ=255)                                                                                                                      tgGCGTCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAAc                    <  1:788330/140‑1 (MQ=255)                                                                                                                           tCCGGCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCtt               >  1:1499088/1‑140 (MQ=255)                                                                                                                              ggCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTcc            <  1:1271418/140‑1 (MQ=255)                                                                                                                              ggCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTcc            <  1:1256210/140‑1 (MQ=255)                                                                                                                              ggCGCTTTACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTcc            <  1:246394/140‑1 (MQ=255)                                                                                                                                     tACCAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCcggcgg     <  1:443751/140‑1 (MQ=255)                                                                                                                                        cAGAGCC‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGTGGATCTGCCCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGaa  <  1:777519/140‑1 (MQ=255)                                                                                                                                             |                                                                                                                                         CAGCCAGATGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAG‑C‑CCCC‑CACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAA  >  NC_000913/2165428‑2165701

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