Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 575,007 G→A 7.0% G15D (GGC→GAC)  fdoG → formate dehydrogenase‑O major subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947575,0070GA 7.0% 165.0 / 2.2 57G15D (GGC→GAC) fdoGformate dehydrogenase‑O major subunit
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base A (2/2);  total (28/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 = 3.11e-01

GACTAGTCTTAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTTCGGTCGGCTGCGGCCTGATC  >  NC_002947/574863‑575146
                                                                                                                                                  |                                                                                                                                           
gACTAGTCTTAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCt                                                                                                                                          <  2:230565/148‑1 (MQ=255)
 aCTAGTCTTAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAACAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCtt                                                                                                                                         <  2:48857/148‑1 (MQ=255)
     gTCTTAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAgg                                                                                                                                    >  5:114595/1‑149 (MQ=255)
      tCTTAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAgg                                                                                                                                    >  4:112942/1‑148 (MQ=255)
        ttAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTc                                                                                                                                 >  5:245396/1‑149 (MQ=255)
           aCCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGa                                                                                                                               <  4:83448/148‑1 (MQ=255)
              tttGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGc                                                                                                                             <  8:87030/147‑1 (MQ=255)
              tttGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGc                                                                                                                             >  7:87030/1‑147 (MQ=255)
              tttGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCt                                                                                                                           >  6:135157/1‑149 (MQ=255)
                tGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCtt                                                                                                                          <  1:110864/148‑1 (MQ=255)
                      ccATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGt                                                                                                                   <  6:245396/149‑1 (MQ=255)
                       cATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCTCGACCACCACCTGATCAGGAGTGCACGATGGACCTCAACC‑‑GTCGGCAGTTCTTCAAGCTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGtt                                                                                                                  >  3:171182/1‑149 (MQ=255)
                           aaTTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTggg                                                                                                               <  6:114595/148‑1 (MQ=255)
                           aaTTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGc                                                                                                              <  4:228235/149‑1 (MQ=255)
                            aTTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCa                                                                                                             >  5:120116/1‑149 (MQ=255)
                                  cgccACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCATGCGGCGTTGGGCATGGccc                                                                                                       >  2:74384/1‑149 (MQ=255)
                                           gACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGa                                                                                                                               <  5:33459/116‑1 (MQ=255)
                                           gACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGa                                                                                                                               >  6:33459/1‑116 (MQ=255)
                                                ccAAGAGGCGCCGGCAATGTTCCGCGCGCCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTGTCGGCAGTTATTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCAt                                                                                           >  1:215929/1‑149 (MQ=255)
                                                 cAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCg                                                                                        >  7:102245/1‑149 (MQ=255)
                                                  aaGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGc                                                                                       <  1:105347/149‑1 (MQ=255)
                                                         gcCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCCACGCCGGCATTCGCCGAGCAg                                                                                >  5:5097/1‑149 (MQ=255)
                                                            ggCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGc                                                                                                         <  6:165729/121‑1 (MQ=255)
                                                            ggCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTAAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGc                                                                                                         >  5:165729/1‑121 (MQ=255)
                                                                        cgcgTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTAGGAGGCTCGAGCCTTGCGGCGTAGGGCATGGCCCCACCGCCCGCATTCGCCCAGCAGGTGCGTCCCTTCAAg                                                                 >  8:135740/1‑149 (MQ=255)
                                                                           gTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTg                                                              <  6:233876/149‑1 (MQ=255)
                                                                               accaccTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGgcgc                                                          <  6:120116/149‑1 (MQ=255)
                                                                               accaccTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGgcg                                                           >  6:95546/1‑148 (MQ=255)
                                                                               accaccTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCCGCGTTGGTCATGGCCCCCACGCCGGCATTCGCCGAGCAGGTGAGTCACTTAAAGCTGgcg                                                           >  8:230599/1‑148 (MQ=255)
                                                                                     tGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTc                                                                         >  2:98720/1‑128 (MQ=255)
                                                                                     tGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTc                                                                         <  1:98720/128‑1 (MQ=255)
                                                                                          aGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCTCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCAATTCAAGCTGGAGCATACCAAAGa                                                >  6:270597/1‑148 (MQ=255)
                                                                                              gTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGa                                                <  5:210927/144‑1 (MQ=255)
                                                                                              gTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGa                                                >  6:210927/1‑144 (MQ=255)
                                                                                                 cACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGt                                         >  5:226627/1‑148 (MQ=255)
                                                                                                 cACCATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTCGGCATGGCCCCAACGCCGGCATTCCCCGAGCACCTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGt                                         >  5:226625/1‑148 (MQ=255)
                                                                                                     aTGGACCTGAACC‑‑GTCGGGAGTTCTTCTAGGTCGCCGCCGTGGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAagac                                    <  6:226625/149‑1 (MQ=255)
                                                                                                     aTGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAacac                                    <  6:226627/149‑1 (MQ=255)
                                                                                                         aCCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGc                                <  5:135157/149‑1 (MQ=255)
                                                                                                           cTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCCCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGccc                              >  1:24000/1‑149 (MQ=255)
                                                                                                              aaCC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGcccc                             <  3:112942/147‑1 (MQ=255)
                                                                                                              aaCC‑‑GTCGGCAGTTCTTCAAGGTCGACGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATTGACCCAACGCCGTCATTCGCCGAGCAGGTGCGTCACTTCAATCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTa                           <  5:270597/149‑1 (MQ=255)
                                                                                                                     tCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAAcc                                            <  4:126246/127‑1 (MQ=255)
                                                                                                                     tCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAAcc                                            >  3:126246/1‑127 (MQ=255)
                                                                                                                         cAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTtcggt                  <  2:215929/149‑1 (MQ=255)
                                                                                                                              ctccAAGGCCGACGGCGGCGACCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGccc                              >  7:719/4‑132 (MQ=255)
                                                                                                                              ctccAAGGCCGACGGCGGCGACCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGccc                              <  8:719/129‑1 (MQ=255)
                                                                                                                               ttcAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTTCGGTCGGCTg            >  2:12402/1‑149 (MQ=255)
                                                                                                                               tccAAGGCCGACGGCGGCGACCTTGAAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCaa                                                                  >  6:192858/3‑95 (MQ=37)
                                                                                                                               tccAAGGCCGACGGCGGCGACCTTGAAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCaa                                                                  <  5:192858/93‑1 (MQ=37)
                                                                                                                                 cAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAg                                                                 >  4:47060/1‑94 (MQ=255)
                                                                                                                                 cAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAg                                                                 <  3:47060/94‑1 (MQ=255)
                                                                                                                                 cAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTTCGGTCGGCTGCg          >  5:40525/1‑149 (MQ=255)
                                                                                                                                   aGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTTCGGTCGGCTGCGGc        <  2:24000/149‑1 (MQ=255)
                                                                                                                                   aGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTTCGGTCGGCTGCGGc        <  3:173174/149‑1 (MQ=255)
                                                                                                                                     gTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTTCGGTCGGCTGCGGCCt      <  6:40525/149‑1 (MQ=255)
                                                                                                                                         ccgccgGCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTTCGGTCGGCTGCGGCCTGATc  <  7:230599/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
GACTAGTCTTAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACC‑‑GTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTTCGGTCGGCTGCGGCCTGATC  >  NC_002947/574863‑575146

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: