Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,676,674	A→G	7.0%	Y634H (TAC→CAC)	ldcC ←	lysine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,676,674	0	A	G	7.0%	165.1 / 2.1	57	Y634H (TAC→CAC)	ldcC	lysine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (26/27);  new base G (2/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACACGCCCCAGCAATGCCTCGATCGGCACCGCCTCCACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTCATAGTGGCGCTTGAACTCCAGCAGCTCGGTGAGCAGG  >  NC_002947/4676532‑4676819                                                                                                                                               |                                                                                                                                                  acacGCCCCAGCAATGCCTCGATCGGCACCGCCTCCACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGt                                                                                                                                              >  5:10144/1‑148 (MQ=255)          aGCAATGCCTCGATCGGCACCGCCTCCACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTg                                                                                                                                    <  4:221206/149‑1 (MQ=255)           gCAATGCCTCGATCGGCACCGCCTCCACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTgg                                                                                                                                   >  1:66/1‑149 (MQ=255)            cAATGCCTCGATCGGCACCGCCTCCACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTg                                                                                                                                    >  4:50756/1‑147 (MQ=255)            cAATGCCTCGATCGGCACCGCCTCCACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTg                                                                                                                                    <  3:50756/147‑1 (MQ=255)                  cTCGATCGGCACCGCCTCCACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGcac                                                                                                                             <  8:118142/148‑1 (MQ=255)                  cTCGATCGGCACCGCCTCCACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGcac                                                                                                                             >  7:118142/1‑148 (MQ=255)                                ccTCCACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGcc                                                                                                               >  3:134929/1‑148 (MQ=255)                                 cTCCACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCa                                                                                                             <  6:10144/149‑1 (MQ=255)                                   ccACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATg                                                                                                           >  1:51785/1‑149 (MQ=255)                                    cACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGc                                                                                                                <  1:261741/143‑1 (MQ=255)                                    cACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGc                                                                                                                >  2:261741/1‑143 (MQ=255)                                     aCTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGcc                                                                                                         >  2:36674/1‑149 (MQ=255)                                     aCTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGcc                                                                                                         <  4:134929/149‑1 (MQ=255)                                      cTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGccc                                                                                                        <  4:63955/149‑1 (MQ=255)                                         cGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCt                                                                                                       <  3:18465/147‑1 (MQ=255)                                             aCGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAgc                                                                                                 <  3:88548/149‑1 (MQ=255)                                                 aCCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAgcgc                                                                                               >  2:233384/1‑147 (MQ=255)                                                  ccATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAg                                                                                            <  5:125368/149‑1 (MQ=255)                                                         gTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCgg                                                                                     >  3:143753/1‑149 (MQ=255)                                                          tCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGc                                                                                    >  6:3603/1‑149 (MQ=255)                                                            aTAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAg                                                                                                                                               <  4:234865/88‑1 (MQ=255)                                                            aTAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAg                                                                                  >  3:109543/1‑149 (MQ=255)                                                             tAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGc                                                                                 >  8:278144/1‑149 (MQ=255)                                                                gCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGcccc                                                                              >  7:241131/1‑149 (MQ=255)                                                                gCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGcccc                                                                              <  2:66/149‑1 (MQ=255)                                                                  cTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCa                                                                                                                   >  6:41461/1‑110 (MQ=255)                                                                  cTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCa                                                                                                                   <  5:41461/110‑1 (MQ=255)                                                                       ccGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTg                                                                       <  7:278144/149‑1 (MQ=255)                                                                        cGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTgg                                                                      <  8:241131/149‑1 (MQ=255)                                                                          ggCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGggc                                                                    <  2:51785/149‑1 (MQ=255)                                                                             tGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGggcagg                                                                 >  7:128067/1‑149 (MQ=255)                                                                                   ccACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTgg                                                                                                     <  5:277597/107‑1 (MQ=255)                                                                                   ccACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTgg                                                                                                     >  6:277597/1‑107 (MQ=255)                                                                                       cTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAgctgct                                                         >  4:59678/1‑147 (MQ=255)                                                                                          cGGCAAACGGGTGATCAGCCGCTTCAGCTGCTTTGCCGTGGCATTTGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGGCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTg                                                    <  4:85459/149‑1 (MQ=255)                                                                                          cGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTg                                                    >  6:45538/1‑149 (MQ=255)                                                                                          cGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCCCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTg                                                    <  3:268763/149‑1 (MQ=255)                                                                                            gCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTggg                                                  <  5:3603/149‑1 (MQ=255)                                                                                                   ggTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGc                                           >  6:242591/1‑149 (MQ=255)                                                                                                    gTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGcc                                          <  1:233384/149‑1 (MQ=255)                                                                                                          aGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTCa                                      >  7:177383/1‑147 (MQ=255)                                                                                                                ttCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTCATAGTGgcg                              >  5:165107/1‑149 (MQ=255)                                                                                                                  cAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTCATAGTGgcgc                             <  6:165107/148‑1 (MQ=255)                                                                                                                   aGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTCATAGTGGCGCtt                           <  8:128067/149‑1 (MQ=255)                                                                                                                   aGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTCATAGTGGCGCt                            <  8:177383/148‑1 (MQ=255)                                                                                                                    gctgctTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTCATAGTGGCGCtt                           >  5:190501/1‑148 (MQ=255)                                                                                                                            gCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGGGGGGTGTTGCCGTCATAGTGGCGCTTGAACTCcag                  >  2:219758/1‑149 (MQ=255)                                                                                                                              cGTGGCATTGGCGCGGTGGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTTCTGAGTGGGGTGTTGCCGTc                                       >  5:281378/1‑126 (MQ=255)                                                                                                                              cGTGGCATTGGCGCGGTGGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTc                                       >  5:285869/1‑126 (MQ=255)                                                                                                                              cGTGGCATTGGCGCGGTGGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTc                                       <  6:285869/126‑1 (MQ=255)                                                                                                                              cGTGGCATTGGCGCGGTGGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTACTGAGTGGGGTGTTGCCGTc                                       <  6:281378/126‑1 (MQ=255)                                                                                                                                   cATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTCATAGTGGCGCTTGAACTCCAGCAGCTCg           >  5:223037/1‑149 (MQ=255)                                                                                                                                     ttGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTCATAGTGGCGCTTGAACTCCAGCAGCTCGGt         <  8:253080/149‑1 (MQ=255)                                                                                                                                     ttGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTCATAGTGGCGCTTGAACTCCAGCAGCTCGGt         <  6:223037/149‑1 (MQ=255)                                                                                                                                          gcgGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTCATAGTGGCGCTTGAACTCCAGCAGCTCGGTGAg      >  4:100006/1‑147 (MQ=255)                                                                                                                                            ggTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTCATAGTGGCGCTTGAACTCCAGCAGCTCGGTGAGCAgg  <  3:59678/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                  ACACGCCCCAGCAATGCCTCGATCGGCACCGCCTCCACTTCGCCACGCACCATCTTGTCATAGGCCTGGGCCGGGCTGACGGCCACCTCCGGCAAACGGGTGAACAGCCGCTTCAGCTGCTTGGCCGTGGCATTGGCGCGGTAGCAGTCATGTAACTGGTCGCACAGGTCGCGCAGGCCCATGCCCTGGTAGCGCAAGGCGTCGGCAGCCCCCACACTGGGCAGGCAGCTGCTGAGTGGGGTGTTGCCGTCATAGTGGCGCTTGAACTCCAGCAGCTCGGTGAGCAGG  >  NC_002947/4676532‑4676819

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