Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,829,328	C→G	7.2%	G136A (GGC→GCC)	astD ←	aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,829,328	0	C	G	7.2%	151.3 / 4.9	56	G136A (GGC→GCC)	astD	aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (29/23);  new base G (2/2);  total (31/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 = 4.52e-01

AACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTT  >  NC_000913/1829188‑1829470                                                                                                                                             |                                                                                                                                               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)                                                                                                                                            gCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCCTCCCAGCGCGGCt   >  2:777971/1‑143 (MQ=255)                                                                                                                                            gCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGGACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGCTCATCACCGTCACCTCGGTTGCCGCTTCCCAGTGCGGCt   >  2:305978/1‑143 (MQ=255)                                                                                                                                             ccGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCtt  <  1:124746/143‑1 (MQ=255)                                                                                                                                             ccGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTCCTCGCCGGCACGCCCGCGATACGCCTTCCTTGATATCGCGCTCTTCTTGACCATCGCCGCCACTTCGGCTGCCGCTTCCTCGCGCCGCtt  >  1:354959/1‑143 (MQ=255)                                                                                                                                             |                                                                                                                                               AACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTT  >  NC_000913/1829188‑1829470

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