Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 6,022,319 A→G 6.4% L326P (CTG→CCG)  puuC_2 ← NADP/NAD‑dependent aldehyde dehydrogenase PuuC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24406,022,3190AG 6.4% 193.1 / 2.4 62L326P (CTG→CCG) puuC_2NADP/NAD‑dependent aldehyde dehydrogenase PuuC
Reads supporting (aligned to +/- strand):  ref base A (30/26);  new base G (2/2);  total (32/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 = 1.00e+00

CCGGCCAGCAGCTTGGCGCCGTCCTTGTGGCCAGCCTCGATGTACGACAGCACGGTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTGGTTGAAGGCGATGGCGCTGGCTGCGGCCTCGGCGGCTGCTT  >  KT2440/6022187‑6022445
                                                                                                                                    |                                                                                                                              
ccgGCCAGCAGCTTGGCGCCGTCCTTGTGGCCAGCCTCGATGTACGACAGCACGGTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAgg                                                                                                                              >  7:79507/1‑135 (MQ=255)
           cTTGGCGCCGGCCTTGTGGCCAGCCTCGATGTACGACAGCACGGTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCTGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTaccac                                                                                                                   <  7:227996/135‑1 (MQ=255)
               gcgcCGTCCTTGTGGCCAGCCTCGATGTACGACAGCACGGTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATc                                                                                                               <  5:199582/135‑1 (MQ=255)
                cgcCGTCCTTGTGGCCAGCCTCGATGTACGACAGCACGGTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCg                                                                                                              <  3:142088/135‑1 (MQ=255)
                     tCCTTGTGGCCAGCCTCGATGTACGACAGCACGGTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAgg                                                                                                         <  2:110868/135‑1 (MQ=255)
                        ttGTGGCCAGCCTCGATGTACGACAGCACGGTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAAc                                                                                                      >  2:1084/1‑135 (MQ=255)
                                     cGATGTACGACAGCACGGTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCCGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCg                                                                                         >  8:119619/1‑135 (MQ=255)
                                     cGATGTACGACAGCACGGTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATc                                                                                          >  6:100943/1‑134 (MQ=255)
                                             gACAGCACGGTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTc                                                                                  <  2:171317/134‑1 (MQ=255)
                                                   aCGGTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACca                                                                             >  7:5845/1‑133 (MQ=255)
                                                     ggTGTTCATCTGCTGGGTTTCCCCCAAGGCCCCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACcagcag                                                                         <  8:238256/135‑1 (MQ=255)
                                                     ggTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACcagcag                                                                         >  8:49154/1‑135 (MQ=255)
                                                     ggTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACcagcag                                                                         <  1:170746/135‑1 (MQ=255)
                                                     ggTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCGCTCCCACCATCGGCAGGACCTTGTCCTTGCTCGAACGCTCGACcagccg                                                                         >  5:18376/1‑133 (MQ=255)
                                                     ggTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGGCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACcagcag                                                                         <  1:109921/135‑1 (MQ=255)
                                                        gtTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCg                                                                      <  3:278975/135‑1 (MQ=255)
                                                          tCATCTGCTGGGTATCGCCCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCGACCACCATCGGCAGGAACTTGTCCTTGATCGACCGCACGCCCAGCAGGCg                                                                      >  7:74063/1‑133 (MQ=255)
                                                          tCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCggg                                                                    <  8:201389/135‑1 (MQ=255)
                                                          tCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCggg                                                                    <  8:10365/135‑1 (MQ=255)
                                                          tCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCg                                                                      >  2:305528/1‑133 (MQ=255)
                                                          tCATCTGCTGGGTATCGACCAAGGCACCGACAGGGGTCTGCGGGTCCAGCGGGGTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCggg                                                                    <  7:267797/135‑1 (MQ=255)
                                                            ctCTGCTGGGTAGCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAg                                                                  <  8:166959/134‑1 (MQ=255)
                                                             tCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGc                                                                 <  1:1084/135‑1 (MQ=255)
                                                             tCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGc                                                                 <  4:160144/135‑1 (MQ=255)
                                                              ctggtgGGTATCGTACATTGCTCCGTGAGTTGTCTGCGGGTCCAGCGGGTTGGCTGGCTTCCAGCCGTTCAGGGCCTCTTCCACCATCGGCTGGAACTTGTCCTTGTTCGATGGCTGGTCCAGCAGGCTGGatc                                                                 <  6:327924/134‑3 (MQ=255)
                                                               tgctgGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCt                                                               >  8:98152/1‑135 (MQ=255)
                                                                    ggTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCGGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACcagc                                                                           <  7:180388/118‑1 (MQ=255)
                                                                    ggTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCGGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACcagc                                                                           >  8:180390/1‑118 (MQ=255)
                                                                    ggTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCGGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACcagc                                                                           >  8:180388/1‑118 (MQ=255)
                                                                    ggTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCGGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACcagc                                                                           <  7:180390/118‑1 (MQ=255)
                                                                    ggTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGt                                                          >  3:303266/1‑135 (MQ=255)
                                                                    ggTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGt                                                          >  7:204188/1‑135 (MQ=255)
                                                                      tATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTACTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGc                                                        >  7:17686/1‑135 (MQ=255)
                                                                         cGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGcaca                                                     <  7:98152/135‑1 (MQ=255)
                                                                         cGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGcaca                                                     <  4:331008/135‑1 (MQ=255)
                                                                         cGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTACTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGcac                                                      >  5:233428/1‑134 (MQ=255)
                                                                            ccAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGCCAGGAACTTGTACTTGATCGAGCGCTCGACCACCCGGCGGGAGCCTGCAGTGCACACtt                                                  >  8:88450/1‑135 (MQ=255)
                                                                             caagtaaCCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCCGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTc                                                 <  2:301090/129‑1 (MQ=255)
                                                                             cAAGGCACCGACCGTGGGCTGCGGGTCCAGGGGGTTGCCTGGCTTCCAGCCCTTCCGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCCGCAGGCGGGAGCCTGCAGTGCACACTTc                                                 <  8:128366/135‑1 (MQ=255)
                                                                             cAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTc                                                 <  4:294784/135‑1 (MQ=255)
                                                                              aaGGCACCGACAGTGTTCTGCGGGTCCAGCGGGTTGCCTGGCTGCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTTTCCTTGATCGAACGCTCGACCAGCAGGCGGTAGCCTGCAGTGCACACTTCg                                                <  3:97732/135‑1 (MQ=255)
                                                                              aaGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGCGCCTGCAGTGCACACTTCg                                                >  1:96170/1‑135 (MQ=255)
                                                                                  cACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACATCGCCCt                                            >  5:200834/1‑135 (MQ=255)
                                                                                   aCCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTg                                           >  7:94871/1‑135 (MQ=255)
                                                                                   aCCGACAGTGGTCTGCGGGTCAAGCTAGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTg                                           >  8:308339/1‑135 (MQ=255)
                                                                                    ccGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTgg                                          <  4:165884/135‑1 (MQ=255)
                                                                                        cAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTGGtt                                        >  4:122617/1‑133 (MQ=255)
                                                                                        cAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTGGtt                                        >  1:92894/1‑133 (MQ=255)
                                                                                              tCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTGGTTGAAGGCGa                                >  2:88270/1‑135 (MQ=255)
                                                                                                tGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTGGTTGAAGGCGCTg                              >  5:262090/1‑135 (MQ=255)
                                                                                                   gggTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGGTCGCCCAGCAGGCCGGAGCCTCCAGTGCACACTTCGCCCTCGGAGAAGGCGATGgcg                           >  1:161180/1‑135 (MQ=255)
                                                                                                   gggTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTGGTTGAAGGCGATGgcg                           >  1:204014/1‑135 (MQ=255)
                                                                                                      tCCAGCGGGTGGCCTGGCTTCCAGCCCTTCAGGGCCGCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTGGTTGAAGGCGATGGCgctg                        <  1:97623/135‑1 (MQ=255)
                                                                                                          gCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCCGTGCACACTTCGCCCTGGTTGAAGGCGATGGCgctggctg                    >  4:308118/1‑135 (MQ=255)
                                                                                                                     ggCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGca                                                       >  6:136709/1‑89 (MQ=255)
                                                                                                                     ggCTTCCAGCCCTTCAGGGCATCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGca                                                       <  5:136709/89‑1 (MQ=255)
                                                                                                                       cTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTGGTTGAAGGCGATGGCGCTGGCTGCGGCCTcggcg         <  6:200834/133‑1 (MQ=255)
                                                                                                                       cTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTGGTTGAAGGCGATGGCGCTGGCTGCGGCCTcggcg         <  8:79507/133‑1 (MQ=255)
                                                                                                                        ttCCAGCCATTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGCCCGATCGCCCAGCAGTAGGGAGCCTGCAGTGCACACTTCGCCCTGGTTGAAGTCGATGGAGATGGCTGCGGCCTCGGCGgct      >  6:145664/1‑135 (MQ=255)
                                                                                                                          ccAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTGGTTGAAGGCGATGGCGCTGGCTGCGGCCTCGGCGgctgc    >  7:196201/1‑135 (MQ=255)
                                                                                                                            aGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTGGTTGAAGGCGATGGCGCTGGCTGCGGCCTCGGCGgctgct   <  4:303266/134‑1 (MQ=255)
                                                                                                                             gCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTGGTTGAAGGCGATGGCGCTGGCTGCGGCCTCGGCGGCTGCtt  >  6:20241/1‑134 (MQ=255)
                                                                                                                                    |                                                                                                                              
CCGGCCAGCAGCTTGGCGCCGTCCTTGTGGCCAGCCTCGATGTACGACAGCACGGTGTTCATCTGCTGGGTATCGACCAAGGCACCGACAGTGGTCTGCGGGTCCAGCGGGTTGCCTGGCTTCCAGCCCTTCAGGGCCTCTACCACCATCGGCAGGAACTTGTCCTTGATCGAACGCTCGACCAGCAGGCGGGAGCCTGCAGTGCACACTTCGCCCTGGTTGAAGGCGATGGCGCTGGCTGCGGCCTCGGCGGCTGCTT  >  KT2440/6022187‑6022445

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