Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 635,842 A→C 6.9% E22D (GAA→GAC) ‡ mpl → UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947635,8420AC 6.9% 174.7 / 2.6 56E22D (GAA→GAC) ‡mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base C (2/2);  total (27/29)
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.00e+00

CCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGT  >  NC_002947/635697‑635990
                                                                                                                                                 |                                                                                                                                                     
ccAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACtt                                                                                                                                                    >  1:170530/1‑149 (MQ=255)
 cAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTg                                                                                                                                                   <  5:83811/149‑1 (MQ=255)
 cAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTg                                                                                                                                                   >  1:77171/1‑149 (MQ=255)
       cGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCAcc                                                                                                                                             <  4:24415/149‑1 (MQ=255)
        gACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCCTTGGCCagca                                                                                                                                            <  4:90427/149‑4 (MQ=255)
        gACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACcg                                                                                                                                            >  8:61542/1‑149 (MQ=255)
             aGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCa                                                                                                                                       >  6:125384/1‑149 (MQ=255)
                gAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCg                                                                                                                                    >  2:80433/1‑149 (MQ=255)
                 aaaCAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCgg                                                                                                                                   <  3:203548/149‑1 (MQ=255)
                 aaaCAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCgg                                                                                                                                   <  6:192411/149‑1 (MQ=255)
                       cGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGa                                                                                                                             >  2:71618/1‑149 (MQ=255)
                       cGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGa                                                                                                                             <  2:170530/149‑1 (MQ=255)
                        gggTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGAc                                                                                                                            >  3:89987/1‑149 (MQ=255)
                                       cgcCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATccc                                                                                                             <  2:246013/149‑1 (MQ=255)
                                          ccTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATccccc                                                                                                           >  3:69729/1‑148 (MQ=255)
                                               gCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAg                                                                                                     >  3:67667/1‑149 (MQ=255)
                                                 tACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGAACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCATCGTCTATCGCCCGATGAGCa                                                                                                   >  8:225554/1‑149 (MQ=255)
                                                   cccTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCAcc                                                                                                 <  6:254534/149‑1 (MQ=255)
                                                    ccTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGa                                                                                                      >  6:244993/1‑143 (MQ=255)
                                                    ccTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGa                                                                                                      <  5:244993/143‑1 (MQ=255)
                                                            gATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGaa                                                                                        <  2:236837/149‑1 (MQ=255)
                                                            gATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGaa                                                                                        <  2:77171/149‑1 (MQ=255)
                                                               gTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATc                                                                                                               >  3:235690/1‑123 (MQ=255)
                                                               gTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATc                                                                                                               <  4:235690/123‑1 (MQ=255)
                                                                    gCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAgg                                                                                 <  7:61542/148‑1 (MQ=255)
                                                                    gCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAgg                                                                                 <  1:159227/148‑1 (MQ=255)
                                                                    gCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAgg                                                                                 <  1:121853/148‑1 (MQ=255)
                                                                                     tatTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAg                                                              >  1:138184/1‑149 (MQ=255)
                                                                                      atTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCAGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAgg                                                             >  5:122143/1‑149 (MQ=255)
                                                                                       ttCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGc                                                            >  4:28421/1‑149 (MQ=255)
                                                                                          acaTTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTAt                                                         <  2:138184/149‑1 (MQ=255)
                                                                                            aTTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGa                                                       <  4:69729/149‑1 (MQ=255)
                                                                                             ttCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGAc                                                      >  3:123984/1‑149 (MQ=255)
                                                                                              tctcGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGAcc                                                     >  2:293676/1‑149 (MQ=255)
                                                                                                  ggTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCgg                                                  <  5:125384/148‑1 (MQ=255)
                                                                                                           ggCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGa                                        <  3:28421/149‑1 (MQ=255)
                                                                                                                tttCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGCCTATGACCCGGCCCAGCTGGATCCgc                                   >  5:271229/1‑148 (MQ=255)
                                                                                                                   cATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcg                                 <  4:67667/148‑1 (MQ=255)
                                                                                                                    aTGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcgc                                <  4:89987/148‑1 (MQ=255)
                                                                                                                     tGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGCATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCa                              <  1:71618/149‑1 (MQ=255)
                                                                                                                      gggTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTCTCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcgc                                >  5:86171/1‑146 (MQ=255)
                                                                                                                      gggTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcgc                                >  6:19394/1‑146 (MQ=255)
                                                                                                                      gggTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCa                              <  5:19394/148‑1 (MQ=255)
                                                                                                                       ggTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAg                             <  3:54718/148‑1 (MQ=255)
                                                                                                                        gTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGc                                                                                      >  6:105644/1‑91 (MQ=255)
                                                                                                                        gTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGc                                                                                      <  5:105644/91‑1 (MQ=255)
                                                                                                                            gCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtgg                       >  5:91277/1‑149 (MQ=255)
                                                                                                                               ggCGGGGCTGGCCAAGGAACTTGGCCACCTCGTCAACGGCTCTGCCGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCCAGTCTATGACCCGGCCCATCTGGCTCCGGCGCCAGACCtggtgg                    <  2:247023/149‑1 (MQ=255)
                                                                                                                                gCGGTGGTGGCCAAGGAACTTGGCCAACGGGTCACCGGGTCTGACCCCAAACGCTATCCCCCGATGAGCACCCATCTCGAAGCCCATGGTATCTAAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggtg                     <  1:80433/148‑1 (MQ=255)
                                                                                                                                  ggTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCa                 >  6:142947/1‑149 (MQ=255)
                                                                                                                                   gcggtGGCCAAGTTCCTTGGCCAGCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgc                                  <  5:280572/127‑1 (MQ=255)
                                                                                                                                   gcggtGGCCAAGTTCCTTGGCCAGCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgc                                  >  6:280572/5‑131 (MQ=255)
                                                                                                                                     ggtGGCCAAGTTCCTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGTCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCg              >  4:243231/3‑149 (MQ=255)
                                                                                                                                         gCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGc            <  5:142947/147‑1 (MQ=255)
                                                                                                                                         gCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCaa          >  1:278580/1‑149 (MQ=255)
                                                                                                                                             aGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCAACGGCAACGc       >  4:97903/1‑148 (MQ=255)
                                                                                                                                                 aCTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGt  >  6:129980/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTTTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATC‑GAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCATGT  >  NC_002947/635697‑635990

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