Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,279,197 A→C R145S (AGA→AGC agaA → putative N‑acetylgalactosamine‑6‑phosphate deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,279,1970AC99.4% 435.4 / NA 165R145S (AGA→AGCagaAputative N‑acetylgalactosamine‑6‑phosphate deacetylase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (67/96);  total (68/97)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

CGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGAGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAG  >  NC_000913/3279124‑3279269
                                                                         |                                                                        
cGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGt                                                                        >  1:3477568/1‑76 (MQ=255)
 gTCTCTGCATCGTGCGCGATTGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTc                                                                       <  1:5414069/76‑1 (MQ=255)
 gTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTc                                                                       >  2:5429651/1‑76 (MQ=255)
 gTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTc                                                                       <  1:3025789/76‑1 (MQ=255)
 gTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTc                                                                       >  1:3253538/1‑76 (MQ=255)
 gTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTc                                                                       >  1:2671132/1‑76 (MQ=255)
 gTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTc                                                                       >  1:3836685/1‑76 (MQ=255)
 gTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTc                                                                       <  1:89589/76‑1 (MQ=255)
 gTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTc                                                                       >  2:2928461/1‑76 (MQ=255)
 gTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTc                                                                       >  2:3184477/1‑76 (MQ=255)
 gTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTc                                                                       >  2:5100407/1‑76 (MQ=255)
  tCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCg                                                                      <  1:5497898/76‑1 (MQ=255)
   cgatGCATCCGGCGCGCATGCTGGGTGTTGATGGTGTTCTGGGAGCGATTAAACCGGGCCAACGCGCCAGcgtcgt                                                                     >  1:4931495/4‑76 (MQ=255)
    gCTGCATCCGGCGCGAATGCTGGGTGATGATGGTGTTCTGGGTGCGCTTAAACGGGGCACACGCGCGAGAGTCCtt                                                                    >  1:5639484/1‑76 (MQ=255)
      tGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTgc                                                                  >  1:3801678/1‑76 (MQ=255)
        cATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTgcgc                                                                >  1:1654284/1‑76 (MQ=255)
         aTCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCt                                                               <  2:4753236/76‑1 (MQ=255)
         aTCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCt                                                               >  1:3275899/1‑76 (MQ=255)
          tCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTCgcgccg                                                              >  2:838764/1‑74 (MQ=255)
           ccGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTgg                                                             >  1:929979/1‑76 (MQ=255)
            cGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGa                                                            <  1:5139855/76‑1 (MQ=255)
            cGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGa                                                            <  1:279167/76‑1 (MQ=255)
             ggCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGAt                                                           <  1:434838/76‑1 (MQ=255)
             ggCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGAt                                                           >  2:3420016/1‑76 (MQ=255)
              gcgcGAATGATGGGTGTTGATGTTGGTCTGTGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATa                                                          <  1:3714540/76‑1 (MQ=255)
               cgcgAATGCTGGGTGTTGATGGTGTTCTGGGATCTCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAg                                                         <  1:3707636/76‑1 (MQ=255)
               cgcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAg                                                         <  1:3593969/76‑1 (MQ=255)
                tgtaaTGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACAGGGCAAACTCGCCAGCGTCGTTGCGCTGGATAGc                                                        <  1:1674981/73‑1 (MQ=255)
                gcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGc                                                        <  1:3686834/76‑1 (MQ=255)
                gcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGc                                                        <  2:1209479/76‑1 (MQ=255)
                gcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGGCGTTGCGCTGGATAGc                                                        >  2:1871400/1‑76 (MQ=255)
                gcgAATGCTGGGCGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGc                                                        >  2:1627998/1‑76 (MQ=255)
                 cgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGCATAGCg                                                       >  2:1427422/1‑76 (MQ=255)
                  gAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCgg                                                      <  2:4693082/76‑1 (MQ=255)
                    aTGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGc                                                    >  2:3885449/1‑76 (MQ=255)
                    aTGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGc                                                    <  2:1680500/76‑1 (MQ=255)
                    aTGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGc                                                    <  1:1078500/76‑1 (MQ=255)
                      gCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTa                                                  >  2:3544001/1‑76 (MQ=255)
                      gCTGGGTGTTGATGGTGTACTGGTATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTa                                                  <  2:4188872/76‑1 (MQ=255)
                      gCTGGGTGTTGATGGGGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTa                                                  >  2:4664713/1‑76 (MQ=255)
                       cTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTAc                                                 <  2:4145309/76‑1 (MQ=255)
                       cTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTAc                                                 <  1:3719832/76‑1 (MQ=255)
                         gggTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACAt                                               <  2:5590522/76‑1 (MQ=255)
                          ggTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACAtg                                              >  1:5567163/1‑76 (MQ=255)
                             gtTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGc                                           >  1:2611957/1‑76 (MQ=255)
                               tGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTCGCGCTGGATAGCGGGCTACATGTGcaa                                         >  2:3033058/1‑76 (MQ=255)
                                gATGGTGTGCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaac                                        <  2:2603756/76‑1 (MQ=255)
                                 aTGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaaca                                       <  1:3166573/76‑1 (MQ=255)
                                 aTGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaaca                                       <  1:3002413/76‑1 (MQ=255)
                                 aTGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaaca                                       >  1:2226288/1‑76 (MQ=255)
                                  tGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGAt                                                           <  2:5302266/55‑1 (MQ=255)
                                    gtgtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCCACAAAt                                    >  2:3685696/1‑76 (MQ=255)
                                    gtgtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGGCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAAt                                    >  2:4093684/1‑76 (MQ=255)
                                     tgtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATc                                   <  1:4181399/76‑1 (MQ=255)
                                      gtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCt                                  <  1:3838725/76‑1 (MQ=255)
                                      gtTCTGGGATCGCTTAAACCGAGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCGACATGCGCAACAAATCt                                  >  2:1072791/1‑76 (MQ=255)
                                      ggtCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCGACATGGGCAACAAATCt                                  >  2:4888162/3‑76 (MQ=255)
                                       ttCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTg                                 <  2:1474577/76‑1 (MQ=255)
                                       ttCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTg                                 <  2:4651902/76‑1 (MQ=255)
                                        tCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTgg                                <  1:527255/76‑1 (MQ=255)
                                         cgggggTCGCTGAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGa                               <  2:63716/74‑1 (MQ=255)
                                         cTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGa                               >  1:582298/1‑76 (MQ=255)
                                          tGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGAt                              >  2:2890778/1‑76 (MQ=255)
                                           gggATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGAtt                             >  2:2876660/1‑76 (MQ=255)
                                            ggATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTc                            <  1:2442766/76‑1 (MQ=255)
                                            ggATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTc                            <  2:4476564/76‑1 (MQ=255)
                                             gATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCa                           <  1:1229457/76‑1 (MQ=255)
                                             gATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCa                           >  1:5450494/1‑76 (MQ=255)
                                               tCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAgg                         <  2:4773428/76‑1 (MQ=255)
                                                cGCTTGAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAggg                        <  1:700035/76‑1 (MQ=255)
                                                cGCTTAAACCGGGCAAACGCGCCGGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAggg                        <  2:2226288/76‑1 (MQ=255)
                                                 gCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGt                       <  1:2984527/76‑1 (MQ=255)
                                                  catAACCCGGGGAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTc                      <  1:2867463/74‑1 (MQ=255)
                                                  cTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGTATTCAGGGTc                      <  1:1145487/76‑1 (MQ=255)
                                                  cTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTc                      <  2:2487981/76‑1 (MQ=255)
                                                  cTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTc                      <  2:5575086/76‑1 (MQ=255)
                                                   ttAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCa                     <  1:457604/76‑1 (MQ=255)
                                                   ttAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCa                     <  2:4923652/76‑1 (MQ=255)
                                                   ttAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCa                     >  2:3804061/1‑76 (MQ=255)
                                                       aCCGGGGCAACTCGCGAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCCGGGTCAATTa                 <  1:1281257/76‑1 (MQ=255)
                                                       aCCGGGCACACGCGCCAGCGTCGTTGCGCTGGCTCGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAAttc                 <  2:3808900/76‑2 (MQ=255)
                                                       aCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTa                 >  1:4087760/1‑76 (MQ=255)
                                                       aCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTa                 <  1:1871400/76‑1 (MQ=255)
                                                          gggCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCt              >  1:3205976/1‑76 (MQ=255)
                                                          gggCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCt              >  1:2266445/1‑76 (MQ=255)
                                                           ggCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCCGGGTCAATTAGCtt             >  2:2047010/1‑76 (MQ=255)
                                                            gCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTc            >  1:1081060/1‑76 (MQ=255)
                                                            gCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTc            >  1:5698540/1‑76 (MQ=255)
                                                            gCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTc            >  2:5229802/1‑76 (MQ=255)
                                                             cAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCg           >  1:2180747/1‑76 (MQ=255)
                                                             cAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCg           <  2:4957286/76‑1 (MQ=255)
                                                             cAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCg           <  1:1235415/76‑1 (MQ=255)
                                                              aaaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGt          >  2:215149/1‑76 (MQ=255)
                                                              aaaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGt          <  2:3450277/76‑1 (MQ=255)
                                                               aaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGtt         >  2:4286987/1‑76 (MQ=255)
                                                                aCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGttt        <  1:4666356/76‑1 (MQ=255)
                                                                 cgcgCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGtttt       >  2:461603/1‑76 (MQ=255)
                                                                    gcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGAt    >  1:4534236/1‑76 (MQ=255)
                                                                     ccAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATa   <  2:3274781/76‑1 (MQ=255)
                                                                      cAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAg  <  2:2804083/76‑1 (MQ=255)
                                                                         |                                                                        
CGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGAGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAG  >  NC_000913/3279124‑3279269

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