Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,194 A→T 8.8% intergenic (+28/+25) yicH → / ← yicI AsmA famly protein YicH/alpha‑D‑xyloside xylohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,1940AT 8.8% 146.8 / 6.9 60intergenic (+28/+25)yicH/yicIAsmA famly protein YicH/alpha‑D‑xyloside xylohydrolase
Reads supporting (aligned to +/- strand):  ref base A (27/25);  new base T (3/2);  total (32/28)
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.68e-01

ACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTC  >  NC_000913/3832053‑3832336
                                                                                                                                             |                                                                                                                                              
aCGGTTAGCGGGCAATTGCATGCCGTGGACGCTGCAAAGCAGCAAGTGACTTAAAACTAAAAACCAATCATCCCGAGCATCGGGGTAGTGACATCGACGGAGCCGGTTCGGTAAACATCAATCTCGCCGGGTCCGACTCacac                                                                                                                                               >  2:188509/1‑143 (MQ=255)
           gCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAGGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGGCTCAGACGTACCCgggct                                                                                                                                    >  2:365017/1‑140 (MQ=255)
             acttGCATGCCTTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcgaga                                                                                                                                  <  2:503609/141‑5 (MQ=255)
             aaTTGCATGCCGTGAACGCCGCGAAGCAGGAAGGGACTCAAACCATGAATGCGGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcgaga                                                                                                                                  <  2:406939/143‑5 (MQ=255)
               ttGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGCCTCCGACGTACCCGGCATCCCa                                                                                                                                >  1:444684/1‑143 (MQ=255)
                tGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                                                                               <  1:788519/143‑1 (MQ=255)
                tGCAGGGCGGGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                                                                               <  1:406203/143‑1 (MQ=255)
                  cATGCCGTGCACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCa                                                                                                                             <  1:522190/143‑1 (MQ=255)
                     gCCGCGAACGCGGCACTGCAGTATGTGACTCAAACCATCAATGCTGGCATCGTTTCCAGCGGTGAAGTTACAGCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATa                                                                                                                          <  1:159317/143‑1 (MQ=255)
                            aCGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACa                                                                                                                   >  1:463062/1‑143 (MQ=255)
                                cgcAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTCCCCGGCATTCCATCAATCAGTTACAACGt                                                                                                               >  2:221803/1‑143 (MQ=255)
                                cgcAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCAGGCAGTCCATCAATAGATTACAACGt                                                                                                               >  2:221843/1‑143 (MQ=255)
                                           aaGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACCTCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCCTCAATAGATTACAACGTAATTGTCAgcg                                                                                                    >  2:367692/1‑143 (MQ=255)
                                           aaGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcg                                                                                                    >  2:681343/1‑143 (MQ=255)
                                           aaGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcg                                                                                                    <  1:414179/143‑1 (MQ=255)
                                              tGACTCACACCATGAATGCTGGCATCGTTTCAAGCGCTGAAGGTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                                                                                 <  2:865595/143‑1 (MQ=255)
                                              tGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTACACATCAATCTCGCCCGGTACGACTCAGACGTACCCGGCATTCCATCAATCGATTACAACCTACTTGTCACCGCAt                                                                                                 >  2:427296/1‑143 (MQ=255)
                                                 cTCAAACCATGAAGGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTcc                                                                                              <  2:271355/143‑1 (MQ=255)
                                                  tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCATTCGCGCCGGGTACGACGCAGACGTAGCGGCCATTCCGTCAATAGGTTCCAAGGGAATTGTCAGCGCATTccc                                                                                             >  2:250464/1‑143 (MQ=255)
                                                  tCAAACCATGAATGCTGGCATCGGTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACAGCAATCTCGCCGGGTCGCACTCAAGCGTACGCGACATTCAATCAATAGAATGCAACGACCAAGTCAGCGCAGTccc                                                                                             >  2:300621/1‑143 (MQ=255)
                                                    aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                                                                           <  1:24396/143‑1 (MQ=255)
                                                    aaaCCATCACTGCTGGCCTCGCTTCCAGCGGTGAAGTTACCACCACGTAGCTTTTGCGGTCAACAACCATCTCGCCTGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                                                                           <  1:608908/143‑1 (MQ=255)
                                                     aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                                                                          >  1:413287/1‑143 (MQ=255)
                                                      aCCATGAATGCTGGCATCGTTTCCAGCTGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                                                                         <  1:431299/143‑1 (MQ=255)
                                                           gAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAACCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCtt                                                                                    <  1:449251/143‑1 (MQ=255)
                                                             aTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                                                                  <  1:428544/143‑1 (MQ=255)
                                                              tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcac                                                                                 <  1:811760/143‑1 (MQ=255)
                                                              tGCTGGCATCGTTTCCAGCGGTGAAGCTACATCGACGGCGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcac                                                                                 <  1:86887/143‑1 (MQ=255)
                                                                cTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGAc                                                                                                                                               >  2:228284/1‑79 (MQ=255)
                                                                cTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGAc                                                                                                                                               <  1:228284/79‑1 (MQ=255)
                                                                   gCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATCCCCTTGAGGCTTCAccacca                                                                            >  1:453943/1‑143 (MQ=255)
                                                                              ggCGGTGAAGTTACATCGACGGGGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCa                                                                 <  2:752668/142‑1 (MQ=255)
                                                                              aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATACATCGATTCCAACGTAATTGTCAGCGCATTCCCTCGCGGCTTCACCACCACCCCCTGcccc                                                                 >  2:279434/1‑140 (MQ=255)
                                                                                      agtttgaTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGc                                                         <  2:585163/137‑1 (MQ=255)
                                                                                         tACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTg                                                      <  1:607743/143‑1 (MQ=255)
                                                                                          acttCGACGGAGCCGGTGCGGGAAACATCATGCTCGCCGGGTACGACTCAGGCGTACCCGGCATTCCATCAATCGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc                                                     <  2:682169/140‑1 (MQ=255)
                                                                                            aTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGa                                                   >  1:514747/1‑143 (MQ=255)
                                                                                             tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCCGACCCTGCGCACTTTCAGCCTGCGaa                                                  >  2:279258/1‑143 (MQ=255)
                                                                                                      ccGGTGCGGTAAACATCGATGTCTCCGGGTACGTCTCATACGTGCCGTGCAATCCATCCATAGAGTTCAACTGATTTTGCAGCGCATTCCCTTGGTGCTGCACCACCAGCCCCTGCTCACTTTCATCCTGGGCACCGTCTTGc                                         >  2:854687/1‑143 (MQ=255)
                                                                                                         gcgCGGTAAACATCCATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGTTTCACTACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGa                                      <  1:37579/141‑1 (MQ=255)
                                                                                                            cGGTAAACATCAATCTCGCCGGGTACGACTAAGCCGTACCCGCCCTGCCAACAATAGATTACAACGTAATTGTCAGCGCATGCCCTTGAGGCTTCACCACCCGCGCCTGCTCACTTTCAGCCAGCGCCCCGTCTCGCCGACCa                                   >  1:378266/1‑143 (MQ=255)
                                                                                                              gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCACACTTGCGGCCTGCGAACCGTCTTGCGGACCAtt                                 >  1:691859/1‑143 (MQ=255)
                                                                                                                     tCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCACTAGATTACAACGTCATTCTCCGCGCACTCCCTTCATGCTTCACCACCACCCCCTGCTCACTTTCAGCCTGCGAACCCTCTTGCAGACCACTCACttgc                          >  2:173146/1‑141 (MQ=255)
                                                                                                                       aaTCTCGCCGGGTCCGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTcac                        <  1:188509/143‑1 (MQ=255)
                                                                                                                       aaTCTCCTCGGGTACTACCCACACGTACCCGCCATTCCTTCATTAGATTGCATCGCAATTGCAAGCTTTATCTCTTGAGGCGTCACCACCAGCCCCTGTTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTcac                        <  2:58439/143‑1 (MQ=255)
                                                                                                                           tcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCacaaca                    >  2:764094/1‑143 (MQ=255)
                                                                                                                           tcGCCGGGGACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCCCCAGCACCTGCTCACTCCCAGCCTGCGCACCGTCTTGCAAACCATTTACTTTCccacca                    >  2:645718/1‑143 (MQ=255)
                                                                                                                             gCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTGTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACAtt                  <  1:783697/143‑1 (MQ=255)
                                                                                                                               cGGGTACGTCTGAGTCGTACCCGGCCTTCCATCAATAGATTACAACGTACTTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACAttcc                >  2:206253/1‑141 (MQ=255)
                                                                                                                               cGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTTCAGACCATTTAATTTCACAACATTgc                >  1:423124/1‑143 (MQ=255)
                                                                                                                               cGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACCTGCAGCCTGCGAAGCGACCTCCACACCATTTACTTTCACACCATTgc                >  2:669885/1‑143 (MQ=255)
                                                                                                                                gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTACTTGTCACCGCATTCCCTCGAGGCTTCACCACCCGCCCCTGCCCACCTTCAGCCTGCGCACCGCCTAGCCGACCCCTTACTTTCACAGCAtcccc               >  2:533984/1‑139 (MQ=255)
                                                                                                                                gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCGCTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTGTCCCAGCATTgcg               >  1:392759/1‑143 (MQ=255)
                                                                                                                                 tgTACGACGCAGACGTCCCTGACATTCCATGAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGATGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTgcgc              <  2:817519/142‑1 (MQ=255)
                                                                                                                                   tACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAg            >  1:784852/1‑143 (MQ=255)
                                                                                                                                     cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGCACCGTCTTGCCGACCACTTACTTTCACTACTTTGCGCATGc          >  1:273164/1‑143 (MQ=255)
                                                                                                                                        cTCAGACGTACCCGGCATTCCACCAACAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGcaca       <  1:487838/143‑1 (MQ=255)
                                                                                                                                         tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCtt                                                                                    >  2:889195/1‑65 (MQ=255)
                                                                                                                                         tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCtt                                                                                    <  1:889195/65‑1 (MQ=255)
                                                                                                                                         tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAg      >  1:400514/1‑143 (MQ=255)
                                                                                                                                         tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTACCCCGAAGCGTCACCACCAGCACCTGCTCACTCTCAGCCTGCGAACCGTCTTGCAGACCATCTACTTTCACAACATTGCGCAGGcacaa      >  2:397163/1‑142 (MQ=255)
                                                                                                                                             aCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTc  >  1:55472/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                              
ACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTC  >  NC_000913/3832053‑3832336

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