Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,587,112	C→T	A461V (GCA→GTA)	narQ →	sensor histidine kinase NarQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,587,112	0	C	T	100.0%	110.9 / NA	36	A461V (GCA→GTA)	narQ	sensor histidine kinase NarQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/22);  total (14/22)

CTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTT  >  NC_000913/2586867‑2587352                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  cttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAATGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGCGCTTTTGACTTCTTTCCGCCTTACCCTTCCGCAGGCTGATCTCCCCTCCGCATTTAGGGACATTCTGGGTTCGTTACAAAATCCAACCATCGCCAAACTGCACCTCCACTGCCGTCTTCCCACCCTGGTACTgg                                                                                                                                                                                                                                              >  2:312728/1‑251 (MQ=255)    acGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGAt                                                                                                                                                                                                                                            <  2:180103/250‑1 (MQ=255)      gTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATgcg                                                                                                                                                                                                                                         >  1:260469/1‑251 (MQ=255)      gTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATgcg                                                                                                                                                                                                                                         <  1:342088/251‑1 (MQ=255)                      aCTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCAtt                                                                                                                                                                                                                         <  1:355886/251‑1 (MQ=255)                                 cGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCGCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAatt                                                                                                                                                                                                              <  2:320901/251‑1 (MQ=255)                                     aCCGGAAGATAACGCCACCTCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAattatt                                                                                                                                                                                                           <  1:174055/250‑1 (MQ=255)                                      ccGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGGGCATTTGTTGCAAATTATTcg                                                                                                                                                                                                         >  2:63343/1‑251 (MQ=255)                                                    cACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGGGCATTTGTTGCAAATTATTCGCGAAGCGGGGCTg                                                                                                                                                                                            >  1:5005/1‑250 (MQ=255)                                                       cGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATg                                                                                                                                                                                        <  1:87138/251‑1 (MQ=255)                                                            aaaGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGTAGTTACGCGATCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCAGTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATg                                                                                                                                                                                   <  2:369204/251‑1 (MQ=255)                                                              aGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGaa                                                                                                                                                                                 <  2:31032/251‑1 (MQ=255)                                                                  tCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGCGCTGTTGACTACTTTCCGCCTGACGCTGCCGCAGGCGGATCTCCCCTCCGCATTGAGGGCACTGCTGGATACGTTACAAAATCCAACACGCGCCAAACTGCCCCTCGACTGCCGTCTGCCAACCCTTGTACTGGATTCACAAATGCCGGTGCATTTGTTGCCAATTATTCGCGACGCGGTGCTTCATGCGATGAAGCac                                                                                                                                                                             >  2:21020/1‑251 (MQ=255)                                                                       gCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAGCTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacg‑ccaa                                                                                                                                                                       <  1:187320/251‑1 (MQ=255)                                                                                          tGAATGATGCTGATCGGCAGTTACGCGAGCTGTTGACTACGTGCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGt                                                                                                                                                     <  2:112197/250‑1 (MQ=255)                                                                                                                cctcgAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACg                                                                                                                              <  2:5005/248‑1 (MQ=255)                                                                                                                                            aCGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCAGGATACGTTACAAAATCAACCCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTAGCGAAGCGGTGCTTCATGCGATGAAGCAAGACCAA‑GCCAGCGGAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGATTATATCCGTGATa                                                                                                   >  2:296170/1‑250 (MQ=255)                                                                                                                                                       ggCGGATTCCCCCACCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACGGACCCGCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAAGGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTa                                                                                       <  2:96864/251‑1 (MQ=255)                                                                                                                                                                 ccccccGCATTGAGGGAAATGCTTGATTGTTTAAAAAATCAAACCAGCGCCAAACTGGCCCTCGACTGGCGTCTGCCAAACCTGGTACTGGATGCTCAAATGCAGGTGCATTTGTTGCGAATTATTCGCGATGCGGTGCTGAATGCGATGAAGCTCG‑CCAACGCCAGCGAAATCGCAGGCAGTTGTGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAAcc                                                                             <  2:211159/247‑1 (MQ=255)                                                                                                                                                                   cTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCg                                                                            >  2:378222/1‑250 (MQ=255)                                                                                                                                                                         aTTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAAc                                                                     >  1:185786/1‑251 (MQ=255)                                                                                                                                                                            gAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCATCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGCTTGGAATCGGTGAACCGAAAGAACCCg                                                                  >  2:257408/1‑251 (MQ=255)                                                                                                                                                                                        ggATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATg                                                      <  1:353376/251‑1 (MQ=255)                                                                                                                                                                                         gatgcGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCTCAAATGCAGGTGCAGTTGTTGCAAATTATGCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTGTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATgg                                                     <  1:30495/247‑1 (MQ=255)                                                                                                                                                                                           tACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTc                                                   >  1:110992/1‑251 (MQ=255)                                                                                                                                                                                                tACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGaa                                               <  1:63343/250‑1 (MQ=255)                                                                                                                                                                                                 aCAAAATCAATACAGCGCCAAACTGACCCTCGACGGCCGTCTGCCACCCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGACAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGGGAACCGAAAGATCCCGAAGGTCATTATGGTCTGAat                                              <  2:295407/250‑1 (MQ=255)                                                                                                                                                                                                    aaaTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCTCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATc                                           >  1:24121/1‑250 (MQ=255)                                                                                                                                                                                                         aaaCCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATgcgc                                     >  1:95564/1‑251 (MQ=255)                                                                                                                                                                                                           aCCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGa                                   <  1:64006/251‑1 (MQ=255)                                                                                                                                                                                                                           cccTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTa                   <  1:126225/251‑1 (MQ=255)                                                                                                                                                                                                                               cGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCTGTTGTTGCAAATTATCCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑ACAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGGGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAgg                 <  2:95564/249‑1 (MQ=255)                                                                                                                                                                                                                               cGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTg               <  1:85619/251‑1 (MQ=255)                                                                                                                                                                                                                                            ccAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACttt  >  1:111776/1‑251 (MQ=255)                                                                                                                                                                                                                                               aCCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTCTTCTCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAAccc                                                                   >  2:203251/1‑183 (MQ=255)                                                                                                                                                                                                                                               aCCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAAccc                                                                   <  1:203251/183‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  CTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACG‑CCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTT  >  NC_000913/2586867‑2587352

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