Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,436,301	C→A	7.1%	T79K (ACG→AAG)	gpI →	tail formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,436,301	0	C	A	7.1%	179.3 / 3.7	56	T79K (ACG→AAG)	gpI	tail formation protein
Reads supporting (aligned to +/- strand):  ref base C (24/28);  new base A (2/2);  total (26/30)
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.89e-01

GTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTT  >  NC_002947/3436156‑3436447                                                                                                                                                   |                                                                                                                                                   gTGGCCGACTCGTCACGCTGCCCGCCGAACTACC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGt                                                                                                                                                 >  7:318710/1‑149 (MQ=255)     ccGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgc                                                                                                                                             <  7:252962/149‑1 (MQ=255)         cTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGATGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                                                                                                                                         <  3:115716/149‑1 (MQ=255)                 gCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGTCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                                                                                                                                 <  1:77313/149‑1 (MQ=255)                 gCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                                                                                                                                 <  3:147945/149‑1 (MQ=255)                 gCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                                                                                                                                 >  2:276336/1‑149 (MQ=255)                      ccgccgAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGCGCACGGTCGGCGCGATTCGGCGGGTGc                                                                                                                             >  2:332337/1‑148 (MQ=255)                       cgccgAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTc                                                                                                                           <  1:319493/149‑1 (MQ=255)                         ccgAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTc                                                                                                                           >  4:278080/1‑147 (MQ=255)                          cgACCTTCCATGCCCAGGCCAGCCAGGGCA‑GGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                                                                                                                                         <  8:118937/132‑1 (MQ=37)                          cgACCTTCCATGCCCAGGCCAGCCAGGGCA‑GGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                                                                                                                                         >  7:118937/1‑132 (MQ=37)                          cgAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAg                                                                                                                        <  1:276336/149‑1 (MQ=255)                           gAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAgg                                                                                                                       <  8:318710/149‑1 (MQ=255)                             aCTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCg                                                                                                                     <  5:33727/149‑1 (MQ=255)                              cTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCCGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCAAGGTCGGCGCGATTCGGCGGGTGCGCAcggcgg                                                                                                                    >  6:26275/1‑149 (MQ=255)                              cTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCAAGGTCGGCGCGATTCGGCGGGTGCTCAAGGCgg                                                                                                                    >  6:21046/1‑149 (MQ=255)                                  c‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCg                                                                                                                >  7:179699/1‑149 (MQ=255)                                    tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg                                                                                                               >  4:177577/1‑149 (MQ=255)                                      cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCg                                                                                                                >  3:212791/1‑146 (MQ=255)                                      cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCg                                                                                                                <  4:212791/146‑1 (MQ=255)                                                 cTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGAtt                                                                                                  <  2:313209/149‑1 (MQ=255)                                                 cTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGAt                                                                                                   >  2:29842/1‑148 (MQ=255)                                                     gCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACaa                                                                                              <  3:177577/149‑1 (MQ=255)                                                      cATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGt                                                                                                             <  6:354393/133‑1 (MQ=255)                                                      cATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAg                                                                                             <  2:233464/149‑1 (MQ=255)                                                      cATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCAGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGt                                                                                                             >  5:354393/1‑133 (MQ=255)                                                       aTGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAg                                                                                             >  8:273647/1‑148 (MQ=255)                                                         ggAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCGCAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGaa                                                                                          >  3:329045/1‑149 (MQ=255)                                                         ggAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGaa                                                                                          >  6:322974/1‑149 (MQ=255)                                                          gAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGGTTACAAGGAAt                                                                                         >  5:156950/1‑149 (MQ=255)                                                           aaGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATg                                                                                        <  7:9681/149‑1 (MQ=255)                                                               tCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                                                                                                                                         >  5:103243/1‑96 (MQ=255)                                                               tCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                                                                                                                                         <  6:103243/96‑1 (MQ=255)                                                               tCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGgcgc                                                                                    <  8:179699/149‑1 (MQ=255)                                                                                  ttACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGgcgc                                                                 >  1:290287/1‑149 (MQ=255)                                                                                       ccGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGATTGGCGCGAGATACCCAACGCGGCGCCGTac                                                            <  1:29842/149‑1 (MQ=255)                                                                                         gacgacCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTacac                                                          <  3:278080/149‑1 (MQ=255)                                                                                               cAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTc                                                      >  5:215852/1‑147 (MQ=255)                                                                                                aGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACACGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAg                                                   <  5:112842/149‑1 (MQ=255)                                                                                                aGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGTCGTGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAg                                                   <  2:116248/149‑1 (MQ=255)                                                                                                     cgcgcgCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCAc                                              >  7:193036/1‑149 (MQ=255)                                                                                                     cgcgcgCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCAc                                              <  2:290287/149‑1 (MQ=255)                                                                                                        gcgcTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGc                                           >  1:8726/1‑149 (MQ=255)                                                                                                            tGATCCGCGATGACATTCCGGTTCACAAGACCAACGGCACGGTCGGCGCGATTCGGGGGGTGCGTAAGGCGGTTCGGGTCAACTCGGATTACAAGGTATGTCGCGAGATACACAACTCGGCGCCGTACACGTTCCAGGTCAc                                              <  2:260296/142‑1 (MQ=255)                                                                                                            tGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGTCGCGAGATAACCAACGCGGCGCAGTACACTTTCCAGGTCACgagatag                                       >  1:260296/1‑143 (MQ=255)                                                                                                              aTCCGCGAGGCCATTCCGGTTCCCCAGACCAAGGGCAAGGGCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGGCACGGCTTggg                                      <  5:26275/148‑1 (MQ=255)                                                                                                              aTCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGGTTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTggg                                      <  6:156950/148‑1 (MQ=255)                                                                                                              aTCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTgg                                       <  8:193036/147‑1 (MQ=255)                                                                                                              aTCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCAAGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTggg                                      <  5:21046/148‑1 (MQ=255)                                                                                                                   cgAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACg                                <  2:8726/149‑1 (MQ=255)                                                                                                                     aGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACgag                              <  7:273647/149‑1 (MQ=255)                                                                                                                                      aaGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAgg             >  5:197584/1‑149 (MQ=255)                                                                                                                                           cAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGa        >  5:73687/1‑149 (MQ=255)                                                                                                                                                gCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAAccggccgg                    <  4:48140/132‑1 (MQ=255)                                                                                                                                                gCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAAccggccgg                    >  3:48140/1‑132 (MQ=255)                                                                                                                                                gCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCAtt   >  5:314522/1‑149 (MQ=255)                                                                                                                                                 ctcGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGGTCGGGTCAACGCGGATTACAAGGATTGGCGCTTGATACCCAACGCGGCGCCGTACACGTTCCAGGTAACGTCTTGGGCCTACGAGAACCGGCCGGGTGAAGGCTCGATCAttt  <  6:207316/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTT  >  NC_002947/3436156‑3436447

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: ‑