Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,297,595 T→C 8.9% A153A (GCA→GCG)  PP_3768 ← shikimate 5‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,297,5950TC 8.9% 152.4 / 4.8 56A153A (GCA→GCG) PP_3768shikimate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (30/21);  new base C (3/2);  total (33/23)
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.46e-01

CGGGAAACGCCAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGCCAGGTAATCGGT  >  NC_002947/4297447‑4297741
                                                                                                                                                    |                                                                                                                                                  
cGGGAAACGCCAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGt                                                                                                                                                    >  4:104211/1‑149 (MQ=255)
 gggAAACGCCAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGt                                                                                                                                                    <  3:104211/148‑1 (MQ=255)
        gCCAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGAACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACgat                                                                                                                                            >  7:305038/1‑149 (MQ=255)
         ccAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAATTTTCGCCTGGCTCAGGCACCCAGCGATACCAACAGCAATCAGCCAACTGACCCCCGGCCTGCACGTTGCACGCAACgag                                                                                                                                            >  4:235964/1‑147 (MQ=255)
         ccAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCACGCAACgagc                                                                                                                                           <  3:235964/149‑3 (MQ=255)
          cAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACgatga                                                                                                                                          >  1:15917/1‑149 (MQ=255)
           aGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACgatgat                                                                                                                                         >  1:152532/1‑149 (MQ=255)
              tcttctGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGcc                                                                                                                                      >  6:102044/1‑149 (MQ=255)
                ttctGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCtt                                                                                                                                    >  5:317052/1‑149 (MQ=255)
                     ccTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCgg                                                                                                                                 <  6:317052/147‑1 (MQ=255)
                     ccTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCgg                                                                                                                                 <  8:121477/147‑1 (MQ=255)
                       tCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCaaa                                                                                                                             >  1:87453/1‑149 (MQ=255)
                        ccAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAAc                                                                                                                            >  1:330562/1‑149 (MQ=255)
                            cccGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCgg                                                                                                                        <  2:15917/149‑1 (MQ=255)
                                        aCCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAgg                                                                                                            >  7:288048/1‑149 (MQ=255)
                                                gCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGcc                                                                                                    >  3:39112/1‑149 (MQ=255)
                                                gCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGcc                                                                                                    <  6:212101/149‑1 (MQ=255)
                                                     aCGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCACGCAACgagca                                                                                                                                          <  4:218651/106‑4 (MQ=255)
                                                     aCGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCACGCAACgagca                                                                                                                                          >  3:218651/1‑103 (MQ=255)
                                                          cACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTg                                                                                           <  2:152532/148‑1 (MQ=255)
                                                            ccAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGc                                                                                         >  3:253268/1‑148 (MQ=255)
                                                                     gcgggcAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGc                                                                                               <  3:222314/133‑1 (MQ=255)
                                                                     gcgggcAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGc                                                                                               >  4:222314/1‑133 (MQ=255)
                                                                        ggcAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTg                                                                            >  2:219520/1‑149 (MQ=255)
                                                                         gcAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGc                                                                           <  2:330562/149‑1 (MQ=255)
                                                                           aaaTGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCt                                                                         >  5:257829/1‑149 (MQ=255)
                                                                              tGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCg                                                                       >  7:288253/1‑148 (MQ=255)
                                                                                       gCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCtt                                                                                            >  4:191009/1‑118 (MQ=255)
                                                                                       gCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCtt                                                                                            <  3:191009/118‑1 (MQ=255)
                                                                                           aGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCg                                                         >  8:154246/1‑149 (MQ=255)
                                                                                               aCCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGcc                                                     >  4:262482/1‑149 (MQ=255)
                                                                                               aCCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGcc                                                     >  2:152907/1‑149 (MQ=255)
                                                                                               aCCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGcc                                                     >  5:107168/1‑149 (MQ=255)
                                                                                                cccAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGccc                                                    <  7:154246/149‑1 (MQ=255)
                                                                                                cccAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGccc                                                    <  6:257829/149‑1 (MQ=255)
                                                                                                   aGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGa                                                 <  1:219520/149‑1 (MQ=255)
                                                                                                   aGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGa                                                 <  8:288253/149‑1 (MQ=255)
                                                                                                           ccccaGACCTCCTCCAACTGACGCCCGGCCGGCTCGTGGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTg                                         <  8:305038/146‑1 (MQ=255)
                                                                                                            cacaGACATCCGCCAACTGCCGCCCGGCCTGCTCGTTGCCTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGa                                        <  8:289823/149‑1 (MQ=255)
                                                                                                            cacaGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGc                                                           >  4:127163/1‑130 (MQ=255)
                                                                                                            cacaGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGc                                                           <  3:127163/130‑1 (MQ=255)
                                                                                                            cacaGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGa                                        >  2:173761/1‑149 (MQ=255)
                                                                                                            cacaGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGa                                        >  7:77339/1‑149 (MQ=255)
                                                                                                              caGACAACCGCCAACTGAAGCCCGGCCTGCTCGATGATTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGACTTCGCCAAGCCACCGCTTACTAGCAGGGCAAAGTCCGAGCAAGGATCGACCTGAt                                       >  2:320094/1‑148 (MQ=255)
                                                                                                               aGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGc                                     <  8:288048/149‑1 (MQ=255)
                                                                                                                   aTCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAg                                 >  7:155240/1‑149 (MQ=255)
                                                                                                                          aaCTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGAc                                            <  6:12190/131‑1 (MQ=255)
                                                                                                                          aaCTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGAc                                            >  5:12190/1‑131 (MQ=255)
                                                                                                                            cTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAAtt                        <  4:253268/149‑1 (MQ=255)
                                                                                                                               aCGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTg                       <  8:77339/147‑1 (MQ=255)
                                                                                                                                 gCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGgcg                    >  1:212409/1‑148 (MQ=255)
                                                                                                                                 gCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGgcg                    >  8:267706/1‑148 (MQ=255)
                                                                                                                                 gCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGgcg                    >  4:86690/1‑148 (MQ=255)
                                                                                                                                 gCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGgcg                    >  5:225559/1‑148 (MQ=255)
                                                                                                                                   ccGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGACACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCAc                 >  5:153269/1‑149 (MQ=255)
                                                                                                                                     ggCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACgg               <  3:321466/149‑1 (MQ=255)
                                                                                                                                        atcatCGTTGCACGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGc              >  8:61716/5‑147 (MQ=255)
                                                                                                                                                  cGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGCCAGGTAATCGgt  <  7:267706/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGGGAAACGCCAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGCCAGGTAATCGGT  >  NC_002947/4297447‑4297741

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