Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,103,473	A→T	6.9%	A220A (GCA→GCT)  M1M (ATG→TTG) †	PP_2713 → PP_2714 →	two‑component system response regulator sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,103,473	0	A	T	6.9%	170.8 / 3.2	58	A220A (GCA→GCT)  M1M (ATG→TTG) †	PP_2713 PP_2714	two‑component system response regulator sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base T (2/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

GAGCGCCTGACTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGAACGCGCGTTTGC  >  NC_002947/3103327‑3103620                                                                                                                                                   |                                                                                                                                                    gAGCGCCTGACTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCAt                                                                                                                                                    >  2:289488/1‑148 (MQ=255)  aGCGCCTGACTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGa                                                                                                                                                  <  1:7352/149‑1 (MQ=255)      ccTGACTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTc                                                                                                                                              <  1:289488/149‑1 (MQ=255)        tGACTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGc                                                                                                                                            >  3:266875/1‑149 (MQ=255)          aCTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCATGCCTTGACCTCgac                                                                                                                                           <  7:174388/148‑3 (MQ=255)          aCTCAACTGCTCTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCATGCCTTGACCTCga                                                                                                                                            <  6:307457/147‑2 (MQ=255)               aCTGCTGTATGGCTGGCTTGAAGAGCCGGAGAGCAAGACGCTGGAGGTGAAGATCTACCATTTGCCCATGATGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGGCGAGGTCAAGGGATGACCTCGCTTCGCCa                                                                                                                                     <  2:93644/149‑1 (MQ=255)                 tgctgTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAgc                                                                                                                                   >  7:284428/1‑149 (MQ=255)                    tgTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCa                                                                                                                                >  8:121105/1‑149 (MQ=255)                    tgTATGCCTGGGATGAAGAGCCGGAGAGCAACACGCTTGAGGTAAACATCTACCATTTGCTAAAGAAGCTTTTCAACGGACTGATTCGCACGGTGCGCGGAAATCGATATCTGTTAGAGGTCAAGGCATGACCTAACTTAGCCAgagaa                                                                                                                                >  8:246316/1‑145 (MQ=255)                          gCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCtgt                                                                                                                          <  7:121105/149‑1 (MQ=255)                           cTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCtgtg                                                                                                                         >  2:86375/1‑149 (MQ=255)                            tGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTgg                                                                                                                        >  4:279609/1‑149 (MQ=255)                              ggATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGgcg                                                                                                                      >  1:118458/1‑149 (MQ=255)                                 tGAAGAGCCGGAGAGCAACAAGCTGCAGCTGACCAGCTACCATTTGCGCAAGAAGCTGTTCAACCGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGt                                                                                                                   <  4:266875/149‑1 (MQ=255)                                        ccGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGctc                                                                                                            >  4:340568/1‑149 (MQ=255)                                         cGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTg                                                                                                                  <  5:185438/142‑1 (MQ=255)                                         cGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTg                                                                                                                  >  6:185438/1‑142 (MQ=255)                                             gagCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTtgct                                                                                                       <  2:225764/149‑1 (MQ=255)                                              agCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTtgctg                                                                                                      <  4:330151/149‑1 (MQ=255)                                               gCAACACGCTGGAGGTGAACATCTACCATTTGCGAAAGAAGCTTTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTtgctgc                                                                                                     <  7:246316/149‑1 (MQ=255)                                                    acGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTc                                                                                                <  8:284428/149‑1 (MQ=255)                                                    acGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTc                                                                                                <  3:340568/149‑1 (MQ=255)                                                      gCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGGGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCgg                                                                                              <  4:61070/149‑1 (MQ=255)                                                             gTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGc                                                                                       <  3:279609/149‑1 (MQ=255)                                                                      tACCATTTGCGCGAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTcgcg                                                                              >  2:258298/1‑149 (MQ=255)                                                                      tACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTcgcg                                                                              >  8:222418/1‑149 (MQ=255)                                                                       aCCATTTGCGCAAGAAGCTGTTCAACTGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGc                                                                                       <  4:151932/139‑1 (MQ=255)                                                                       aCCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTcgcg                                                                              <  3:83912/148‑1 (MQ=255)                                                                       aCCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTcgcg                                                                              >  4:83912/1‑148 (MQ=255)                                                                       aCCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCgcgc                                                                             >  1:315214/1‑149 (MQ=255)                                                                       aCCATTTGCGCAAGAAGCTGTTCAACCGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGc                                                                                       >  3:151932/1‑139 (MQ=255)                                                                        ccATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCgcgc                                                                             <  1:47056/148‑1 (MQ=255)                                                                                gcAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAAc                                                                     >  8:221423/1‑148 (MQ=255)                                                                                   agaagCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTAcc                                                                 <  6:59085/149‑1 (MQ=255)                                                                                     cagCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCAc                                                               <  4:34952/148‑1 (MQ=255)                                                                                     aagCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCa                                                                <  2:325715/148‑1 (MQ=255)                                                                                     aagCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCa                                                                >  1:325715/1‑148 (MQ=255)                                                                                        cTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGACGCTCtt                                                                                                          <  2:8098/103‑1 (MQ=255)                                                                                        cTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGAAGCTCtt                                                                                                          >  1:8098/1‑103 (MQ=255)                                                                                                          cGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTaa                                           <  7:221423/148‑1 (MQ=255)                                                                                                                 tGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCGCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCCCGACAGCCGCCACGAAATTAACGAGGta                                    >  6:46220/1‑148 (MQ=255)                                                                                                                          ttGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGAcgcg                           >  5:125239/1‑148 (MQ=255)                                                                                                                             gCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATaa                                                                      >  8:130673/1‑102 (MQ=255)                                                                                                                             gCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATaa                                                                      <  7:130673/102‑1 (MQ=255)                                                                                                                               tATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGa                                        >  6:150128/1‑130 (MQ=255)                                                                                                                               tATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGa                                        <  5:150128/130‑1 (MQ=255)                                                                                                                               tATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTg                     >  4:156089/1‑149 (MQ=255)                                                                                                                               tATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGAACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTg                     >  8:235265/1‑149 (MQ=255)                                                                                                                                 tCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGc                   >  3:275329/1‑149 (MQ=255)                                                                                                                                      gcGAGGTCATGCCTTGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAGCGAGGTATATGACGCGCACCTGGCCCAg               >  2:115004/2‑148 (MQ=255)                                                                                                                                       cGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGa              <  7:256463/148‑1 (MQ=255)                                                                                                                                       cGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGa              >  8:256463/1‑148 (MQ=255)                                                                                                                                        gAGGTCATGCCTTGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGAAc            >  8:63849/1‑149 (MQ=255)                                                                                                                                        gAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGaa             <  2:118458/148‑1 (MQ=255)                                                                                                                                           gTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGAAcgcg         >  5:274331/1‑149 (MQ=255)                                                                                                                                              aaGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGAACGCGCGt      >  7:138311/1‑149 (MQ=255)                                                                                                                                                  cATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGAACGCGCGTTtgc  <  1:86375/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GAGCGCCTGACTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGAACGCGCGTTTGC  >  NC_002947/3103327‑3103620

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑