Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,711,857 G→A 8.6% E232E (GAG→GAA ubiE → ubiquinone/menaquinone biosynthesis C‑methyltransferase UbiE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,711,8570GA 8.6% 162.3 / 4.3 58E232E (GAG→GAAubiEubiquinone/menaquinone biosynthesis C‑methyltransferase UbiE
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base A (3/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.49e-01

TCCAAGGCCTACGACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGC  >  NC_002947/5711714‑5712002
                                                                                                                                               |                                                                                                                                                  
tCCAAGGCCTACGACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCgg                                                                                                                                               >  5:137702/1‑149 (MQ=255)
  cAAGGCCTACGACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAAGCCGGc                                                                                                                                              >  3:194811/1‑148 (MQ=255)
   aaGGCCTACGACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCtt                                                                                                                                            <  6:137702/149‑1 (MQ=255)
      gCCTACGACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGCCCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGa                                                                                                                                         >  5:85799/1‑149 (MQ=255)
       ccTACGACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGAc                                                                                                                                        <  4:79464/149‑1 (MQ=255)
          acgacgCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACcgc                                                                                                                                     >  1:129874/1‑149 (MQ=255)
           cgacgCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACcgc                                                                                                                                     <  6:122833/148‑1 (MQ=255)
            gacgCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAAGCCGGCCTCGAccatca                                                                                                                                   <  4:25346/149‑5 (MQ=255)
             acgCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTc                                                                                                                                  <  4:277462/149‑1 (MQ=255)
             acgCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAAGCCGGCCTCGAccatcat                                                                                                                                  <  4:194811/149‑6 (MQ=255)
              cgCCTACTCGTTCGCCTTCATGCCGGTGGCTGGCAAGCTGATCACCTATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAACCCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTct                                                                                                                                 <  6:309057/149‑2 (MQ=255)
               gCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCAc                                                                                                                                <  3:128138/149‑1 (MQ=255)
                 cTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTc                                                                                                                                           <  7:270844/136‑1 (MQ=255)
                 cTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTc                                                                                                                                           >  8:270844/1‑136 (MQ=255)
                     tCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACca                                                                                                                          <  3:125160/149‑1 (MQ=255)
                       gTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAACCCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGCCCGCGTCACCTAcccc                                                                                                                         >  6:55002/1‑146 (MQ=255)
                           gcactCATGCCGCTGGCTGGCAAGCTGATCACCACTCACGCCGAGACCTACCGTTACCTCCCCGAGTCGATCCGCATGCACCCGGACCAGGAACCCCTGATGAGCAGGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACAt                                                                                                                    <  6:215201/145‑1 (MQ=255)
                                  tGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGc                                                                                                            >  7:67748/1‑149 (MQ=255)
                                    ccGCTGGATGGCAAGCTGATCACCAATGACGCAGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATACACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCg                                                                                                           <  1:49788/148‑1 (MQ=255)
                                            tggcAAGCTGATCACCAATGACGCAGAGAGCTACCGTTACCTCGCCGAGTCGCTCCGCATGCACCCGGACCAGGAACCCCTGAAGAGCATCATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGAGGCATCCTTg                                                                                                  >  6:22096/1‑149 (MQ=255)
                                               cAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCg                                                                                               >  6:316378/1‑149 (MQ=255)
                                                 aGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTg                                                                                             <  2:56696/149‑1 (MQ=255)
                                                  gCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGc                                                                                            >  7:133993/1‑149 (MQ=255)
                                                         aCCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCgggg                                                                                     >  5:78101/1‑149 (MQ=255)
                                                         aCCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCggg                                                                                      <  6:297909/148‑1 (MQ=255)
                                                          ccAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCTGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCgggg                                                                                     >  2:287939/1‑148 (MQ=255)
                                                          ccAAAGACTCCAATAGCTACCCATACCTCGCCGACTCCAACCGCAAGCACCCGGACAAGGAAACCCTCAATACCAACATGCTCAAGGCCAGCTTCCACAGCGTCACCTACCACAACATGACCAGCGGCATAGATTCCGAACACCgcgc                                                                                     >  4:3092/1‑145 (MQ=255)
                                                                 cGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTTAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGc                                                                             <  8:211053/149‑1 (MQ=255)
                                                                  gCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGcc                                                                            >  2:93683/1‑149 (MQ=255)
                                                                    cGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCt                                                                          >  8:34536/1‑149 (MQ=255)
                                                                             ccGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCt                                                                 <  6:78101/149‑1 (MQ=255)
                                                                                     tCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGTAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCgg                                                                                       <  1:104405/119‑1 (MQ=255)
                                                                                     tCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCgg                                                                                       >  2:104405/1‑119 (MQ=255)
                                                                                        ccGAGTCGATCCGCATGCACCCGGACCAGTAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACAACAACATCTCCCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGgctg                                                       <  1:93683/149‑1 (MQ=255)
                                                                                               gATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAg                                               >  1:191389/1‑149 (MQ=255)
                                                                                                 tCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCg                                             >  2:174449/1‑149 (MQ=255)
                                                                                                      aTGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGaa                                        <  7:34536/149‑1 (MQ=255)
                                                                                                       tGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGaa                                        >  8:263054/1‑148 (MQ=255)
                                                                                                        gCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACa                                      >  7:260320/1‑149 (MQ=255)
                                                                                                          aCCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATg                                    >  5:317246/1‑149 (MQ=255)
                                                                                                          aCCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATg                                    >  2:158285/1‑149 (MQ=255)
                                                                                                          aCCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATg                                    >  2:179121/1‑149 (MQ=255)
                                                                                                                        aCCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTc                      <  1:174449/149‑1 (MQ=255)
                                                                                                                          ccTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCt                    >  2:268090/1‑149 (MQ=255)
                                                                                                                              aaGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTgc                  <  6:85799/147‑1 (MQ=255)
                                                                                                                                 accATGAGGGGCGAGGCCGGCTGCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTg             <  8:99510/147‑1 (MQ=255)
                                                                                                                                   cttgatgGGCGAGTGCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGAGCAAGCCCTGATGTGGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGa           <  5:22096/147‑1 (MQ=255)
                                                                                                                                   cATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGTGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGa           >  4:261075/1‑149 (MQ=255)
                                                                                                                                   cATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGa           >  7:203984/1‑149 (MQ=255)
                                                                                                                                        tgGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTgc                  <  6:287709/137‑1 (MQ=255)
                                                                                                                                        tgGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTgc                  >  5:287709/1‑137 (MQ=255)
                                                                                                                                        cggTCGAAGCCGGCCTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCa      >  7:279601/2‑149 (MQ=255)
                                                                                                                                         ggTCGAAGCCGGCCTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCAc     >  6:90261/1‑149 (MQ=255)
                                                                                                                                           tCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACgg   >  4:38846/1‑149 (MQ=255)
                                                                                                                                           tCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACgg   <  1:268090/149‑1 (MQ=255)
                                                                                                                                           tCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACgg   <  1:179121/149‑1 (MQ=255)
                                                                                                                                           tCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACgg   >  6:47818/1‑149 (MQ=255)
                                                                                                                                            cGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTGGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGc  <  5:55002/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TCCAAGGCCTACGACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGC  >  NC_002947/5711714‑5712002

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