Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,158,586	A→G	S33G (AGC→GGC)	oxyR →	oxidative and nitrosative stress transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,158,586	0	A	G	100.0%	232.5 / NA	64	S33G (AGC→GGC)	oxyR	oxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/31);  total (33/31)

CCTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTC  >  NC_000913/4158451‑4158720                                                                                                                                        |                                                                                                                                       ccTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCg                                                                                                                                     >  1:299901/1‑139 (MQ=255) ccTATCGCCATGAACTATCGTGGCGATCGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCg                                                                                                                                     >  1:851484/1‑139 (MQ=255)    aTCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGc                                                                                                                                  >  2:1247052/1‑139 (MQ=255)    aTCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGc                                                                                                                                  >  2:1206602/1‑139 (MQ=255)      cGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCaa                                                                                                                                >  2:500785/1‑139 (MQ=255)      cGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCaa                                                                                                                                >  1:254033/1‑139 (MQ=255)      cGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCaa                                                                                                                                >  2:667610/1‑139 (MQ=255)           tGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCg                                                                                                                           >  2:173100/1‑139 (MQ=255)              aCTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTaa                                                                                                                        <  2:988560/139‑1 (MQ=255)                     tGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGaa                                                                                                                 >  2:1053762/1‑139 (MQ=255)                     tGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGaa                                                                                                                 >  1:1046552/1‑139 (MQ=255)                       gCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTaa                                                                                                                        <  2:883458/130‑1 (MQ=255)                       gCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTaa                                                                                                                        >  1:883458/1‑130 (MQ=255)                       gCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGa                                                                                                               >  1:1286393/1‑139 (MQ=255)                           tggtggATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTGCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAg                                                                                                           <  1:847046/139‑1 (MQ=255)                            ggaggaTGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGc                                                                                                          >  1:359052/1‑139 (MQ=255)                            ggaggaTGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGc                                                                                                          >  2:1307611/1‑139 (MQ=255)                                  tggatAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCg                                                                                                    >  1:77098/1‑139 (MQ=255)                                   ggatAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGt                                                                                                   <  1:173100/139‑1 (MQ=255)                                   ggatAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGt                                                                                                   <  1:622602/139‑1 (MQ=255)                                     atAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGa                                                                                                 >  2:1045946/1‑139 (MQ=255)                                           agtatTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCCTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGGTGGAAGATGAGCTGGGCGTGATGTTGc                                                                                           <  1:696412/137‑1 (MQ=255)                                           aaTATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGc                                                                                           <  2:713254/139‑1 (MQ=255)                                           aaTATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGc                                                                                           <  2:299901/139‑1 (MQ=255)                                           aaTATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGc                                                                                           <  1:667610/139‑1 (MQ=255)                                                tCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAg                                                                                      <  2:133443/139‑1 (MQ=255)                                                tCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAg                                                                                      <  2:1286393/139‑1 (MQ=255)                                                  gTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGc                                                                                                                      >  1:971475/1‑105 (MQ=255)                                                  gTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGc                                                                                                                      <  2:971475/105‑1 (MQ=255)                                                   tGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCgg                                                                                   <  1:990621/139‑1 (MQ=255)                                                      tCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGAcc                                                                                <  2:359052/139‑1 (MQ=255)                                                        ttGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAg                                                                              <  2:254033/139‑1 (MQ=255)                                                          gAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGcc                                                                            >  2:780321/1‑139 (MQ=255)                                                                cTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAg                                                                      >  2:152303/1‑139 (MQ=255)                                                                  ggtggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGAc                                                                                 >  2:652954/1‑126 (MQ=255)                                                                  ggtggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGAc                                                                                 <  1:652954/126‑1 (MQ=255)                                                                  ggtggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAgtg                                                                    >  2:512173/1‑139 (MQ=255)                                                                   gtggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAgtgt                                                                   >  2:132813/1‑139 (MQ=255)                                                                    tggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCa                                                                               <  1:1299316/126‑1 (MQ=255)                                                                    tggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCa                                                                               >  2:1299316/1‑126 (MQ=255)                                                                    tggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTaaa                                                                       <  2:109772/134‑1 (MQ=255)                                                                    tggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTaaa                                                                       >  1:109772/1‑134 (MQ=255)                                                                       cATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGAAGCGGACCAGCCGTAAAg                                                                      >  1:848693/1‑132 (MQ=255)                                                                       cATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGACAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAg                                                                      <  2:848693/132‑1 (MQ=255)                                                                          tGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCAc                                                            <  1:152303/139‑1 (MQ=255)                                                                            gCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCg                                                                                                                           >  1:64753/1‑74 (MQ=255)                                                                            gCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCg                                                                                                                           <  2:64753/74‑1 (MQ=255)                                                                                     cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCt                                                                                                         >  2:403773/1‑83 (MQ=255)                                                                                     cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCt                                                                                                         <  1:403773/83‑1 (MQ=255)                                                                                     cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTa                                                                         >  1:17595/1‑115 (MQ=255)                                                                                     cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTa                                                                         <  2:17595/115‑1 (MQ=255)                                                                                      gCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGt                                                                                                   <  1:947067/88‑1 (MQ=255)                                                                                      gCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGt                                                                                                   >  2:947067/1‑88 (MQ=255)                                                                                       ccATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAg                                                                                      <  2:25861/100‑1 (MQ=255)                                                                                       ccATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAg                                                                                      >  1:25861/1‑100 (MQ=255)                                                                                          ttttCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAAtgctg                                            >  2:887689/1‑139 (MQ=255)                                                                                              cGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCtgg                                        <  1:887689/139‑1 (MQ=255)                                                                                                          gATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGgcgc                            <  1:132813/139‑1 (MQ=255)                                                                                                                       ttAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCg               >  2:386807/1‑139 (MQ=255)                                                                                                                         agccagccGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTg             <  2:1046552/139‑1 (MQ=255)                                                                                                                            cagccGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAgg          <  1:1053762/139‑1 (MQ=255)                                                                                                                               ccGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGa       >  2:129132/1‑139 (MQ=255)                                                                                                                                 gACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGaaa     <  1:1206602/139‑1 (MQ=255)                                                                                                                                    gCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTc  <  2:77098/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       CCTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTC  >  NC_000913/4158451‑4158720

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