Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,829,972	T→G	9.3%	H192Q (CAT→CAG)	PP_2482 →	molybdenum cofactor biosynthesis protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,829,972	0	T	G	9.3%	152.0 / 7.1	54	H192Q (CAT→CAG)	PP_2482	molybdenum cofactor biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base T (25/24);  new base G (2/3);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

GCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAGCGCGGGGCGT  >  NC_002947/2829832‑2830120                                                                                                                                             |                                                                                                                                                       gcgcgcCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGcc                                                                                                                                                  >  1:213060/1‑147 (MQ=255) gcgcgcCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGcc                                                                                                                                                  >  1:36888/1‑147 (MQ=255) gcgcgcCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCa                                                                                                                                                <  2:72054/149‑1 (MQ=255)    cgcCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCAGATGGCCCatgc                                                                                                                                             <  1:47554/149‑4 (MQ=255)       cGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATg                                                                                                                                           >  8:58415/1‑148 (MQ=255)       cGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGc                                                                                                                                          <  7:58415/149‑1 (MQ=255)              gCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCt                                                                                                                                   >  2:306634/1‑149 (MQ=255)                      aTCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑Ac                                                                                                                          >  1:72936/1‑149 (MQ=255)                        cAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑Acc                                                                                                                         >  2:112449/1‑148 (MQ=255)                          aGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGGGCTTGTTGGCGTTCTGGCTTCTGGAGGGCTAGGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAg                                                                                                                      <  4:134141/149‑1 (MQ=255)                           gATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCa                                                                                                                       >  7:166749/1‑147 (MQ=255)                           gATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCa                                                                                                                       <  8:166749/147‑1 (MQ=255)                                aCATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTtcgtc                                                                                                                <  1:306634/149‑1 (MQ=255)                                       gCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGc                                                                                                                                          <  3:205491/117‑1 (MQ=255)                                       gCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGc                                                                                                                                          >  4:205491/1‑117 (MQ=255)                                        ccGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAgcg                                                                                                        >  7:149864/1‑149 (MQ=255)                                         cGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAgcgc                                                                                                       <  2:36888/149‑1 (MQ=255)                                          gATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTAGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAgcgc                                                                                                       >  5:150610/1‑148 (MQ=255)                                               gTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAgctg                                                                                                <  8:77387/149‑1 (MQ=255)                                                  gTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGctgctg                                                                                             <  2:72936/149‑1 (MQ=255)                                                     gCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAg                                                                                          <  2:19748/149‑1 (MQ=255)                                                      cAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGc                                                                                         <  1:204440/149‑1 (MQ=255)                                                               tCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGc                                                                                 >  3:120309/1‑148 (MQ=255)                                                               tCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCAGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCg                                                                                <  1:112449/149‑1 (MQ=255)                                                                      ttGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGtt                                                                          <  2:213060/148‑1 (MQ=255)                                                                        gCCCATGTTGGCGTTCTGCCTTGAGCAGGGGTACGAGTTGCGTTTTATCGAACTTATGTGCATGGGGCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACGCAGTTCGTCGGCCAGCGCAGAGCTGCTGGAGCTTATCGGCGAGCGTTa                                                                         <  1:178832/148‑1 (MQ=255)                                                                         cccATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGc                                                                      >  8:79561/1‑149 (MQ=255)                                                                          ccATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGcc                                                                     <  3:91430/149‑1 (MQ=255)                                                                          ccATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGcc                                                                     >  4:91430/1‑149 (MQ=255)                                                                                  gCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCa                                                             <  8:149864/149‑1 (MQ=255)                                                                                    gTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAgg                                                           >  7:11048/1‑149 (MQ=255)                                                                                       cTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCt                                                         >  7:79417/1‑148 (MQ=255)                                                                                          ccTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGa                                                       <  4:120309/147‑1 (MQ=255)                                                                                          ccTTCAGCAGGCCTACGAGTTGCGTTTTATCGAACTTATGCGCATTGGCCATCTGGCCCAGCATTCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAgg                                                     <  8:132765/149‑1 (MQ=255)                                                                                           cTTCAGCAGGGCAACGAGTTGCGTTTTATGGAACTTATGCGCATTGGCCATCTGGCACAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTTCTGTAGCTTTTCGGCGAGCGTTATGCCTTCACCCAGGCTGAGgc                                                    <  7:111119/149‑1 (MQ=255)                                                                                           attCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATTCGCATGCGCAATCTGGACCAGCATGCACAGGCATTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGgc                                                    >  8:111119/2‑149 (MQ=255)                                                                                            ttCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGgcg                                                   >  8:274869/1‑149 (MQ=255)                                                                                                gcagGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCgg                                               >  2:136274/1‑149 (MQ=255)                                                                                                 cagGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGt                                              <  8:11048/149‑1 (MQ=255)                                                                                                    ggCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGa                                           <  7:79561/149‑1 (MQ=255)                                                                                                       tACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTc                                        >  1:31102/1‑149 (MQ=255)                                                                                                                         gAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTtata                       <  1:136274/148‑1 (MQ=255)                                                                                                                                 gcgcATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGc               >  2:59766/1‑148 (MQ=255)                                                                                                                                 gcgcATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGc               >  8:39930/1‑148 (MQ=255)                                                                                                                                  atgctGGGCCAGATGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCAc                                                                <  4:7020/94‑1 (MQ=255)                                                                                                                                   tgctGGGCCAGATGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAg            >  1:100606/4‑149 (MQ=255)                                                                                                                                   gcATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAg            <  1:59766/149‑1 (MQ=255)                                                                                                                                   gcATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAg            >  5:37801/1‑149 (MQ=255)                                                                                                                                     ctGGGCCAGATGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTgg                >  4:216491/2‑143 (MQ=255)                                                                                                                                     ctGGGCCAGATGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTgg                <  3:216491/142‑1 (MQ=255)                                                                                                                                     aTGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAgcg          >  1:96606/1‑149 (MQ=255)                                                                                                                                        ggCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCTGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAGCGCgg       <  7:274869/149‑1 (MQ=255)                                                                                                                                        ggCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAGCGCgg       >  5:134955/1‑149 (MQ=255)                                                                                                                                         gCCATCTGGCCCAGCATGCACAGGCTTTTCGTCACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAGCGCgg       >  1:35796/1‑149 (MQ=255)                                                                                                                                             tCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAGCGCGGGGCGt  >  4:15062/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                       GCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGT‑ACCCAGTTCGTCGGCCAGCGC‑GAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAGCGCGGGGCGT  >  NC_002947/2829832‑2830120

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