Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 124,938 A→G I15V (ATT→GTT)  pdhR → transcriptional regulator of pyruvate dehydrogenase complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606124,9380AG97.3% 777.6 / 6.6 188I15V (ATT→GTT) pdhRtranscriptional regulator of pyruvate dehydrogenase complex
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (104/79);  total (106/82)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

Only 100 of 190 total aligned reads displayed.

GACTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAAC  >  REL606/124865‑125013
                                                                         |                                                                           
gACTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGtt                                                                           <  2:2780521/76‑1 (MQ=255)
 aCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTg                                                                          <  1:4870839/76‑1 (MQ=255)
 aCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTg                                                                          <  1:1740572/76‑1 (MQ=255)
 aCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTg                                                                          <  1:6048035/76‑1 (MQ=255)
  cTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGa                                                                         >  1:2934733/1‑76 (MQ=255)
   ttCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGag                                                                        >  2:2363235/1‑76 (MQ=255)
   ttCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGag                                                                        <  2:4813221/76‑1 (MQ=255)
   ttCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGag                                                                        >  1:2024518/1‑76 (MQ=255)
   ttCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGag                                                                        <  2:4982009/76‑1 (MQ=255)
   ttCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGag                                                                        >  2:5885/1‑76 (MQ=255)
    tCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGagc                                                                       <  1:398036/76‑1 (MQ=255)
       gCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGagcagc                                                                    >  1:1382326/1‑76 (MQ=255)
       gCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGagcagc                                                                    >  2:5487632/1‑76 (MQ=255)
        cAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAgcagca                                                                   >  2:1867645/1‑76 (MQ=255)
        cAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAgcagca                                                                   <  1:2129806/76‑1 (MQ=255)
         aaGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCaa                                                                  >  1:5306337/1‑76 (MQ=255)
          aGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAAc                                                                 >  1:3451675/1‑76 (MQ=255)
           gTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACt                                                                >  1:692060/1‑76 (MQ=255)
           gTGGCTAAAGACAGGAACTCAGGCCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACt                                                                <  2:2207620/76‑1 (MQ=255)
             ggCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTgg                                                              <  2:6761496/76‑1 (MQ=255)
              gCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGa                                                             >  2:5671484/1‑76 (MQ=255)
                ttAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGt                                                           >  2:1633654/1‑76 (MQ=255)
                ttAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGt                                                           >  1:6764898/1‑76 (MQ=255)
                ttAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGt                                                           >  1:4464286/1‑76 (MQ=255)
                ttAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGt                                                           >  1:366167/1‑76 (MQ=255)
                 tAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGtt                                                          >  2:292515/1‑76 (MQ=255)
                   aGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGtttt                                                        >  1:5338287/1‑76 (MQ=255)
                    gACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGCGCAGCAACTGGAGttttt                                                       >  1:2843175/1‑76 (MQ=255)
                     aCAGGAACTCAGGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTg                                                      <  1:2717978/76‑1 (MQ=255)
                        ggACCTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATc                                                   <  2:4596852/76‑1 (MQ=255)
                         gAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATcc                                                  <  2:6687626/76‑1 (MQ=255)
                          aaCTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTGATCCt                                                 <  2:5995077/76‑1 (MQ=255)
                           aCTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTc                                                >  2:6878461/1‑76 (MQ=255)
                           aCTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTc                                                <  2:6362957/76‑1 (MQ=255)
                            cTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCg                                               <  1:4183438/76‑1 (MQ=255)
                            cTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCg                                               <  2:4753376/76‑1 (MQ=255)
                             tCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGa                                              >  1:3758439/1‑76 (MQ=255)
                             tCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGa                                              >  2:5641827/1‑76 (MQ=255)
                             tCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGa                                              >  1:1498513/1‑76 (MQ=255)
                              cATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGaa                                             >  1:4734365/1‑76 (MQ=255)
                              cATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGaa                                             <  2:1397699/76‑1 (MQ=255)
                               aTGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAg                                            >  1:2519583/1‑76 (MQ=255)
                                tGGCCTACAGAAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAgg                                           <  2:2524568/76‑1 (MQ=255)
                                 ggCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGc                                          >  2:3180238/1‑76 (MQ=255)
                                 ggCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGc                                          >  1:1292227/1‑76 (MQ=255)
                                  gCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCa                                         >  2:6623011/1‑76 (MQ=255)
                                   ccTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCAc                                        >  2:1066153/1‑76 (MQ=255)
                                   ccTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCAc                                        >  2:4251256/1‑76 (MQ=255)
                                   ccTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCAc                                        >  1:2904100/1‑76 (MQ=255)
                                    cTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCAct                                       >  1:5514774/1‑76 (MQ=255)
                                    cTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCAct                                       <  1:2778562/76‑1 (MQ=255)
                                     tACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCActc                                      >  1:4717871/1‑76 (MQ=255)
                                      aCAGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCActct                                     >  2:3797248/1‑76 (MQ=255)
                                      aCAGCAAAACCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCActct                                     <  2:4903031/76‑1 (MQ=255)
                                        aGCAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTcc                                   >  2:1692634/1‑76 (MQ=255)
                                         gCAAAACCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCg                                  <  1:6576406/76‑1 (MQ=255)
                                          cAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGc                                 <  1:4370989/76‑1 (MQ=255)
                                          cAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGc                                 <  2:3129386/76‑1 (MQ=255)
                                          cAAAATCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGc                                 <  2:2689396/76‑1 (MQ=255)
                                           aaaaTCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGcc                                >  2:2601895/1‑76 (MQ=255)
                                             aacccGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCg                              <  2:22764/73‑1 (MQ=255)
                                             aaTCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCg                              >  1:232940/1‑76 (MQ=255)
                                              aTCCGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCgg                             >  1:1219782/1‑76 (MQ=255)
                                              aTCCCCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCgg                             <  1:649194/76‑1 (MQ=255)
                                                 cGCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCg                          <  1:948597/76‑1 (MQ=255)
                                                  gCCAACCAAAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGa                         >  2:3298360/1‑76 (MQ=255)
                                                   ccaaccaaACCTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGaa                        <  2:3577728/76‑1 (MQ=255)
                                                   ccaaccaaAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGaa                        >  2:1437232/1‑76 (MQ=255)
                                                    caaccaaAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGaaa                       <  2:418236/76‑1 (MQ=255)
                                                     aaccaaAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGaaaa                      <  2:1589243/76‑1 (MQ=255)
                                                      accaaAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGaaaaa                     <  2:2343885/76‑1 (MQ=255)
                                                       ccaaACCTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCCCTCTCCGCCCGGGCGAAAAAc                    <  1:5560395/76‑1 (MQ=255)
                                                        caaAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACt                   >  2:6440113/1‑76 (MQ=255)
                                                        caaAACTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACt                   >  1:3651871/1‑76 (MQ=255)
                                                         aaaaCTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTc                  >  1:5919824/1‑76 (MQ=255)
                                                          aaaCTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTcc                 >  1:1807545/1‑76 (MQ=255)
                                                            aCTCTCCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCa               >  2:1985725/1‑76 (MQ=255)
                                                              tctcCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCAcc             >  2:4039765/1‑76 (MQ=255)
                                                              tctcCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCAcc             <  2:1768058/76‑1 (MQ=255)
                                                               ctcCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCg            >  1:3935833/1‑76 (MQ=255)
                                                               ctcCGATGTGGTTGACCAGCAACTGGAGTTTTTGATCCTCGAAGACACTCTCCGCCCGGGCGAAAAACTCCCACCg            <  2:1211003/76‑1 (MQ=255)
                                                                tcCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCgg           <  2:3656459/76‑1 (MQ=255)
                                                                tcCGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCgg           >  2:1203805/1‑76 (MQ=255)
                                                                tcCGATGGGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCgg           <  1:3578003/76‑1 (MQ=255)
                                                                 ccGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGa          <  2:5846397/76‑1 (MQ=255)
                                                                 ccGATGTGATTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGa          >  1:164688/1‑76 (MQ=255)
                                                                  cGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGaa         >  1:229959/1‑76 (MQ=255)
                                                                  cGATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGaa         >  2:45963/1‑76 (MQ=255)
                                                                   gATGTGGTTGAGCAGCAACTTGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAAc        <  2:3437969/76‑1 (MQ=255)
                                                                   gATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAAc        >  2:5191273/1‑76 (MQ=255)
                                                                   gATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAAc        <  1:5835143/76‑1 (MQ=255)
                                                                   gATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAAc        >  1:816581/1‑76 (MQ=255)
                                                                   gATGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAAc        <  1:1959772/76‑1 (MQ=255)
                                                                    aTGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAGCTCCCACCGGAAcg       <  2:932739/76‑1 (MQ=255)
                                                                    aTGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAAcg       <  2:1838377/76‑1 (MQ=255)
                                                                    aTGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAAcg       >  2:4002143/1‑76 (MQ=255)
                                                                    aTGTGGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAAcg       >  1:4567653/1‑76 (MQ=255)
                                                                      gtgGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAAcgcg     >  1:3867411/1‑76 (MQ=255)
                                                                      gtgGTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAAcgcg     >  2:4456804/1‑76 (MQ=255)
                                                                         aTTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAAc  <  1:5008715/76‑1 (MQ=255)
                                                                         |                                                                           
GACTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAAC  >  REL606/124865‑125013

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