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.3% 353.1 / NA 136R145S (AGA→AGCagaAputative N‑acetylgalactosamine‑6‑phosphate deacetylase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (57/77);  total (58/78)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

GTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGAGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATA  >  NC_000913/3279125‑3279268
                                                                        |                                                                       
 tCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCg                                                                     <  1:4814918/76‑1 (MQ=255)
  cGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTAGGATCGCTTAAACCGGGCAAACGCGCCAGcgtcgt                                                                    <  1:3755498/76‑1 (MQ=255)
   gCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGtt                                                                   >  1:2540/1‑76 (MQ=255)
   gCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGtt                                                                   >  1:3305460/1‑76 (MQ=255)
    cTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTg                                                                  <  1:2604821/76‑1 (MQ=255)
    cTGCAGCCGGCGCGAATGCTGTGTGTTGATGGTGTTCTGGGATCGCTTAAACAGGGCAAACGCGCCAGCGTCGTTg                                                                  <  1:2609335/76‑1 (MQ=255)
     tGCATCCGGCGCGAATGCTGGTTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTgc                                                                 <  1:1800969/76‑1 (MQ=255)
        aTCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCt                                                              <  1:1813749/76‑1 (MQ=255)
        aTCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCt                                                              <  2:1891114/76‑1 (MQ=255)
        aTCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCt                                                              >  2:1521195/1‑76 (MQ=255)
        aTCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCt                                                              <  1:4156425/76‑1 (MQ=255)
        aTCCGGCGCGAATGCTGGGTGTTGATGATGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCt                                                              >  2:2175281/1‑76 (MQ=255)
         tCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTg                                                             <  2:406074/76‑1 (MQ=255)
         tCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTg                                                             >  1:4060027/1‑76 (MQ=255)
          ccGGCGCGAATGCTGGGTGTTGACGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTgg                                                            >  1:922380/1‑76 (MQ=255)
           cGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGa                                                           <  1:4802330/76‑1 (MQ=255)
           cGGCGCGAATGCTGGGCGTTGATGGTGTTCTGGGATCGCGTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGa                                                           >  1:4471962/1‑76 (MQ=255)
             gcgcGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATa                                                         <  2:4696279/76‑1 (MQ=255)
               gcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAgg                                                       >  1:298707/1‑75 (MQ=255)
               gcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGc                                                       <  1:618533/76‑1 (MQ=255)
               gcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGc                                                       >  2:3051057/1‑76 (MQ=255)
               gcgAAGGCGGGGAGTTGATGGAGTTCGGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTGGAGCTGGATaac                                                       >  2:3631847/1‑74 (MQ=255)
                 gAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCgg                                                     <  1:2901759/76‑1 (MQ=255)
                 gAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCgg                                                     <  1:2168611/76‑1 (MQ=255)
                 gAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCgg                                                     <  1:1692866/76‑1 (MQ=255)
                  aaTGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCggg                                                    <  1:88952/76‑1 (MQ=255)
                  aaTGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCggg                                                    <  1:385505/76‑1 (MQ=255)
                   aTGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGc                                                   <  1:4831579/76‑1 (MQ=255)
                    tGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCt                                                  <  1:1240600/76‑1 (MQ=255)
                     gCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCggg                                                    <  2:3641495/73‑1 (MQ=255)
                     gCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTa                                                 >  1:345859/1‑76 (MQ=255)
                     gCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTa                                                 <  2:4529809/76‑1 (MQ=255)
                     gCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTa                                                 >  1:4739/1‑76 (MQ=255)
                      cTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTAc                                                <  1:2921697/76‑1 (MQ=255)
                       tGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCcaca                                               >  2:649851/1‑76 (MQ=255)
                        gggTGTTGATGGTGTTCTGGGATCGCTCAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACAt                                              <  2:1646244/76‑1 (MQ=255)
                         ggTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACAtg                                             >  1:3003404/1‑76 (MQ=255)
                         ggTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACAtg                                             <  1:729146/76‑1 (MQ=255)
                          gtgtTGCTGGCGTTCTGGGATCGCTTAAACCGGGCAAACGCGCAAGAGGCGTGGCGCGGGTTAGCGGGTTGCAtgt                                            >  1:4738228/1‑76 (MQ=255)
                          gtgtTAATGATGTTCTGGAATCGCTTCGCCCGGTGAAGCGGTCGAGCGTCGTTGCGCTGGATAGCGGGCTACAtgt                                            <  1:2102190/76‑1 (MQ=255)
                           ggtTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACAtgtg                                           <  2:1913503/75‑1 (MQ=255)
                            gtTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGc                                          >  2:3984158/1‑76 (MQ=255)
                             ttGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATGGGGGGCCACTTGTGca                                         >  2:4838946/1‑76 (MQ=255)
                             ttGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGca                                         <  1:2171133/76‑1 (MQ=255)
                             gtGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGca                                         <  1:469795/75‑1 (MQ=255)
                               gATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaac                                       >  1:1702613/1‑76 (MQ=255)
                               gATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaac                                       <  2:561842/76‑1 (MQ=255)
                                gtGGTGTTCTGGGGTCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaaca                                      <  2:1476898/75‑1 (MQ=255)
                                aTGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAg                                                        <  1:4615106/58‑1 (MQ=255)
                                 gggTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaacaa                                     <  2:4042796/75‑1 (MQ=255)
                                  ggTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACaaa                                    >  1:4732284/1‑76 (MQ=255)
                                   gtgtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAAt                                   <  1:163841/76‑1 (MQ=255)
                                   gtgtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAAt                                   <  1:3161566/76‑1 (MQ=255)
                                     gtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCt                                 <  1:4651526/76‑1 (MQ=255)
                                     gtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCt                                 >  2:3808621/1‑76 (MQ=255)
                                      ttCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTg                                <  1:2945014/76‑1 (MQ=255)
                                       tCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATCGCGGGCTACATGTGCAACAAATCTgg                               >  2:4329370/1‑76 (MQ=255)
                                       tCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTgg                               >  1:2890488/1‑76 (MQ=255)
                                       tCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTgg                               >  2:4324823/1‑76 (MQ=255)
                                        cTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGa                              >  2:4577257/1‑76 (MQ=255)
                                         tGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGAt                             <  1:3708584/76‑1 (MQ=255)
                                          gggATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGAtt                            <  1:4682478/76‑1 (MQ=255)
                                          gggATCGCTGAACCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGAtt                            <  2:148686/76‑1 (MQ=255)
                                            gATCGCTTAAACCGGGGAGACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCa                          <  1:2366533/76‑1 (MQ=255)
                                            gATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTcg                          >  2:1840613/1‑75 (MQ=255)
                                            gATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCa                          >  1:943695/1‑76 (MQ=255)
                                            gATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCa                          <  2:881403/76‑1 (MQ=255)
                                            gATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCa                          <  1:775848/76‑1 (MQ=255)
                                            gATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCa                          >  2:2879016/1‑76 (MQ=255)
                                               cGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAggg                       >  1:4323990/1‑76 (MQ=255)
                                                gCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGt                      <  1:2870376/76‑1 (MQ=255)
                                                gCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGt                      <  1:4146055/76‑1 (MQ=255)
                                                  ttAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCa                    >  2:3017932/1‑76 (MQ=255)
                                                  ttAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCa                    <  1:670321/76‑1 (MQ=255)
                                                  ttAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCa                    <  2:2146969/76‑1 (MQ=255)
                                                   tAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCaa                   >  2:1994824/1‑76 (MQ=255)
                                                    aaaCCGGGCAAACGCGCCAGCTGCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAAt                  <  2:3539728/76‑1 (MQ=255)
                                                    aaaCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAGTCTGGATTCAGGGTCAAt                  >  1:1393103/1‑76 (MQ=255)
                                                    aaaCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAAt                  <  2:3496543/76‑1 (MQ=255)
                                                    aaaCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAAt                  >  2:1620250/1‑76 (MQ=255)
                                                      cgcGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTa                <  2:2762580/74‑1 (MQ=255)
                                                       ccGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAg               <  1:1069047/76‑1 (MQ=255)
                                                       ccGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAg               <  1:2749891/76‑1 (MQ=255)
                                                         gggCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCt             <  1:2629859/76‑1 (MQ=255)
                                                         gggAAAGCGCGCGAGGGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCt             <  2:5184/76‑1 (MQ=255)
                                                          ggCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTTGCtt            >  2:904248/1‑76 (MQ=255)
                                                           gCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTc           <  2:3905167/76‑1 (MQ=255)
                                                            cAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCg          >  1:1738151/1‑76 (MQ=255)
                                                              aaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGtt        >  2:1891261/1‑76 (MQ=255)
                                                               aCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGttt       >  2:3112095/1‑76 (MQ=255)
                                                               aCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGttt       <  1:2101681/76‑1 (MQ=255)
                                                                 gcgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTg     >  1:3451732/1‑76 (MQ=255)
                                                                  cgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTgg    >  2:4451310/1‑75 (MQ=255)
                                                                  cgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGa    <  2:3118772/76‑1 (MQ=255)
                                                                   gcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGAt   <  1:4590115/76‑1 (MQ=255)
                                                                   gcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGAt   >  1:392712/1‑76 (MQ=255)
                                                                   gcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGAt   >  1:3885883/1‑76 (MQ=255)
                                                                    ccAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATa  >  1:4318303/1‑76 (MQ=255)
                                                                    ccAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATa  <  2:4666305/76‑1 (MQ=255)
                                                                    ccAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATa  >  2:1104091/1‑76 (MQ=255)
                                                                        |                                                                       
GTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGAGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATA  >  NC_000913/3279125‑3279268

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