Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,282,577 C→T 24.9% P905L (CCG→CTG)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,282,5770CT24.9% 137.8 / 43.8 82P905L (CCG→CTG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base T (10/10);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGG  >  NC_000913/1282480‑1282672
                                                                                                 |                                                                                               
ggCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGccctc                                                                                               >  2:1357619/1‑98 (MQ=255)
ggCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGc                                                                                               >  2:1971266/1‑100 (MQ=255)
 gCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCt                                                                                              <  2:1108399/100‑1 (MQ=255)
 gCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCt                                                                                              <  2:1805790/100‑1 (MQ=255)
     gTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg                                                                                          <  2:2128741/100‑1 (MQ=255)
      tGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCTGCTGATgg                                                                                         <  1:1752267/100‑1 (MQ=255)
       gCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATgg                                                                                         <  2:868141/99‑1 (MQ=255)
        cGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGgag                                                                                       <  1:1641472/100‑1 (MQ=255)
        cGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGACGACTTACGAACGCTTTACCTCTATCGGCCAGCTGATGgag                                                                                       >  1:1693952/1‑100 (MQ=255)
         gACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGgaga                                                                                      >  1:980237/1‑100 (MQ=255)
         gACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGgaga                                                                                      >  1:1576706/1‑100 (MQ=255)
         gACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGgaga                                                                                      >  1:2108381/1‑100 (MQ=255)
          aCCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGaa                                                                                     >  1:366896/1‑100 (MQ=255)
            cTGATCCAAGGTAAAACCGCGAAACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCTGCTGATGgag                                                                                       <  1:1371682/96‑1 (MQ=255)
             tGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGaaaa                                                                                   >  2:2069695/1‑99 (MQ=255)
              gATCCCAGTTAAAACCGCTACACACATTCTGGTCGTAGAGCGCGCTTATCCGACGACTTACGAACGCTTTACCTCCATCAGCGGGCCGATagaggtaaag                                                                                 <  1:1568433/100‑7 (MQ=16)
              gATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGaaaa                                                                                   >  1:24690/1‑98 (MQ=255)
                tCCCAGGTAAACCCGCGCCACACATTATGGTCGGAGAGCGCGATTAGCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg                                                                               <  1:2071863/100‑1 (MQ=255)
                tCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg                                                                               <  1:329785/100‑1 (MQ=255)
                 cccAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt                                                                              >  2:2062688/1‑100 (MQ=255)
                  ccAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTa                                                                             >  2:2062259/1‑100 (MQ=255)
                  ccAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGCAAATCgg                                                                               >  2:1212186/1‑98 (MQ=255)
                   cAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                                                            >  2:203130/1‑100 (MQ=255)
                     ggTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg                                                                          >  1:403767/1‑100 (MQ=255)
                       tAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGc                                                                        >  1:1725280/1‑100 (MQ=255)
                         aaaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCgg                                                                      <  1:1990872/100‑1 (MQ=255)
                          aaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt                                                                     <  1:533622/100‑1 (MQ=255)
                           aCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTa                                                                    <  1:144548/100‑1 (MQ=255)
                           aCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTa                                                                    >  1:309647/1‑100 (MQ=255)
                           aCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAAGCGGTAATTTCGGTa                                                                    >  1:1218391/1‑100 (MQ=255)
                            ccGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTa                                                                    <  2:1595611/99‑1 (MQ=255)
                               cgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAgg                                                                <  2:1208720/100‑1 (MQ=255)
                               cgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTaaatg                                                                <  1:1520460/100‑3 (MQ=255)
                                 ccACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGTGTAAAggg                                                               >  2:1424071/1‑98 (MQ=255)
                                 ccACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCGTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGa                                                              <  2:662168/100‑1 (MQ=255)
                                    cacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTg                                                           <  1:1200428/100‑1 (MQ=255)
                                    cacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTg                                                           >  1:1877857/1‑100 (MQ=255)
                                      caTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGAtt                                                            >  1:1820575/1‑97 (MQ=255)
                                      caTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGcc                                                         >  1:1819206/1‑100 (MQ=255)
                                         tATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                                                      <  1:1357625/100‑1 (MQ=255)
                                         tATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                                                      <  2:1818305/100‑1 (MQ=255)
                                          aTGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGa                                                     <  1:1925403/100‑1 (MQ=255)
                                             gTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAaca                                                  >  1:328209/1‑100 (MQ=255)
                                                gTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAccc                                               <  2:2051885/100‑1 (MQ=255)
                                                 tAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCa                                              >  2:1397618/1‑100 (MQ=255)
                                                   gagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCaga                                            >  1:647654/1‑100 (MQ=255)
                                                   gagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCaga                                            <  2:1100616/100‑1 (MQ=255)
                                                     gcgcGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGc                                          >  1:2002759/1‑100 (MQ=255)
                                                       gcgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga                                        <  1:1499744/100‑1 (MQ=255)
                                                        cgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgag                                       >  2:1394289/1‑100 (MQ=255)
                                                        cgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgag                                       >  1:581571/1‑100 (MQ=255)
                                                          aTTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGCTTGCCTGGAACACCCAGAGCGAAAt                                     >  1:33239/1‑100 (MQ=255)
                                                          aTTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGAt                                     <  1:1686193/100‑1 (MQ=255)
                                                           ttATCCGGCGATTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATg                                    <  1:709299/100‑1 (MQ=255)
                                                            tATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATTGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGAt                                     >  2:903009/1‑98 (MQ=255)
                                                               ccGGCGACTTACGAACGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATc                                >  1:759930/1‑100 (MQ=255)
                                                               ccGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATTGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACATCCAGAGCGAGATGGATc                                <  2:573827/100‑1 (MQ=255)
                                                               ccGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATc                                <  2:573757/100‑1 (MQ=255)
                                                                    gACTTACGAACGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATctgctg                           <  1:205031/100‑1 (MQ=255)
                                                                    gACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATctgctg                           <  1:1481129/100‑1 (MQ=255)
                                                                     aCTTACGAACGCTTTACCTCGATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCtgctgc                          <  1:990073/100‑1 (MQ=255)
                                                                        tACGAACGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTa                       >  1:1039516/1‑100 (MQ=255)
                                                                        tACGAACGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTa                       >  2:1544925/1‑100 (MQ=255)
                                                                        tACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTa                       >  1:1297964/1‑100 (MQ=255)
                                                                        tACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTa                       <  1:482305/100‑1 (MQ=255)
                                                                          cGAACGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAg                     >  1:42501/1‑100 (MQ=255)
                                                                           gAACGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAg                     <  1:343604/99‑1 (MQ=255)
                                                                           gAACGCTTTACCTCTATCGGCCCGCTTATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGc                    <  2:1986840/100‑1 (MQ=255)
                                                                           gAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGc                    <  1:354223/100‑1 (MQ=255)
                                                                            aaCGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCt                   >  1:630610/1‑100 (MQ=255)
                                                                             aCGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTc                  >  1:1878288/1‑100 (MQ=255)
                                                                              cGCTTTACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCa                 <  2:257021/100‑1 (MQ=255)
                                                                                 tttACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACt              <  1:530176/100‑1 (MQ=255)
                                                                                   tACCTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTac            >  1:1767752/1‑100 (MQ=255)
                                                                                    aCCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCATAGCGAGATGTATCTGCTGCGTAAGCTCAAc               <  2:1259740/96‑1 (MQ=255)
                                                                                    aCCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTaca           <  2:309648/100‑1 (MQ=255)
                                                                                     cgatcttCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTacac          <  2:528063/94‑1 (MQ=255)
                                                                                     ccTCTATCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTacac          <  1:1015412/100‑1 (MQ=255)
                                                                                     ccTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTacac          >  1:2097707/1‑100 (MQ=255)
                                                                                      ctctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAcc         >  1:421384/1‑100 (MQ=255)
                                                                                           tCGGCCTGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTacacagatcg    >  1:528062/1‑94 (MQ=255)
                                                                                           tCGGCCCGCTGATGGAGAAAAACGGAAATGGCGGAAAAGGGATTGCCTGGAACACAAAGAACGAGATGGATCTGCTCCGCAAGCTCAACAACACCaa       >  2:729739/1‑97 (MQ=255)
                                                                                             ggCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCgg  >  1:1456965/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGG  >  NC_000913/1282480‑1282672

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