Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 635,840 G→T 11.6% E22* (GAA→TAA) ‡ 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,8400GT11.6% 143.8 / 6.3 51E22* (GAA→TAA) ‡mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base T (3/3);  total (25/26)
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

ATTACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCAT  >  NC_002947/635693‑635988
                                                                                                                                                    |                                                                                                                                                     
aTTACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAgg                                                                                                                                                       >  5:159197/1‑148 (MQ=255)
 ttACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGaa                                                                                                                                                     <  6:314499/149‑1 (MQ=255)
  tACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTCCCCTGCCTCGAGGT‑TTGGCGAGGTCGACGCTGCATATTCACGTTGGGGGTATTTGCGGGAGTTTGATGGGTTCGCTGGCGGTGCTGGCCAAGGacc                                                                                                                                                    >  4:183476/1‑147 (MQ=255)
     cAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTg                                                                                                                                                <  3:32418/149‑1 (MQ=255)
       ggggCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGc                                                                                                                                              >  5:271242/1‑149 (MQ=255)
         ggCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCC‑TTGGCCa                                                                                                                                            <  5:168568/149‑1 (MQ=255)
            gACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTATTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCAAAGGAAC‑TTGGCCAcc                                                                                                                                          >  1:88689/1‑149 (MQ=255)
                 aGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCC‑TTGGCCag                                                                                                                                           <  5:62615/142‑2 (MQ=255)
                 aGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACcgcg                                                                                                                                       <  6:112858/146‑1 (MQ=255)
                 aGCGAAACAGCGGGTGCGATTTGTCGCGACCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTATCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCC‑TTGGCCag                                                                                                                                           >  6:62615/1‑141 (MQ=255)
                  gCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCAc                                                                                                                                   >  6:310909/1‑149 (MQ=255)
                    gAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCg                                                                                                                                 >  6:283425/1‑149 (MQ=255)
                        cAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGctc                                                                                                                             <  6:163266/149‑1 (MQ=255)
                         aGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGctct                                                                                                                            >  1:69895/1‑149 (MQ=255)
                           cGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGa                                                                                                                          <  6:76040/149‑1 (MQ=255)
                            gggTGCGATTTGTCGCGCCCTGGGCTACCCTGTCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGAc                                                                                                                         >  1:263302/1‑149 (MQ=255)
                            gggTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGAc                                                                                                                         <  6:159197/149‑1 (MQ=255)
                            gggTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGAc                                                                                                                         <  5:191829/149‑1 (MQ=255)
                               tGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGcg                                                                                                                      <  4:108704/149‑2 (MQ=255)
                                        tCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTAt                                                                                                             >  5:109718/1‑149 (MQ=255)
                                               cTGGGCTACCCTGTCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGa                                                                                                      <  2:263302/149‑1 (MQ=255)
                                               cTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGa                                                                                                      <  3:5404/149‑1 (MQ=255)
                                                  ggCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGa                                                                                                   <  6:109718/149‑1 (MQ=255)
                                                  ggCTACCCTGCCTCGATGT‑TTTGCGAGGTCGAACATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGa                                                                                                   <  5:177716/149‑1 (MQ=255)
                                                      aCCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCTGTATTTGCTGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGTCTCTGACGCCAACGTCTATCCCCCGATGAGCa                                                                                                <  5:310909/148‑1 (MQ=255)
                                                          tGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAg                                                                                           >  2:308302/1‑149 (MQ=255)
                                                            ccTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGc                                                                                          <  2:69895/148‑1 (MQ=255)
                                                                gATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCg                                                                                       >  5:78403/1‑147 (MQ=255)
                                                                          cGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAggg                                                                             <  8:141063/148‑1 (MQ=255)
                                                                           gAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCa                                                                                >  1:276896/1‑144 (MQ=255)
                                                                           gAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCa                                                                                <  2:276896/144‑1 (MQ=255)
                                                                             ggTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATc                                                                         <  6:271242/149‑1 (MQ=255)
                                                                                             tCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGaa                                                                                     >  4:261159/1‑121 (MQ=255)
                                                                                             tCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGaa                                                                                     <  3:261159/121‑1 (MQ=255)
                                                                                                caTTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTAt                                                       >  7:84967/1‑148 (MQ=255)
                                                                                                   tctcGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATACCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGAcc                                                   >  7:124818/1‑149 (MQ=255)
                                                                                                     tcGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACGTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGAccc                                                  >  1:244205/1‑149 (MQ=255)
                                                                                                     tcGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCg                                                 >  5:144440/1‑149 (MQ=255)
                                                                                                       ggTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGc                                               <  6:144440/149‑1 (MQ=255)
                                                                                                         tATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGcc                                              >  4:41741/1‑148 (MQ=255)
                                                                                                                    cTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCg                                  <  8:284038/149‑1 (MQ=255)
                                                                                                                         aTGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGGACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcgc                              <  6:78403/148‑1 (MQ=255)
                                                                                                                                cGCTGGCGGTGCCGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtg                      >  5:285066/1‑149 (MQ=255)
                                                                                                                                     gCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGGCGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTggtggt                 >  5:291189/1‑149 (MQ=255)
                                                                                                                                     gCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggtg                   <  2:88689/147‑1 (MQ=255)
                                                                                                                                        gcggtGGCCAAGTTCC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATc                                    <  2:57015/123‑1 (MQ=255)
                                                                                                                                        gcggtGGCCAAGTTCC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATc                                    >  1:57015/5‑127 (MQ=255)
                                                                                                                                            tGGCCAAGTTCC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCgg           >  4:297791/1‑148 (MQ=255)
                                                                                                                                             ggCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGc          >  7:43030/1‑148 (MQ=255)
                                                                                                                                              gCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGc          <  1:308302/147‑1 (MQ=255)
                                                                                                                                                cAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACACAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACg      >  2:260433/1‑149 (MQ=255)
                                                                                                                                                  aGGAAC‑TTGGCCCCCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGCACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGc     >  6:280691/1‑148 (MQ=255)
                                                                                                                                                   ggAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCa   >  5:121348/1‑149 (MQ=255)
                                                                                                                                                    gAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCt                       <  6:251637/128‑1 (MQ=255)
                                                                                                                                                    gAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCt                       >  5:251637/1‑128 (MQ=255)
                                                                                                                                                    gAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCAt  <  3:183476/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCAT  >  NC_002947/635693‑635988

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