Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,142,680 C→T 8.9% I271I (ATC→ATT)  PP_3646 → aldehyde dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,142,6800CT 8.9% 158.0 / 4.5 56I271I (ATC→ATT) PP_3646aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base C (31/20);  new base T (2/3);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

CGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGAC  >  NC_002947/4142537‑4142827
                                                                                                                                               |                                                                                                                                                    
cgcAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATg                                                                                                                                                  <  1:19590/148‑1 (MQ=255)
   aaGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATTGATGgcgc                                                                                                                                              <  4:71317/149‑1 (MQ=255)
            ttCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCgg                                                                                                                                     >  8:82914/1‑149 (MQ=255)
              cACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCgg                                                                                                                                     <  6:224948/147‑1 (MQ=255)
                 cggcggTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATc                                                                                                                                 >  6:59412/1‑148 (MQ=255)
                 cggcggTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATc                                                                                                                                 >  6:185105/1‑148 (MQ=255)
                      gTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCACTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGcc                                                                                                                           <  8:76033/149‑1 (MQ=255)
                           gccACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTc                                                                                                                      >  2:189711/1‑149 (MQ=255)
                            ccACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTcc                                                                                                                     >  7:7323/1‑149 (MQ=255)
                                   ccgccACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCa                                                                                                                >  1:74836/1‑147 (MQ=255)
                                        aCGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCAGCGTCCGGGCAAAGCTGc                                                                                                         >  4:98469/1‑149 (MQ=255)
                                             gTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGc                                                                                                    <  2:137722/149‑1 (MQ=255)
                                                cGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCGAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTCCGTGGCCgg                                                                                                 >  3:51753/1‑149 (MQ=255)
                                                 gcagcaGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGc                                                                                                >  6:155399/1‑149 (MQ=255)
                                                     cagcGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGGCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTcgc                                                                                            <  3:233445/149‑1 (MQ=255)
                                                                     gCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGa                                                                            >  2:320848/1‑149 (MQ=255)
                                                                     gCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGa                                                                            <  3:341870/149‑1 (MQ=255)
                                                                       cAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGAt                                                                           >  3:38687/1‑148 (MQ=255)
                                                                            tGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATc                                                                     >  2:128768/1‑149 (MQ=255)
                                                                                 cTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTcga                                                                <  8:7323/149‑1 (MQ=255)
                                                                                    gAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTcgac                                                               <  2:133271/147‑1 (MQ=255)
                                                                                    gAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTcgac                                                               >  8:52950/1‑147 (MQ=255)
                                                                                      aCTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGt                                                           <  5:155399/149‑1 (MQ=255)
                                                                                      aCTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGt                                                           >  4:332229/1‑149 (MQ=255)
                                                                                       cTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGtt                                                          <  5:161398/149‑1 (MQ=255)
                                                                                          ggcggcAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTt‑c                                                        <  1:189711/147‑1 (MQ=255)
                                                                                           gcggcAAGTCGCCGAACATCATCTTCGCCGACGCTGACCTGGACAGCGCCATCAATGGCGCCGTGGCCTGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTt‑cgtc                                                     <  1:160998/149‑1 (MQ=255)
                                                                                            cggcAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTt‑cgtc                                                     >  4:160244/1‑148 (MQ=255)
                                                                                             ggcAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTt‑cgtc                                                     >  3:150715/1‑147 (MQ=255)
                                                                                               cAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAgc                                                 <  5:185105/149‑1 (MQ=255)
                                                                                                  gTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGcc                                              >  8:150910/1‑149 (MQ=255)
                                                                                                  gTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGACATCAATGGCGACGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGcc                                              >  5:181270/1‑149 (MQ=255)
                                                                                                   tCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCt                                             >  6:114843/1‑149 (MQ=255)
                                                                                                      ccGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGAt                                          >  6:115030/1‑149 (MQ=255)
                                                                                                       cGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGAtt                                         >  5:242530/1‑149 (MQ=255)
                                                                                                       cGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGAtt                                         >  8:368975/1‑149 (MQ=255)
                                                                                                             tcatcTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCAGCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGCGTT‑CGTCGAGCGCCTGATTGCCcgc                                   <  5:272059/149‑1 (MQ=255)
                                                                                                              catcTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCcgcg                                  >  7:170300/1‑149 (MQ=255)
                                                                                                               atcTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCgcgc                                 >  2:306015/1‑149 (MQ=255)
                                                                                                                 cTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGcc                                >  6:174989/1‑148 (MQ=255)
                                                                                                                     gccgacgccgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAgcg                           >  7:114420/1‑149 (MQ=255)
                                                                                                                     gccgacgccgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAgc                            <  5:174989/148‑1 (MQ=255)
                                                                                                                       cgacgccgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGcgca                         >  3:248199/1‑149 (MQ=255)
                                                                                                                          cgccgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatcc                      <  6:242530/149‑1 (MQ=255)
                                                                                                                              gacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatccgcat                  <  4:248199/149‑1 (MQ=255)
                                                                                                                               acCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatccgcatc                 <  1:320848/149‑1 (MQ=255)
                                                                                                                                 cTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCgg               <  1:168355/149‑1 (MQ=255)
                                                                                                                                 cTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCgg               >  3:115936/1‑149 (MQ=255)
                                                                                                                                       ggcgcCATTGATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatc                       <  1:252001/135‑1 (MQ=255)
                                                                                                                                       ggcgcCATTGATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatc                       <  1:246568/135‑1 (MQ=255)
                                                                                                                                       ggcgcCATTGATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatc                       >  2:252001/2‑135 (MQ=255)
                                                                                                                                       ggcgcCATTGATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatc                       >  2:246568/2‑135 (MQ=255)
                                                                                                                                       aGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAAccc         <  3:160244/149‑1 (MQ=255)
                                                                                                                                         cgcCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCg        >  7:323983/1‑148 (MQ=255)
                                                                                                                                         cgcCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCg        >  6:187811/1‑148 (MQ=255)
                                                                                                                                              tCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGgac  >  6:158682/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTT‑CGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGAC  >  NC_002947/4142537‑4142827

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