Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 428,188 G→T R181L (CGC→CTC)  secD → SecYEG protein translocase auxillary subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913428,1880GT95.1% 272.6 / 3.2 82R181L (CGC→CTC) secDSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (41/37);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGCGGTA  >  NC_000913/428093‑428285
                                                                                               |                                                                                                 
aaaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGAt                                                                                               >  2:1900702/1‑100 (MQ=255)
aaaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGAt                                                                                               <  1:258568/100‑1 (MQ=255)
aaaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGAt                                                                                               >  2:889150/1‑100 (MQ=255)
aaaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGAt                                                                                               >  1:186334/1‑100 (MQ=255)
    atCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTa                                                                                           >  2:1488505/1‑100 (MQ=255)
       gATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAg                                                                                        <  1:613680/100‑1 (MQ=255)
       gATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAg                                                                                        >  2:445188/1‑100 (MQ=255)
        aTAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGc                                                                                       <  2:19981/100‑1 (MQ=255)
          aGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTc                                                                                     <  1:806026/100‑1 (MQ=255)
          aGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTc                                                                                     >  2:453894/1‑100 (MQ=255)
           tccTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCt                                                                                      >  2:1752003/2‑98 (MQ=255)
           gCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCg                                                                                    <  2:2292574/100‑1 (MQ=255)
           gCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCg                                                                                    >  2:310578/1‑100 (MQ=255)
            ccTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGt                                                                                   <  2:1821724/100‑1 (MQ=255)
            ccTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGt                                                                                   <  2:1978799/100‑1 (MQ=255)
            ccTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGt                                                                                   <  1:879887/100‑1 (MQ=255)
            ccTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGt                                                                                   >  2:576997/1‑100 (MQ=255)
             cTACGCAGTGACCCGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtg                                                                                  <  1:1891539/100‑1 (MQ=255)
              tACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtga                                                                                 >  1:1683737/1‑100 (MQ=255)
              tACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtga                                                                                 <  2:726036/100‑1 (MQ=255)
               aCGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtgat                                                                                >  2:935231/1‑100 (MQ=255)
                    gTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGc                                                                           >  2:1872295/1‑100 (MQ=255)
                       aCCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCAt                                                                        >  1:1393697/1‑100 (MQ=255)
                         cTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTg                                                                      <  2:2088012/100‑1 (MQ=255)
                          tGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGc                                                                     >  1:475872/1‑100 (MQ=255)
                           gcgcGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCg                                                                    >  2:273401/1‑100 (MQ=255)
                                aaaaGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATc                                                                >  2:1982704/1‑99 (MQ=255)
                                   aGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAg                                                            <  1:1017636/100‑1 (MQ=255)
                                   aGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCACCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAg                                                            <  2:2031916/100‑1 (MQ=255)
                                    gggCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGc                                                           <  1:453896/100‑1 (MQ=255)
                                      gCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCaa                                                         >  2:27585/1‑100 (MQ=255)
                                        aTCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAGGCCATTGCGTATCTGAGCAAgc                                                       >  1:1909039/1‑100 (MQ=255)
                                         tCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgcg                                                      <  2:121129/100‑1 (MQ=255)
                                         tCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgcg                                                      <  2:1280600/100‑1 (MQ=255)
                                          cccGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgcgc                                                     <  1:1488509/100‑1 (MQ=255)
                                               ataCCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATcc                                                <  2:1470436/100‑1 (MQ=255)
                                                taCCACTGTTCGTAAAGAAAACAACTACGGTCTGATCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCg                                               <  1:1393440/100‑1 (MQ=255)
                                                taCCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCg                                               >  2:1591442/1‑100 (MQ=255)
                                                 aCCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCgg                                              >  2:1231636/1‑100 (MQ=255)
                                                  ccACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGa                                             >  1:1932725/1‑100 (MQ=255)
                                                   cACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAc                                            >  1:2142770/1‑100 (MQ=255)
                                                   cACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAc                                            >  2:557611/1‑100 (MQ=255)
                                                    aCTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAcc                                           <  2:1908876/100‑1 (MQ=255)
                                                    aCTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAcc                                           >  2:2261595/1‑100 (MQ=255)
                                                    aCTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAcc                                           >  1:829245/1‑100 (MQ=255)
                                                     cTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCt                                          >  1:886730/1‑100 (MQ=255)
                                                      tccgatctaaGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTg                                         <  1:1842962/92‑1 (MQ=255)
                                                      tGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTg                                         >  2:2157469/1‑100 (MQ=255)
                                                       gTTCGTAAAGAAAACAACTACGGCCTGAGTATCACTTTCCTCGATGCGAAAGCTCGTGATGAAGCCTTTGCGTATCTGAGCAAGCGCCATCCGGACCTgg                                        <  2:1408638/100‑1 (MQ=255)
                                                        ttCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGt                                       >  2:962993/1‑100 (MQ=255)
                                                        ttCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGt                                       <  2:1393694/100‑1 (MQ=255)
                                                        ttCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGt                                       <  1:576998/100‑1 (MQ=255)
                                                         tCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTg                                      <  2:562355/100‑1 (MQ=255)
                                                         tCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTg                                      <  1:562737/100‑1 (MQ=255)
                                                         tCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTg                                      <  1:801338/100‑1 (MQ=255)
                                                            tAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAtt                                   >  2:264691/1‑100 (MQ=255)
                                                            naaaGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAtt                                   >  1:1131377/2‑100 (MQ=255)
                                                              aaGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTgagatcgga                                 >  2:1842957/1‑92 (MQ=255)
                                                               ttccgatctcTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagc                                <  1:389908/91‑1 (MQ=255)
                                                               aGAAAACAACTATGGACTGAGCATAACATTCCTAGATGCAAAAGCTCGTGATGAAGCAATCGCGTATCTAAGCAAGCGCAATCCAGACCTgagaaacc                                  >  2:2164354/1‑91 (MQ=255)
                                                                gAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagca                               <  2:1131375/100‑1 (MQ=255)
                                                                 aaaaCAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagcag                              >  1:1379247/1‑100 (MQ=255)
                                                                  aaaCAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagcagc                             >  1:269646/1‑100 (MQ=255)
                                                                    acaacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCa                           >  2:2015545/1‑100 (MQ=255)
                                                                     caacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAg                          <  1:2261603/100‑1 (MQ=255)
                                                                     caacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAg                          <  1:22883/100‑1 (MQ=255)
                                                                      aacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAgg                         >  1:307468/1‑100 (MQ=255)
                                                                        cTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagcagatcggaa                       >  2:389906/1‑93 (MQ=255)
                                                                        cTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgc                       >  2:1495788/1‑100 (MQ=255)
                                                                        cTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgc                       <  2:365032/100‑1 (MQ=255)
                                                                        cTACGCCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgc                       <  1:935233/100‑1 (MQ=255)
                                                                         tACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgca                      <  1:1982711/100‑1 (MQ=255)
                                                                         tACGGCCTGAGCATCACGTTCCTCGATGCTAAAGCGCGTGATGAAGCCATAGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgca                      <  1:868637/100‑1 (MQ=255)
                                                                           cGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcagc                    <  2:1955151/100‑1 (MQ=255)
                                                                              ccTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAAc                 <  2:829243/100‑1 (MQ=255)
                                                                                  aGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGc             >  2:795572/1‑100 (MQ=255)
                                                                                   gCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCt            >  2:2009590/1‑100 (MQ=255)
                                                                                       cACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcgt        >  2:740022/1‑100 (MQ=255)
                                                                                       cACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcgt        >  1:641221/1‑100 (MQ=255)
                                                                                         cTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcgtgc      >  2:433633/1‑100 (MQ=255)
                                                                                         cTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcgtgc      <  1:445190/100‑1 (MQ=255)
                                                                                             ccGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGCGGTa  <  2:1379244/100‑1 (MQ=255)
                                                                                               |                                                                                                 
AAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGCGGTA  >  NC_000913/428093‑428285

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