Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,223,638	T→C	intergenic (‑10/‑190)	iclR ← / → metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,223,638	0	T	C	100.0%	106.3 / NA	35	intergenic (‑10/‑190)	iclR/metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (2/33);  total (2/33)

TCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGT  >  NC_000913/4223396‑4223880                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   tCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCt                                                                                                                                                                                                                                             <  1:164821/250‑1 (MQ=255)  ggCCAACCTGATGCCCTAACCCCTGCTGTTGCATCGTGGTTAGCAGGCGGTTGGTCGTTGAATTGGGTAACCCGGCTTGTGGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGGTAAAGACTGAACCTGTCCAGTCGCTGGGGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGGCTCCTTTTTCTGt                                                                                                                                                                                                                                           <  2:546787/250‑1 (MQ=255)    ccAACCTGCCGAACGAAACCCTGCTGTTGCATCGTGTTTAGCAGGCGTGGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGTGAATGGGTGCGTCCATGACAGTCTCCTTTTTCTGTAt                                                                                                                                                                                                                                         <  2:562512/251‑1 (MQ=255)    ccAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAt                                                                                                                                                                                                                                         <  1:571560/251‑1 (MQ=255)        ccTGACTCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGTAATTGTGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCTGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTg                                                                                                                                                                                                                                     <  1:101609/251‑1 (MQ=255)             cGCACTAAACCCGGCTTTGGCCTCGTTGTTAGCAGCCGTTGGGCCGGGGAATTTGGTAATCCGGCTTGTTGCACTAGTTTCGTAAGTGCCCCACTGCCATTGGATTTGGCAATCCACTCCAGTAATTTCAGGCCACGATTTATAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTGTCGCGGTAATGGGTGCGACCATGACTGTCGCCTTTTTTTGTATCGTGGAAAt                                                                                                                                                                                                                                <  2:419517/251‑1 (MQ=255)             cGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGTTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAt                                                                                                                                                                                                                                <  1:313786/251‑1 (MQ=255)                  gAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattt                                                                                                                                                                                                                           <  1:439164/251‑1 (MQ=255)                     aCCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttca                                                                                                                                                                                                                        <  1:633747/251‑1 (MQ=255)                      cccTGCTGTTGTAGCGTGGGTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTCTTTCTGCATCGTGGAAAtcattttcat                                                                                                                                                                                                                       <  1:332163/251‑1 (MQ=255)                       ccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACGCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttcatt                                                                                                                                                                                                                      <  1:75163/251‑1 (MQ=255)                          tctggtGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttcattt                                                                                                                                                                                                                     <  2:533226/244‑1 (MQ=255)                          gctgTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCTTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAAGTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCAttttt                                                                                                                                                                                                                   <  2:423315/251‑1 (MQ=255)                                gCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGtt                                                                                                                                                                                                             <  1:468692/251‑1 (MQ=255)                                  atggTGGTTAGCAGGCGGTGGGTCGTGGAATTGGTGAACCCGGCTTGTTTCGCCAGTTCCGGGAGTGCCACACTGCCATTGGATTTGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTTAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAGCATTTTCATTTCTATTGTTAg                                                                                                                                                                                                           <  2:49072/248‑1 (MQ=255)                                           agcagtgggcgtgccgttgAATTGGGTAACCCGTCGTGTTGCGCCAGTTCCCTGAGTGCCCCTCTGCCATTGGAGTCGGCGATCCACTCCAGTAAGTGCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTTGTGCGGTGGCAACGGCGGGTTTTCGGCCGCGTTTCGCGGGAATGGGTGCGACCTTGACAGTCTCCTTTTGCTGGATCGTGGAAAGCATTTTCAGTTTTATTGTTAGCTAATGCa                                                                                                                                                                                                   <  2:646590/232‑1 (MQ=255)                                                                 tGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGAt                                                                                                                                                                             <  1:9658/250‑1 (MQ=255)                                                                  gggTAACCCGGCTTGGTGCCCCAGTGCCGTGTGTGCCACACTGCCATTGTATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGGTAAAGACTGAACCGGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTTTCCTTTTTCTGTATCGTGGAAATCATTTGCATTTGTATGGTTATCTAATGCAATAGTTACTGAACTGATCCGATg                                                                                                                                                                            <  2:396270/250‑1 (MQ=255)                                                                      aaCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTa                                                                                                                                                                       <  1:184485/251‑1 (MQ=255)                                                                          cGGCTTGTTGCGCCAGGTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCATTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGt                                                                                                                                                                   <  2:533182/251‑1 (MQ=255)                                                                                                         acTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATGGGTGCGGCCATGACAATCTCCTTTTTTTTTATCGTGGGAAACATTTTTATTTTTTTTGTTTGCTAATGCAATAGTTACTTAACTGATCCGATGGGGTTATGTTTGACAAATCTCATGTTGAGTGGTGGTCGCt                                                                                                                                    >  1:460525/1‑251 (MQ=255)                                                                                                             ccATTGGATTCGGAAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCGGTCCCGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGTTAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTa                                                                                                                                <  2:369163/251‑1 (MQ=255)                                                                                                                            aTCCACTCCAGTAATGGCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTTGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTa                                                                                                                 <  2:171756/251‑1 (MQ=255)                                                                                                                                   ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAg                                                                                                           <  1:191856/250‑1 (MQ=255)                                                                                                                                        aaTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCGGTCGCGTGTGGGGGGGCAACGGCGGGTTTGCTGCCGCTTGTAGCGGGAAGGGGTTCGACCATGACAGTCGCCTGTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAAGGTGACTGAACTGATCCGATGAGTTAATGTTGAATAAATCTCATTTTGCGTGGTGGTCGGTTGTACCGCAGATGCGTTTTTGCGAGTATga                                                                                                      <  2:68677/250‑2 (MQ=255)                                                                                                                                             cAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgt                                                                                                 <  1:464587/250‑1 (MQ=255)                                                                                                                                               ggCCACGCGTTAAAGTCTGAACCTGTCCAGGCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTGTTACCACAGATGCTTTTATGCCAGTATGGTTTGTTGa                                                                                              <  2:250788/251‑1 (MQ=255)                                                                                                                                                   aCGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAttt                                                                                          <  1:59164/251‑1 (MQ=255)                                                                                                                                                          aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAg                                                                                                                                                                          <  2:54026/164‑1 (MQ=255)                                                                                                                                                          aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGGGCGGCCATGACAGTCTCCTTTTTTTGTTTCGTGGGAATCATTTTTATTTTTTTTTTTAGCTAATGCAATAGTTACTGAAATGATCCGCTGAg                                                                                                                                                                          >  1:54026/1‑164 (MQ=255)                                                                                                                                                                    ccTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTg                                                                         <  1:467887/251‑1 (MQ=255)                                                                                                                                                                      ggttccgttgggggttcgtgggCAACGGCGGGTTTTCTGTCGCGGTTCGCGTGAATGGGTTCGACCAGGACAGTCTTCTTTTGCTGTATCGTGGAAATCATTTTCATGTTTAGTTTTATCTAATGCAATAGGTACTGAACTGATCCGATGAGTTAATGTTGAACAAAGCTCAGGTTGTGTGGTGGTCGCTTGTACCACAGATGCGTTGATGCCAGTATGGTTTGTTTAATTTTTATGAAATCGGGGTTGa                                                                        <  2:611403/230‑1 (MQ=255)                                                                                                                                                                                     gCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGGTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAg                                                        <  1:384557/251‑1 (MQ=255)                                                                                                                                                                                                                    tCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGGGGAACAACtgcgtgc                          <  1:537300/250‑1 (MQ=255)                                                                                                                                                                                                                                           aCAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGt  <  2:37164/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   TCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGT  >  NC_000913/4223396‑4223880

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