Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	835,968	G→C	6.8%	L37L (CTC→CTG)	PP_0718 ←	sulfate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	835,968	0	G	C	6.8%	186.4 / 3.2	59	L37L (CTC→CTG)	PP_0718	sulfate transporter
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base C (2/2);  total (29/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.15e-01

TCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGCGAGAAGGGCGCTGAACCTC  >  NC_002947/835828‑836113                                                                                                                                              |                                                                                                                                                  tCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAg                                                                                                                                           >  1:143376/1‑149 (MQ=255) tCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAg                                                                                                                                           <  3:350994/149‑1 (MQ=255)   cGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGc                                                                                                                                         >  5:231164/1‑149 (MQ=255)        ccAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCcggcg                                                                                                                                    >  7:47688/1‑149 (MQ=255)           gCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGgcggcg                                                                                                                                 <  1:180724/149‑1 (MQ=255)             aCCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGAt                                                                                                                               <  2:143376/149‑1 (MQ=255)             aCCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGAt                                                                                                                               <  3:123284/149‑1 (MQ=255)                       ggCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTggg                                                                                                                      >  4:218922/1‑148 (MQ=255)                         cTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCtt                                                                                                                         <  3:90790/143‑1 (MQ=255)                         cTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCtt                                                                                                                         >  4:90790/1‑143 (MQ=255)                            ggACCGAACAGGGCATAGGCCAGCAGGGGGATGGTGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGCTCACCACCATGctcgt                                                                                                                                     >  4:295056/1‑124 (MQ=255)                            ggACCGAACAGGGCATAGGCCAGCAGGGGGATGGTGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGCTCACCACCATGctcgt                                                                                                                                     <  3:295056/128‑5 (MQ=255)                             gACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAg                                                                                                               >  8:197997/1‑149 (MQ=255)                                    cAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCt                                                                                                        <  3:112497/149‑1 (MQ=255)                                     aGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAAGCTGCCtt                                                                                                       >  6:317845/1‑149 (MQ=255)                                       ggCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGgcggcg                                                                                                                                 <  4:272445/121‑1 (MQ=255)                                       ggCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGgcggcg                                                                                                                                 >  3:272445/1‑121 (MQ=255)                                       ggCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCCGCGATGTCCTTGGGTAGCCAGGCTGCCTTg                                                                                                      >  7:277418/1‑148 (MQ=255)                                         cATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAg                                                                                                   <  1:21935/149‑1 (MQ=255)                                            aGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGCTCACCACCATGctc                                                                                                                                       >  1:128444/1‑108 (MQ=255)                                            aGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGCTCACCACCATGctc                                                                                                                                       <  2:128444/110‑3 (MQ=255)                                                  gcagGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAgg                                                                                          >  6:237875/1‑149 (MQ=255)                                                        ggATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGAc                                                                                     >  5:195124/1‑148 (MQ=255)                                                                 tAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCggg                                                                            >  3:192480/1‑148 (MQ=255)                                                                  aGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGcc                                                                                                                 <  5:16250/110‑1 (MQ=255)                                                                  aGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGcc                                                                                                                 >  6:16250/1‑110 (MQ=255)                                                                    cGTACAGGCCATAGATACCGGGGACGCCTGAAGC‑CTCGCCGTAGTCGATGGCCACCGGCACGAGCATGTTGGCGAGCACCAGGCCGGCGTCGATGGCCTTGGGTAGCCTGGCTGCCATGTAGTGGAGCAGGGTGACCAGGCCGGGTAAc                                                                        <  4:302156/149‑1 (MQ=255)                                                                    cGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAg                                                                                                                                           >  4:147693/1‑82 (MQ=255)                                                                    cGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAg                                                                                                                                           <  3:147693/82‑1 (MQ=255)                                                                        cAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCa                                                                                            <  3:260990/125‑1 (MQ=255)                                                                        cAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCa                                                                                            >  4:260990/1‑125 (MQ=255)                                                                           gCCATAGATACCGGGGACGCCGGAAGCGCTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGt                                                                  >  1:16088/1‑149 (MQ=255)                                                                            ccATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTg                                                                <  2:16088/149‑1 (MQ=255)                                                                             cATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGc                                                               <  8:47688/149‑1 (MQ=255)                                                                                 gATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGTTAGCCAAGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGt                                                           <  5:317845/149‑1 (MQ=255)                                                                                  aTACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTc                                                          <  7:197997/149‑1 (MQ=255)                                                                                    aCCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAgg                                                                                          >  6:109291/1‑115 (MQ=255)                                                                                    aCCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAgg                                                                                          <  5:109291/115‑1 (MQ=255)                                                                                         ggACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGa                                                   <  3:191122/149‑1 (MQ=255)                                                                                                     c‑cTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTgggg                                        >  6:146774/1‑148 (MQ=255)                                                                                                          ggCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTg                                                                                                 >  2:222059/1‑87 (MQ=255)                                                                                                          ggCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTg                                                                                                 <  1:222059/87‑1 (MQ=255)                                                                                                             gTAGGCGACGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGgtgt                                >  7:167245/1‑149 (MQ=255)                                                                                                                ggCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCgg                             >  2:165958/1‑149 (MQ=255)                                                                                                                ggCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCgg                             >  3:23376/1‑149 (MQ=255)                                                                                                                 gCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGc                            >  7:74995/1‑149 (MQ=255)                                                                                                                 gCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGc                            >  3:330057/1‑149 (MQ=255)                                                                                                                  cGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGct                           <  5:146774/149‑1 (MQ=255)                                                                                                                    acgCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGctct                         <  8:167245/147‑1 (MQ=255)                                                                                                                      gCCCACAGGCACGAGCATGGTGGTGAGCACGAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTggg                                              <  3:248809/126‑1 (MQ=255)                                                                                                                      gCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTggg                                              >  4:248809/1‑126 (MQ=255)                                                                                                                       cccACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTgg                       <  6:195124/148‑1 (MQ=255)                                                                                                                        ccACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAgg                                                                                          >  2:257706/1‑80 (MQ=255)                                                                                                                        ccACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAgg                                                                                          <  1:257706/80‑1 (MQ=255)                                                                                                                        ccACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGc                     <  4:330057/149‑1 (MQ=255)                                                                                                                        ccACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGc                     >  3:146329/1‑149 (MQ=255)                                                                                                                           caGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGCgag                  <  8:74995/149‑1 (MQ=255)                                                                                                                              gCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCTTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGCGAGAAg               <  5:237875/149‑1 (MQ=255)                                                                                                                                 cGAGCATGGTGTTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGCGAGAAGGgc            <  6:231164/149‑1 (MQ=255)                                                                                                                                           ggtgAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGCGAGACGGGCGCTGAACCTc  <  4:192480/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                  TCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGCGTACAGGCCATAGATACCGGGGACGCCGGAAGC‑CTCGGCGTAGGCGATGCCCACAGGCACGAGCATGGTGGTGAGCACCAGGCCGGCGGCGATGTCCTTGGGTAGCCAGGCTGCCTTGTAGTGGAGCAGGGTGACCAGGCCGGGTAACCAGCGTTGCCAGTCGAAACGAGTGGGTTGGGGCATGGTGTCGGCTCTGGGCGAGAAGGGCGCTGAACCTC  >  NC_002947/835828‑836113

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