Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,555,706	(G)6→5	intergenic (‑37/+120)	maeA ← / ← sra	malate dehydrogenase, NAD‑requiring/30S ribosomal subunit protein S22

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,555,701	0	G	.	100.0%	254.4 / NA	59	intergenic (‑32/+125)	maeA/sra	malate dehydrogenase, NAD‑requiring/30S ribosomal subunit protein S22
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (29/30);  total (29/30)

TCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAG  >  NC_000913/1555566‑1555836                                                                                                                                        |                                                                                                                                        tcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTggggg                                                                                                                                     <  2:1204126/140‑1 (MQ=255)       cTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTa                                                                                                                              <  2:795527/140‑1 (MQ=255)        tGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTAt                                                                                                                             <  2:1314779/140‑1 (MQ=255)        tGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTAt                                                                                                                             <  1:996657/140‑1 (MQ=255)          aaGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCg                                                                                                                           >  1:1095569/1‑140 (MQ=255)          aaGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCg                                                                                                                           >  1:166759/1‑140 (MQ=255)                 tGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTa                                                                                                                              >  2:1406083/1‑130 (MQ=255)                 tGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTa                                                                                                                              <  1:1406083/130‑1 (MQ=255)                 tGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGc                                                                                                                    >  1:1388167/1‑140 (MQ=255)                 tGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGc                                                                                                                    >  1:804114/1‑140 (MQ=255)                  gCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCt                                                                                                                   >  1:920069/1‑140 (MQ=255)                   ccTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCtt                                                                                                                  <  1:882121/140‑1 (MQ=255)                    cTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGt                                                                                                                          <  2:399251/131‑1 (MQ=255)                    cTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGt                                                                                                                          >  1:399251/1‑131 (MQ=255)                         ttCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGc                                                                                                            >  1:1369064/1‑140 (MQ=255)                             acaacGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAgg                                                                                                        >  2:287213/1‑140 (MQ=255)                             acaacGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAgg                                                                                                        >  2:740526/1‑140 (MQ=255)                             acaacGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCNTTTATCGTCTTTGCTTTACCGCCAgg                                                                                                        >  2:186172/1‑140 (MQ=255)                              caacGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAggg                                                                                                       >  2:696572/1‑140 (MQ=255)                                 cGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGt                                                                                                    <  2:1369064/140‑1 (MQ=255)                                 cGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACCCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGt                                                                                                    <  1:1417660/140‑1 (MQ=255)                                  ggAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑‑GGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCg                                                                                                  <  1:1293655/140‑1 (MQ=255)                                     aaTTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGc                                                                                                <  1:740526/140‑1 (MQ=255)                                     aaTTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGc                                                                                                <  2:496522/140‑1 (MQ=255)                                             cagTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGta                                                                                        <  1:287213/140‑1 (MQ=255)                                                  cAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAg                                                                                   >  2:789818/1‑140 (MQ=255)                                                     ggCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCAATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAg                                                                                >  1:273346/1‑140 (MQ=255)                                                        cAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATa                                                                             <  1:696572/140‑1 (MQ=255)                                                          gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaa                                                                           >  2:607337/1‑140 (MQ=255)                                                          gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaa                                                                           >  1:166544/1‑140 (MQ=255)                                                          gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaa                                                                           >  2:643782/1‑140 (MQ=255)                                                          gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaa                                                                           >  2:924568/1‑140 (MQ=255)                                                          gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaa                                                                           >  1:188524/1‑140 (MQ=255)                                                          gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaa                                                                           >  1:635454/1‑140 (MQ=255)                                                             tAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaaaac                                                                        <  2:1231730/140‑1 (MQ=255)                                                                  gATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaaaacaaaaa                                                                   <  2:166544/140‑1 (MQ=255)                                                                     ataAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaaaac                                                                        >  2:595521/1‑132 (MQ=255)                                                                     ataAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaaaac                                                                        <  1:595521/132‑1 (MQ=255)                                                                       aaaGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCa                                                              <  1:643782/140‑1 (MQ=255)                                                                           cGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgc                                                          <  2:737535/140‑1 (MQ=255)                                                                           cGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgc                                                          >  2:401652/1‑140 (MQ=255)                                                                           cGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgc                                                          >  1:1147191/1‑140 (MQ=255)                                                                           cGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgc                                                          >  2:205455/1‑140 (MQ=255)                                                                           cGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgc                                                          <  2:188524/140‑1 (MQ=255)                                                                           cGAACGCTGTTTTTTTGTTTTTGGTTCAATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgc                                                          <  2:273346/140‑1 (MQ=255)                                                                            gAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCa                                                              >  2:55618/1‑135 (MQ=255)                                                                            gAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCa                                                              <  1:55618/135‑1 (MQ=255)                                                                             aaCGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgcgc                                                        <  2:920069/140‑1 (MQ=255)                                                                              aCGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCa                                                       >  2:520127/1‑140 (MQ=255)                                                                                   gTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCaa                                                  <  1:924568/140‑1 (MQ=255)                                                                                         ttgtttttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGgtgt                                            <  2:635454/140‑1 (MQ=255)                                                                                                tgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCg                                     <  2:973591/140‑1 (MQ=255)                                                                                                tgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCg                                     <  1:520127/140‑1 (MQ=255)                                                                                                            actcactcTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACt                         <  1:789818/140‑1 (MQ=255)                                                                                                              tcactcTTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGa                       <  1:401652/140‑1 (MQ=255)                                                                                                                 ctctTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATgg                    >  2:802009/1‑140 (MQ=255)                                                                                                                 ctctTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATgg                    >  2:404072/1‑140 (MQ=255)                                                                                                                  tctTTTTTGAATATCAATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGt                   >  2:19101/1‑140 (MQ=255)                                                                                                                   ctTTTTTGAATATCCATCCCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGtt                  <  2:1147191/140‑1 (MQ=255)                                                                                                                                   tccCT‑GGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAg  <  1:607337/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        TCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAG  >  NC_000913/1555566‑1555836

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