Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,637,577	T→C	T294T (ACT→ACC)	creC →	sensory histidine kinase CreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,637,577	0	T	C	100.0%	108.1 / NA	35	T294T (ACT→ACC)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/22);  total (13/22)

CGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATT  >  NC_000913/4637338‑4637818                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  cGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCt                                                                                                                                                                                                                                        >  1:135212/1‑251 (MQ=255)                gTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGcgcg                                                                                                                                                                                                                         <  2:163228/251‑1 (MQ=255)                                              gaactcGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTtacta                                                                                                                                                                                             <  1:140776/246‑1 (MQ=255)                                               aactcGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGc                                                                                                                                                                                          >  1:137208/3‑251 (MQ=255)                                                   cggcagtagCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCGGGAAGGGAAAAACAATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCGAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCGCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAgg                                                                                                                                                                                      <  2:2193/247‑1 (MQ=255)                                                              aGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCGGGAAGGGAAAAACTATATTGAGCAGTATGGTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGGTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTTCTCCACGCAAGGCTGgag                                                                                                                                                                           <  2:120553/251‑1 (MQ=255)                                                                            gcgcATGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTAGATTTAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGTCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCGCGTTTTACCGACAACATTCTGTCGCAAAATGCGCGTAGGCATGTATTGGTAGAAACGTTACTACGCCAGGcgtggcgggttaatcgttaggaag                                                                                                                                                              <  2:52499/250‑24 (MQ=255)                                                                             cgcAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATGCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTc                                                                                                                                                            <  1:39214/251‑1 (MQ=255)                                                                                aGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGtt                                                                                                                                                         >  1:107596/1‑251 (MQ=255)                                                                                  gcgcTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCt                                                                                                                                                       <  2:219291/251‑1 (MQ=255)                                                                                    gcTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGa                                                                                                                                                     <  1:245048/251‑1 (MQ=255)                                                                                       ggAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCAGTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGAc                                                                                                                                                    <  2:82470/249‑1 (MQ=255)                                                                                              aTGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGttgttg                                                                                                                                            >  1:76439/1‑250 (MQ=255)                                                                                                        aGCTGGAAGGGAAACACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTggcggc                                                                                                                                  >  1:49473/1‑250 (MQ=255)                                                                                                           tGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTCAAAAGCCCACTGGCGGCGAATCGTGGCGCGGCGGAAATTTTACGCGAAAGGCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGAGCGTATGAAGGCATTGGTAGAAACGTTACTACGCCACGAAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTggcggccttt                                                                                                                              >  1:269104/1‑247 (MQ=255)                                                                                                            ggAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCatt                                                                                                                               <  1:254356/249‑1 (MQ=255)                                                                                                              aaGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCattattt                                                                                                                           <  1:174931/251‑2 (MQ=255)                                                                                                               aGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCattatttc                                                                                                                          <  2:208605/251‑3 (MQ=255)                                                                                                                              tGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACGACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTTTTGGTGATGTGACGGCATTAGTTCGCCGCGTCAGCGaa                                                                                                            <  2:122595/250‑1 (MQ=255)                                                                                                                                aGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGGAATTTTACGCGAAAGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAACATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGGCGTTCTGACTGTTGTTGATGTGGCGGGattttttgccccgggcagcgagggg                                                                                                         >  2:159855/1‑229 (MQ=255)                                                                                                                                  cAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAgcgcg                                                                                                       <  1:52764/251‑1 (MQ=255)                                                                                                                                              gcgtttacccatgagagaaaaaGCCCACTGGCGTCGATTCGTGGCGCGTCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGAGGTTGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGtt                                                                                           <  2:267396/234‑1 (MQ=255)                                                                                                                                               cgttAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTg                                                                                          <  1:224095/249‑1 (MQ=255)                                                                                                                                                      ccATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACGGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATGATTTCGCCGCGGCAGCGAAGCGCGCAACGTGCAGGTGGCAGaaa                                                                                   <  2:137208/250‑1 (MQ=255)                                                                                                                                                                    gCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCActcta                                                                     <  1:264671/251‑2 (MQ=255)                                                                                                                                                                         cTGGCAGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACCTTCTGCCGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTTGCAGAAAAAAACATCCCTCTACATGt                                                                >  2:249301/1‑251 (MQ=255)                                                                                                                                                                            gcggcgATTCGTGGCGCGGCGGAACTTTTACGCGCAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGGCAGCGGAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACCTGTTat                                                             >  2:232578/1‑250 (MQ=255)                                                                                                                                                                              ggcgATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGAAGAAAAAAACATCACTCTACATGTTatgc                                                           >  1:57102/1‑248 (MQ=255)                                                                                                                                                                                cgATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTTTTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTTACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGccc                                                         <  1:250708/251‑1 (MQ=255)                                                                                                                                                                                  aTTCGTGGCGCGGCGGAAATTTTACGCGACGGTCCGCCACCGGAAGTGGCGGCTCGTTTTCCCGACAACCTTCTGACGCAAAATGCGCGTATGCAGGGATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGGTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTatggccac                                                       >  1:107881/1‑244 (MQ=255)                                                                                                                                                                                  aTTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCTGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCCTTGGTAGAACCGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTTGCGGCATTCTTTCGCCGCGTCAGCGGAGCGCGCACCGTGCAGTTGGCAGAAAACAACATCACTCTACATTTTATGCCCAc                                                       >  1:238148/1‑251 (MQ=255)                                                                                                                                                                                                  aaaTTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTtgct                                       <  1:60691/251‑1 (MQ=255)                                                                                                                                                                                                       ttACGCGAAGGGCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGctgctg                                   <  2:107596/250‑1 (MQ=255)                                                                                                                                                                                                                           ggAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAAACGGCGTTACTGGAGCAgg              >  1:85888/1‑251 (MQ=255)                                                                                                                                                                                                                                        cGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTAGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAAtt  <  1:239517/250‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  CGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATT  >  NC_000913/4637338‑4637818

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