Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,829,322	A→C	6.8%	L138R (CTG→CGG)	astD ←	aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,829,322	0	A	C	6.8%	131.0 / 2.3	59	L138R (CTG→CGG)	astD	aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (32/23);  new base C (2/2);  total (34/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 = 3.03e-01

CTGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCC  >  NC_000913/1829180‑1829459                                                                                                                                               |                                                                                                                                          ctgcCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAACTATACGGCCCAACCACCGCCa                                                                                                                                           >  2:772434/1‑143 (MQ=255)    ccATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGATGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGACATTATACGGCCCCAACACCGCCAGCa                                                                                                                                        >  1:58384/1‑143 (MQ=255)      aTAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGcgcg                                                                                                                                      >  2:812716/1‑143 (MQ=255)       tAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATCCGGCCCAAACACCGCCAGCActc                                                                                                                                     >  1:278513/1‑141 (MQ=255)         aaCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTCTCAAGATAATGGCGTTAAAAGAATGCGATGGCGGAACGATATGTCCGTTCGGCAAATGACCAGGGCAATTATACGGCCCAAACACCGCCAGCACGcgg                                                                                                                                   >  2:844186/1‑141 (MQ=255)          aCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAAGTATACGGCCCAAGCACCGCCAGCACGCCGt                                                                                                                                  >  2:916010/1‑143 (MQ=255)          aCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGGTTAAGGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGCCAAATGACCAGGGAAATTCTACTGCCAAAACACCGCGAGCccgccgc                                                                                                                                  >  2:484840/1‑142 (MQ=255)               tgaCCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGGACGACATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCggg                                                                                                                             <  1:127277/141‑1 (MQ=255)                  ccGCTCCGCCACTCCACGCTGTCAGTTCGCTGGGTTTAAAGATATTGGCGTCACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAACGACCAGGGAAATTATACGGCCCAAACACCGCCCGCACGGCGTGCTGGCgg                                                                                                                          <  1:894137/143‑2 (MQ=255)                   cGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGac                                                                                                                         >  2:677329/1‑142 (MQ=255)                     cTTCGCCACTCGCCGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAACACCCGATAGTACGCCGTGCGGGCGATGt                                                                                                                       <  2:111286/143‑1 (MQ=255)                     cTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGt                                                                                                                       <  1:50162/143‑1 (MQ=255)                        cGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGGACGCCGTGCGGGCGATGTCGc                                                                                                                    >  1:163202/1‑143 (MQ=255)                        cGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGCCGTTCCCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACATCGCCCGCACGGCGTGCTGGCGgtgtttg                                                                                                                    <  1:205019/143‑4 (MQ=255)                          cctcTCCAGGGTGTCAGTGTGCTGGTGTTAAAGATAATGGTGTTGCCTGCCAGCAATGCCGGAACGATATGTCCGTCCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAg                                                                                                                  <  2:859193/140‑1 (MQ=255)                                    gtgtCAGTTCGCTTGGTTTCTATATATTGGTGTCACCTGCCAGCAATGCCGGAACGATATGTCCGGGCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccgcg                                                                                                        <  1:6555/143‑1 (MQ=255)                                       tCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGCCGATGTCCCAGGCTcgccccaccc                                                                                                     >  2:774905/1‑137 (MQ=255)                                           ttCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCCGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGCCCg                                                                                                 >  2:790173/1‑143 (MQ=255)                                           ttCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCg                                                                                                 >  2:124746/1‑143 (MQ=255)                                           ttCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCg                                                                                                 >  2:379518/1‑143 (MQ=255)                                              gCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAACTGACCCGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGGAGGCTcgccgcggcgccaccca                                                                                              >  2:628764/1‑136 (MQ=255)                                                   gTTGAAAGAGAATGGTGTTACCGGCCAGCAATGCCGGAACGATATGTCCGGTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCa                                                                                         <  1:893407/143‑1 (MQ=255)                                                        aaGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATAGGTCCGTTCGGCAAATGACCAGGGAAATTATACGGACGAAACACCGCCAGCACGCCGTGCGGGCGAGGTCGCAGGCTCGCCGCCCCGTCCGTCATTTCACTACg                                                                                    >  1:430802/1‑143 (MQ=255)                                                         aGACAATGGTGTGATTTGCCAGCAATGCCGTAACTATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgc                                                                                   <  2:40560/143‑1 (MQ=255)                                                          gATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCCCAGGCTCGCCCCGCCATCCGCCCTTCCACTACgcc                                                                                  >  2:421149/1‑142 (MQ=255)                                                           aTAATGGTGTTACATGACAGCAATGCCGGAAAAATAGGTCCGTTCGGCAAATGACCAGGGCAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGACGCAGGCCCGCCGCGCCCGCCGGCACTTCAcgaggcgc                                                                                 >  2:109518/1‑136 (MQ=255)                                                               tGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGCCCCGGGATATTATACGTCCCAAACATCGCCAGAACGCAGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCg                                                                             >  1:131159/1‑143 (MQ=255)                                                                   gtgaCCGGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCgg                                                                         <  1:685885/140‑1 (MQ=255)                                                                   gtTCCCTGCCAGCCATGCCGGAACGATCTGGCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCgg                                                                         <  1:40668/143‑1 (MQ=255)                                                                    ttACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCTAAAAACCACCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCCCTACGCTGCTCGcccgc                                                                        >  1:126947/1‑140 (MQ=255)                                                                    ttACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCACACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGt                                                                        <  1:866662/143‑1 (MQ=255)                                                                      aCCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCCTTTCCCTACCCCGCTCGCCGGTAc                                                                      >  1:228728/1‑143 (MQ=255)                                                                         tGCCAGCAATGCCCGACCGACATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGaa                                                                   <  2:162167/143‑1 (MQ=255)                                                                          gCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGACCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTGACCACACTACTCGCCGGTACGacc                                                                  >  2:806088/1‑141 (MQ=255)                                                                                    ccGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATCCGcc                                                        >  1:794458/1‑143 (MQ=255)                                                                                     cGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCAGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTGCGAACGGGATAGGCCt                                                       >  2:129387/1‑143 (MQ=255)                                                                                      ggAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTGCGAACGTGGTACGCCtt                                                      >  1:425043/1‑143 (MQ=255)                                                                                       gAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGAGGCGCGCAGGTACGAACGTGGTACGCCtga                                                     >  1:632225/1‑141 (MQ=255)                                                                                        aaCGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTaa                                                    <  1:80794/143‑1 (MQ=255)                                                                                          cGATATGTCCGCTCGGCAAATGACCAAAAGCAGACGTCTGCCCAAGCACCGCCAGCACGCCGTGCGGGCGATGCCGCAGTCTCGCCGCGCCGTCCGGCATTTCACTACGCCGCTCGTCGGTACGAACGTGATACGCCTTAAtt                                                  <  1:164880/143‑1 (MQ=255)                                                                                          cGATATGCTCGTGCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCGGGGCGGGCGAGGTCGCAGGCGCGCCGCGCCGGTCGGCATTTCACTACGCTGCTCGTCGGTACGAACGTGATACGCCTTAAtt                                                  <  1:843132/143‑1 (MQ=255)                                                                                             tatGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGat                                               >  2:856934/1‑143 (MQ=255)                                                                                                 tCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATcg                                           <  1:629239/143‑1 (MQ=255)                                                                                                   cGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTCCTTGCTATcgcg                                         >  1:289976/1‑143 (MQ=255)                                                                                                    gTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGTTGTCGCAGGCTCGCCGCGCCGTTCAGCATTTCATTACTCTGCTCGCCGGTACGAACGTGATACGCCTTCATTGATATAGCGa                                        >  1:146394/1‑143 (MQ=255)                                                                                                     ttCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAt                                       >  1:731778/1‑143 (MQ=255)                                                                                                        ggCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAtttt                                    >  1:634387/1‑143 (MQ=255)                                                                                                           aaaTGACCGGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCGCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATGTTAtt                                 <  2:48510/143‑1 (MQ=255)                                                                                                           aaaTGACCGGGAACATGACACGGCCCACACACCGCCAGCACGCCGTGCGGGCGATGCCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCTGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAtt                                 <  2:348740/143‑1 (MQ=255)                                                                                                           aaaTGACCAGGTAATTTATCCGGCCCAATCACCGCCAGCACGCCGTGCGGGCCAGTTCGCAGGCTCGCCGCGCCGTCCGTCATTTCTCTACGCTGCTCGCCGGTACGAACTTGATACGCCTTAATTGATATCGCGATTTTAtt                                 <  2:751411/143‑1 (MQ=255)                                                                                                                   aGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTCTTGATCATCg                         >  2:558685/1‑143 (MQ=255)                                                                                                                              aCGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAGTTTATTGATCATCGCCGTCACGTCg              >  1:496722/1‑143 (MQ=255)                                                                                                                               cGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCgg             <  1:718497/143‑1 (MQ=255)                                                                                                                                  cccAGACCCCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTg          <  1:687967/143‑1 (MQ=255)                                                                                                                                   tcatacacCGCCAGCACGCTGTTGGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTATGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGc         <  2:563501/139‑1 (MQ=255)                                                                                                                                      agcacCGCCAGCACGCCGGGCGGGCGAGGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGc      <  1:348956/141‑1 (MQ=255)                                                                                                                                      aaCACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACTGGATACGCCTTACTTGATATCGCGATTCTCTTGATCATCGCCGTCACTACGGTAGccac      >  2:448082/1‑141 (MQ=255)                                                                                                                                          tcGCCCGCACGGCGTGCTGGCGGTGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACGTCGGTTGCCGCTTcc  >  2:671138/2‑143 (MQ=255)                                                                                                                                          tcGCCCGCACGGCGTGCTGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCCGCTCGCCGGTACGAACGCGATACGCCTTAATTGATATCGCGACTTTATTGCTCATCGCCGTCACCTCGGTTGCCGCTTcc  >  2:345286/2‑143 (MQ=255)                                                                                                                                          ccGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTcc  <  2:699971/143‑1 (MQ=255)                                                                                                                                               |                                                                                                                                          CTGCCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCC  >  NC_000913/1829180‑1829459

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