Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,551,799	C→T	I154I (ATC→ATT)	fdnI →	formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,551,799	0	C	T	100.0%	241.5 / NA	74	I154I (ATC→ATT)	fdnI	formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/40);  total (34/40)

CAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCAT  >  NC_000913/1551664‑1551928                                                                                                                                        |                                                                                                                                  cAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCAc                                                                                                                                <  1:716515/139‑1 (MQ=255)    cGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCAcgcg                                                                                                                             >  2:106143/1‑139 (MQ=255)         ggCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGc                                                                                                                        >  1:1225906/1‑139 (MQ=255)         ggCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGc                                                                                                                        >  1:1227972/1‑139 (MQ=255)         ggCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGc                                                                                                                        <  2:1227972/139‑1 (MQ=255)              aaGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTa                                                                                                                   <  2:1225906/139‑1 (MQ=255)                   gatgTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTatcatc                                                                                                              <  1:1442078/139‑1 (MQ=255)                    atgTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCgg                                                                                                                            >  1:1607632/1‑124 (MQ=255)                    atgTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCgg                                                                                                                            <  2:1369753/124‑1 (MQ=255)                    atgTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCgg                                                                                                                            >  1:1369753/1‑124 (MQ=255)                    atgTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCgg                                                                                                                            <  2:1607632/124‑1 (MQ=255)                        tCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGAt                                                                                                         >  1:1019230/1‑139 (MQ=255)                          tGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATcc                                                                                                       >  2:522160/1‑139 (MQ=255)                             tCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACg                                                                                                    <  2:57743/139‑1 (MQ=255)                              cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGc                                                                                                   <  2:1029458/139‑1 (MQ=255)                                 tCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCAt                                                                                                >  1:1055952/1‑139 (MQ=255)                                  cATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATc                                                                                               <  1:106143/139‑1 (MQ=255)                                      gatacagttttCGTGCTGCTTGTGACCGGTGTTATTATCTGGCTTCCGTACTTTGCGCAGTACTTCCCGATGCATGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCCGGTATCATCCTGATCCACGCCATCCTGa                                                                                           <  2:1050382/132‑1 (MQ=255)                                      aGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGa                                                                                           <  2:1019230/139‑1 (MQ=255)                                           gATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCat                                                                                      >  1:636047/1‑139 (MQ=255)                                            aTTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGa                                                                                           <  2:176249/133‑1 (MQ=255)                                            aTTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGa                                                                                           >  1:176249/1‑133 (MQ=255)                                                   tgctgcTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTatat                                                                              <  1:995568/139‑1 (MQ=255)                                                        ctgGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCAt                                                                         >  2:35019/1‑139 (MQ=255)                                                        ctgGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGAAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATAATCCTGATCCACGCCATCCTGATCCATATGTATCTgtctc                                                                         >  1:1050382/1‑135 (MQ=255)                                                                 ggggTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGa                                                                <  1:522160/139‑1 (MQ=255)                                                                  gggTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTa                                                                                                                   <  1:557511/87‑1 (MQ=255)                                                                  gggTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTa                                                                                                                   >  2:557511/1‑87 (MQ=255)                                                                  gggTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATc                                                                                         >  2:843649/1‑113 (MQ=255)                                                                  gggTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATc                                                                                         >  2:274483/1‑113 (MQ=255)                                                                  gggTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATc                                                                                         <  1:843649/113‑1 (MQ=255)                                                                  gggTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATc                                                                                         <  1:274483/113‑1 (MQ=255)                                                                  gggTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGaa                                                               >  2:726141/1‑139 (MQ=255)                                                                       aTTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCgg                                                                                                                      >  1:1321031/1‑79 (MQ=255)                                                                       aTTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCgg                                                                                                                      <  2:1321031/79‑1 (MQ=255)                                                                        ttATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATc                                                         <  2:1055952/139‑1 (MQ=255)                                                                        ttATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATc                                                         <  1:35019/139‑1 (MQ=255)                                                                                  tCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGa                                                                <  2:1660722/122‑1 (MQ=255)                                                                                  tCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGa                                                                >  1:1660722/1‑122 (MQ=255)                                                                                   ccGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGga                                              <  2:636047/139‑1 (MQ=255)                                                                                    cGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTat                                                                                                                  <  1:854135/70‑1 (MQ=255)                                                                                    cGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTat                                                                                                                  >  2:854135/1‑70 (MQ=255)                                                                                    cGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTa                                                                                                                   >  1:500769/1‑69 (MQ=255)                                                                                    cGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTa                                                                                                                   <  2:500769/69‑1 (MQ=255)                                                                                    cGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATcc                                                                                              >  1:326361/1‑90 (MQ=255)                                                                                    cGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATcc                                                                                              <  2:326361/90‑1 (MQ=255)                                                                                          ttGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCt                                                                                             >  1:1481179/1‑85 (MQ=255)                                                                                          ttGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCt                                                                                             <  2:1481179/85‑1 (MQ=255)                                                                                           tGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATc                                                                                         >  2:230831/1‑88 (MQ=255)                                                                                           tGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATc                                                                                         <  1:230831/88‑1 (MQ=255)                                                                                           tGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCgaa                                      <  1:726141/139‑1 (MQ=255)                                                                                             cgcAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGAt                                                                                                         >  1:1388001/1‑70 (MQ=255)                                                                                             cgcAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGAt                                                                                                         <  2:1388001/70‑1 (MQ=255)                                                                                             cgcAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATc                                                                                         >  2:943328/1‑86 (MQ=255)                                                                                             cgcAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATc                                                                                         <  1:943328/86‑1 (MQ=255)                                                                                                  tACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAggg                                               >  1:438976/1‑123 (MQ=255)                                                                                                  tACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAggg                                               <  2:438976/123‑1 (MQ=255)                                                                                                       cccGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTat                                                                                                                  <  2:146103/51‑1 (MQ=255)                                                                                                       cccGATGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTat                                                                                                                  >  1:146103/1‑51 (MQ=255)                                                                                                            tGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACg                                                                                                    <  1:1135752/60‑1 (MQ=255)                                                                                                            tGCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACg                                                                                                    >  2:1135752/1‑60 (MQ=255)                                                                                                             gCAGGTTGTTCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGc                    <  2:876144/139‑1 (MQ=255)                                                                                                                 gttgttCGCTACAGCCTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTggg                <  2:1692365/139‑1 (MQ=255)                                                                                                                               ccTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGa                                                                                           >  2:1374069/1‑50 (MQ=255)                                                                                                                               ccTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGa                                                                                           <  1:1374069/50‑1 (MQ=255)                                                                                                                               ccTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGga                                              >  2:386395/1‑95 (MQ=255)                                                                                                                               ccTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGga                                              <  1:386395/95‑1 (MQ=255)                                                                                                                               ccTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCAt  >  1:502376/1‑139 (MQ=255)                                                                                                                               ccTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGACGTCGCTGGGCgaa            >  1:1337177/1‑129 (MQ=255)                                                                                                                               ccTGCTGATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGACGTCGCTGGGCgaa            <  2:1337177/129‑1 (MQ=255)                                                                                                                                  gctgATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATcc                                                                                              >  2:835417/1‑44 (MQ=255)                                                                                                                                  gctgATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATcc                                                                                              <  1:835417/44‑1 (MQ=255)                                                                                                                                   ctgATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGa                                                       >  1:1244414/1‑82 (MQ=255)                                                                                                                                   ctgATTCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGa                                                       <  2:1244414/82‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  CAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCAT  >  NC_000913/1551664‑1551928

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