Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,904,248:1	+A	7.7%	coding (495/1593 nt)	ilvA‑I ←	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,904,248	1	.	A	7.7%	198.3 / 11.2	52	coding (495/1593 nt)	ilvA‑I	threonine deaminase
Reads supporting (aligned to +/- strand):  ref base . (25/23);  new base A (2/2);  total (27/25)
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.72e-01

GATCCCGGCGATCAAGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGCCTTCATACCC  >  NC_002947/3904107‑3904395                                                                                                                                                |                                                                                                                                                    gATCCCGGCGATCAAGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGc                                                                                                                                              >  5:101530/1‑149 (MQ=255)  aTCCCGGCGATCAAGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCt                                                                                                                                             >  4:292801/1‑149 (MQ=255)            tCAAGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGc                                                                                                                                   <  2:136142/149‑1 (MQ=255)              aaGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGG‑CGCGACCTTCgtg                                                                                                                                             >  1:100609/1‑135 (MQ=255)              aaGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGG‑CGCGACCTTCgtg                                                                                                                                             <  2:100609/137‑3 (MQ=255)              aaGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGG‑CGCGACCTTCgt                                                                                                                                              <  5:110806/136‑2 (MQ=255)              aaGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGG‑CGCGACCTTCgt                                                                                                                                              >  6:110806/1‑135 (MQ=255)                  cgccACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGt                                                                                                                             >  1:276412/1‑149 (MQ=255)                      aCCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGt                                                                                                                         >  6:256380/1‑149 (MQ=255)                      aCCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGt                                                                                                                         <  6:140836/149‑1 (MQ=255)                       ccGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTg                                                                                                                        <  3:79395/149‑1 (MQ=255)                        cGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGTGCGCTGACGGAGGATTTCCATGGCGACGGTTCCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CGTTCGCTGTCGGCCAGCTTCAGTGCGTgg                                                                                                                       <  7:7458/149‑1 (MQ=255)                        cGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTgg                                                                                                                       <  6:101530/149‑1 (MQ=255)                              cGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAgc                                                                                                                 <  3:292801/149‑1 (MQ=255)                               gggCACAAAGCTCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAgcg                                                                                                                <  8:164713/149‑1 (MQ=255)                                 gCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGc                                                                                                                                              >  6:102596/1‑117 (MQ=255)                                 gCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGc                                                                                                                                              <  5:102596/117‑1 (MQ=255)                                   acaAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATg                                                                                                            >  3:118232/1‑149 (MQ=255)                                    caAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGc                                                                                                           <  8:41447/149‑1 (MQ=255)                                     aaaGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCg                                                                                                          >  6:100230/1‑149 (MQ=255)                                          tCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAAGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAg                                                                                                     <  4:6165/149‑1 (MQ=255)                                              gTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGctct                                                                                                 >  1:57730/1‑149 (MQ=255)                                                  aaCGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGACGCGCTCggc                                                                                             >  6:251140/1‑147 (MQ=255)                                                  aaCGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCtctc                                                                                                >  3:169682/1‑146 (MQ=255)                                                  aaCGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCtctc                                                                                                <  4:169682/146‑1 (MQ=255)                                                   aCGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCa                                                                                            >  1:89888/1‑149 (MQ=255)                                                         cAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATgcagc                                                                                      <  2:276412/149‑1 (MQ=255)                                                                     ggaggaTTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCccg                                                                          >  2:253035/1‑149 (MQ=255)                                                                      gaggaTTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCCGGGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTAGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCccg                                                                          >  1:255978/1‑149 (MQ=255)                                                                       aggaTTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCccgcc                                                                        >  7:95894/1‑149 (MQ=255)                                                                                aTGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAacg                                                               >  2:57488/1‑149 (MQ=255)                                                                                 tGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAacg                                                               >  7:302885/1‑148 (MQ=255)                                                                                    ccACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTgg                                                                                                                       >  8:136724/1‑89 (MQ=255)                                                                                    ccACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTgg                                                                                                                       <  7:136724/89‑1 (MQ=255)                                                                                     cACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAacgcacgc                                                          >  1:60268/1‑149 (MQ=255)                                                                                        ggTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAacgcac                                                            >  2:24827/1‑144 (MQ=255)                                                                                        ggTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAacgcac                                                            <  1:24827/144‑1 (MQ=255)                                                                                             ccTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACt                                                  <  2:57730/149‑1 (MQ=255)                                                                                              cTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACtt                                                 >  4:98111/1‑149 (MQ=255)                                                                                                 cccGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTc                                               <  1:57488/148‑1 (MQ=255)                                                                                                     gCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGa                                          >  6:345606/1‑149 (MQ=255)                                                                                                       gATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATg                                        <  2:89888/149‑1 (MQ=255)                                                                                                       gATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATg                                        <  2:60268/149‑1 (MQ=255)                                                                                                            caTCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCgtg                                   <  8:302885/149‑1 (MQ=255)                                                                                                                    gTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGtggtggtgg                            <  4:77467/149‑1 (MQ=255)                                                                                                                    gTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGtggtggtgg                            <  1:253035/149‑1 (MQ=255)                                                                                                                          gAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGa                      >  4:264904/1‑149 (MQ=255)                                                                                                                               gTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATg                 >  2:130102/1‑149 (MQ=255)                                                                                                                                    aCGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGc            >  6:271333/1‑149 (MQ=255)                                                                                                                                      gAAGGGCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGCCt          <  8:126594/149‑1 (MQ=255)                                                                                                                                        aGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGCCtt         >  4:180635/1‑148 (MQ=255)                                                                                                                                         ggTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGCCTTCa       <  3:98111/149‑1 (MQ=255)                                                                                                                                             gcg‑cCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGtggtggtgg                            <  4:136170/124‑1 (MQ=255)                                                                                                                                             gcg‑cCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGtggtggtgg                            >  3:136170/1‑124 (MQ=255)                                                                                                                                              cg‑cCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGCCTTCATAccc  <  8:95894/149‑1 (MQ=255)                                                                                                                                               g‑cCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGcc           >  7:206888/1‑139 (MQ=255)                                                                                                                                               g‑cCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGcc           <  8:206888/139‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    GATCCCGGCGATCAAGCCGCCACCGCCCACGGGCACAAAGATCGCGTCCAACGCCCCAGGGCGCTGACGGAGGATTTCCATGGCCACGGTACCCTGCCCGGCGATCACATCC‑GGGTCATCGAACGGTGGCACGAAGGTCGCG‑CCTTCGCTGTCGGCCAGCTTCAGTGCGTGGGTCAGCGCATGCGGGAAGCTCTCGCCATGCAGCACCACATGCCCGCCGCGCGAACGCACGCCTTCCACTTTCAACGATGGCGTGGTGGTGGGCATGACGATGGTTGCCTTCATACCC  >  NC_002947/3904107‑3904395

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