Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,905,761 G→A G25D (GGT→GAT)  mntP → Mn(2(+)) exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,905,7610GA100.0% 138.1 / NA 37G25D (GGT→GAT) mntPMn(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/19);  total (18/19)

GGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATGAAGATGAAGA  >  NC_000913/1905530‑1905980
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                           
ggCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCg                                                                                                                                                                                                          >  2:598740/1‑251 (MQ=255)
    cGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAAt                                                                                                                                                                                                      >  1:394391/1‑251 (MQ=255)
       aCGGACTGAATACTTTCTGTCAGGCGAGACCATATGCACATAAATCGCTATGCCTGCATGAACGCAAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGGTGGGGGCCATCATGTGTTTTATGGATACCCGGGTCCGGAAATTGTCATGAAGATCAATGCTAATGTTCTTCTTGCGTTTGGTATGTCGATAGATGCAATTGCTCTATCAATCGGAAATGATGCTCCCCACCATAAACCGAAAtttt                                                                                                                                                                                                   >  2:436280/1‑251 (MQ=255)
       aCGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAAtttt                                                                                                                                                                                                   >  1:442476/1‑251 (MQ=255)
         ggACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCt                                                                                                                                                                                                 >  1:115789/1‑251 (MQ=255)
              gAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                                                                                                                                                                                            <  2:349592/251‑1 (MQ=255)
                                  gaCCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAttt                                                                                                                                                                        <  1:2737/251‑1 (MQ=255)
                                    ccATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAttttt                                                                                                                                                                      <  2:101565/251‑1 (MQ=255)
                                                     gCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCt                                                                                                                                                     >  1:115877/1‑251 (MQ=255)
                                                        aTGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGAc                                                                                                                                                  >  1:552629/1‑251 (MQ=255)
                                                           ccTGCATGAACGCGAAAACCATCACGGCTTGTTGGTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGAcgcc                                                                                                                                               <  2:105578/251‑1 (MQ=255)
                                                                      cctaaaaCCATCACTGTTTGTTGTTCATGTATGCATTTTTGCTGGGGGCGATGATGTGTTTTATGGATCCCCCGGTCAGTACATTGTCATGAATATCACTGCTACTGTTGTACTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGTGAACCGGCCTAATTTTTGGTACCGTCGAAACCCGGACGCCGCAGATCGGCt                                                                                                                                    <  2:300109/248‑1 (MQ=255)
                                                                           aaCCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTggggg                                                                                                                               >  2:352723/1‑250 (MQ=255)
                                                                                 cgcgtgttgttgTTCAGTTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCAGGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAATCCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATggg                                                                                                                          <  2:549063/244‑1 (MQ=255)
                                                                                         gttgttCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCTGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTagc                                                                                                                 <  1:525831/251‑1 (MQ=255)
                                                                                              tCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTagccagc                                                                                                             <  2:339555/250‑1 (MQ=255)
                                                                                                 tGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGt                                                                                                         >  1:329174/1‑251 (MQ=255)
                                                                                                                     gCGATGATTTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCTCTGCTCCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCAc                                                                                     <  2:303747/251‑1 (MQ=255)
                                                                                                                            tgtgTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTg                                                                              >  1:46570/1‑251 (MQ=255)
                                                                                                                                    aTGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTtgtg                                                                      <  2:180547/251‑1 (MQ=255)
                                                                                                                                      ggATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGtgct                                                                    >  2:21845/1‑251 (MQ=255)
                                                                                                                                             ccGGTCAGGACCTTGTCATGAATATCACTGCTACTGTTCTTCTTGCGATTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGCTCGGCTGGGGCATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGGTTGCGTTTGTGCTGCTGat                                                              >  2:421601/1‑250 (MQ=255)
                                                                                                                                               ggTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATAtt                                                           >  2:31158/1‑251 (MQ=255)
                                                                                                                                                gTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTc                                                          >  1:255816/1‑251 (MQ=255)
                                                                                                                                                   aGGACATTGTCATGAATATCACGGCTTCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGGGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGATGATATTCCTc                                                       <  2:46570/251‑1 (MQ=255)
                                                                                                                                                    ggACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTcg                                                      >  2:554919/1‑251 (MQ=255)
                                                                                                                                                     gACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTcgg                                                     <  1:547202/251‑1 (MQ=255)
                                                                                                                                                       cATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTcggcg                                                   >  1:440493/1‑251 (MQ=255)
                                                                                                                                                              aTGAATATCACTGCTACTGTTCTTCTGGCGTTTGGTAGGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAAt                                            <  2:43114/251‑1 (MQ=255)
                                                                                                                                                                gAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATg                                           <  1:598740/250‑1 (MQ=255)
                                                                                                                                                                gAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATg                                           <  1:519718/250‑1 (MQ=255)
                                                                                                                                                                      cACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTg                                    <  1:554919/251‑1 (MQ=255)
                                                                                                                                                                                     tcttGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGgc                      <  1:404623/250‑1 (MQ=255)
                                                                                                                                                                                           gTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTg                                                                                                                                           >  2:169746/1‑127 (MQ=255)
                                                                                                                                                                                           gTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTg                                                                                                                                           <  1:169746/127‑1 (MQ=255)
                                                                                                                                                                                                        atggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGGACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGGTTTCGGGGGGCCGGTGatgaagatgaagg  >  2:97534/1‑250 (MQ=255)
                                                                                                                                                                                                        atggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGAtgaagatgaaga  <  1:352723/251‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                           
GGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATGAAGATGAAGA  >  NC_000913/1905530‑1905980

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