Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,079,196 C→T L290L (CTG→CTA cbpA ← curved DNA‑binding protein, DnaJ homologue that functions as a co‑chaperone of DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,079,1960CT99.0% 808.9 / ‑3.0 209L290L (CTG→CTAcbpAcurved DNA‑binding protein, DnaJ homologue that functions as a co‑chaperone of DnaK
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (110/97);  total (111/98)
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.74e-01

Only 100 of 220 total aligned reads displayed.

AAAGTCACCGTAACATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCCAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGA  >  REL606/1079123‑1079268
                                                                         |                                                                        
aaaGTCACCGTAACATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAg                                                                        >  1:1665984/1‑76 (MQ=255)
aaaGTCACCGTAACATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAg                                                                        <  2:5849125/76‑1 (MQ=255)
 aaGTCACCGTAACATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGt                                                                       >  2:4557392/1‑76 (MQ=255)
   gTCACCGTAACATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTtg                                                                     <  2:7602364/76‑1 (MQ=255)
    tCACCGTAACATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTtgc                                                                    >  2:7346258/1‑76 (MQ=255)
    tCACCGTAACATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTtgc                                                                    >  2:4389598/1‑76 (MQ=255)
      aCCGTAACATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTtgctg                                                                  >  1:4554282/1‑76 (MQ=255)
       ccGTAACATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTtgctgc                                                                 <  1:6339418/76‑1 (MQ=255)
       ccGTAACATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTtgctgc                                                                 <  1:1434253/76‑1 (MQ=255)
         gTAACATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCca                                                               >  1:6701187/1‑76 (MQ=255)
          tAACATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCcac                                                              <  1:1933918/76‑1 (MQ=255)
          tAACATTAGCCATTATGCTTCCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCcac                                                              <  2:8299093/76‑1 (MQ=255)
           aaCATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCcaca                                                             >  1:6957506/1‑76 (MQ=255)
           aaCATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCcaca                                                             >  2:7741355/1‑76 (MQ=255)
            aCATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAg                                                            <  2:598797/76‑1 (MQ=255)
            aCATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAg                                                            <  2:3877588/76‑1 (MQ=255)
             cgttAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAgc                                                           >  1:2614462/3‑76 (MQ=255)
              aTTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAgcg                                                          <  1:2958914/76‑1 (MQ=255)
               ttAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAgcgc                                                         <  2:2836328/76‑1 (MQ=255)
                 aGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCgg                                                       <  2:1681543/76‑1 (MQ=255)
                 aGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCgg                                                       <  1:7547686/76‑1 (MQ=255)
                    cATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAg                                                    >  2:8175544/1‑76 (MQ=255)
                    cATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAg                                                    >  2:302027/1‑76 (MQ=255)
                     aTTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAgg                                                   >  2:7847013/1‑75 (MQ=255)
                     aTTATGCTTTCCCCCAATCTTTACGTGGAGCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAgt                                                   <  2:1216071/76‑1 (MQ=255)
                      ttATGCTTTCCCCCAATCTTTACGTGAATCAAAAGACGACTGGGCGTCTGGTAGTTGCTGCCACAGCGCGGCAgtg                                                  <  2:7919828/76‑1 (MQ=255)
                       tATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAgtgt                                                 <  2:3361511/76‑1 (MQ=255)
                         tGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGttt                                               <  1:7953315/76‑1 (MQ=255)
                          tcTTCCCCCCAATCTTTACGGGGAACAAAAGACGACTGGGCGTCTGCCAGTTGCTGCCACAGCGCGGCAGTGtttt                                              <  1:3307158/75‑1 (MQ=255)
                          gCTTTCCCCCATTCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCAACAGCGCCGAAGTGtttt                                              >  2:3971205/1‑76 (MQ=255)
                          gCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGtttt                                              <  1:5128985/76‑1 (MQ=255)
                          gCTTTCCCCAAATTTTTACGTGGATCAAAAGACGACGGGGCGTCTGCTATTTGCTGCCACAGCGCGGCAGTGtttt                                              <  2:8422756/76‑1 (MQ=255)
                           cTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTc                                             >  1:1094500/1‑76 (MQ=255)
                           cTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGCGTTTTc                                             >  2:6095239/1‑76 (MQ=255)
                            tttCCCCCAATCTTTCCGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCa                                            <  1:5606305/76‑1 (MQ=255)
                            tttCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCa                                            >  2:2541843/1‑76 (MQ=255)
                            tttCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCa                                            >  1:1547026/1‑76 (MQ=255)
                             ttCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTTCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCAt                                           <  2:786987/76‑1 (MQ=255)
                             ttCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCAt                                           >  2:3822183/1‑76 (MQ=255)
                             ttCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCAt                                           <  1:6264288/76‑1 (MQ=255)
                             ttCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGGTTTCAt                                           >  2:41487/1‑76 (MQ=255)
                              tCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATc                                          >  2:4729418/1‑76 (MQ=255)
                              tCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATc                                          >  2:8007837/1‑76 (MQ=255)
                               cccccAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATcc                                         >  1:880627/1‑76 (MQ=255)
                               cccccAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGGGTTTTCATcc                                         >  1:4847882/1‑76 (MQ=255)
                                ccccAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCg                                        >  2:8126058/1‑76 (MQ=255)
                                 cccAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATTCgg                                       >  1:5701556/1‑76 (MQ=255)
                                   cAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGtt                                     <  2:7025839/76‑1 (MQ=255)
                                    aaTCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGttt                                    <  2:5286599/76‑1 (MQ=255)
                                    aaTCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGttt                                    <  1:8383142/76‑1 (MQ=255)
                                     aTCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTc                                   >  2:2148631/1‑76 (MQ=255)
                                       cTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTcgg                                 >  2:3173065/1‑76 (MQ=255)
                                       cTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTcgg                                 >  2:1334820/1‑76 (MQ=255)
                                       cTTTACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTcgg                                 >  1:2938516/1‑76 (MQ=255)
                                        tttACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTcggc                                >  1:3079614/1‑76 (MQ=255)
                                        tttACGGGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTcggc                                <  1:4151131/76‑1 (MQ=255)
                                         ttACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTcggcg                               >  1:1853015/1‑76 (MQ=255)
                                         ttACGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTcggcg                               >  1:8176945/1‑76 (MQ=255)
                                           aCGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCggcggc                             <  2:5537630/76‑1 (MQ=255)
                                           aCGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCggcggc                             >  1:545074/1‑76 (MQ=255)
                                            cGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCggcggcc                            >  2:4772835/1‑75 (MQ=255)
                                            cGTGGATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCa                            >  2:1826586/1‑76 (MQ=255)
                                              gggATCAAAAGACGCCTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATc                          <  2:2434165/75‑1 (MQ=255)
                                              gggATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATc                          <  1:5955880/75‑1 (MQ=255)
                                                gATCAAAAGGCGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCAc                        <  1:4019674/76‑1 (MQ=255)
                                                gATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCAc                        <  2:311568/76‑1 (MQ=255)
                                                gATCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCAc                        >  1:30026/1‑76 (MQ=255)
                                                  tCAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGGCTCACGa                      >  2:4853306/1‑76 (MQ=255)
                                                   cAAAAGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGAt                     <  1:987729/76‑1 (MQ=255)
                                                    aaagGACGCCCGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGAtt                    <  2:4266779/76‑1 (MQ=255)
                                                    aaaaGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGAtt                    >  2:3703326/1‑76 (MQ=255)
                                                     aaaGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGAttt                   >  1:2025377/1‑76 (MQ=255)
                                                     aaaGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCCGCGGCATCACGAttt                   >  2:7538712/1‑76 (MQ=255)
                                                      aaGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGAtttt                  <  1:4287786/76‑1 (MQ=255)
                                                       aGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTc                 <  1:4400404/76‑1 (MQ=255)
                                                       aGACGACTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTc                 >  1:6178378/1‑76 (MQ=255)
                                                          cgacTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCCGt              >  2:2279484/1‑76 (MQ=255)
                                                          cgacTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGt              <  2:1286486/76‑1 (MQ=255)
                                                          cgacTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGt              <  2:659736/76‑1 (MQ=255)
                                                          cgacTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGt              <  1:2633210/76‑1 (MQ=255)
                                                            acTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCCTCACGATTTTCAGTAc            >  2:1770159/1‑76 (MQ=255)
                                                            acTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTAc            >  1:2747211/1‑76 (MQ=255)
                                                             cTGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCCTCACGATTTTCAGTAcc           >  2:6647836/1‑76 (MQ=255)
                                                             cTGGGCGTCTGCTAGGTGCTGCCACAGCGCGGCAGTGATTTCATCCGGGTTCGGCGGCATCACGATTTTCAGTAcc           >  2:6405361/1‑76 (MQ=255)
                                                             cTGGGCGTCTGCCAGTCGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTAca           <  2:2342968/76‑2 (MQ=255)
                                                              tGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCg          <  1:6836858/76‑1 (MQ=255)
                                                              tGGGCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGGACCg          >  1:489645/1‑76 (MQ=255)
                                                              tGGGCGTCTGCCAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCg          >  2:1606817/1‑76 (MQ=255)
                                                              ggggCTTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCg          <  2:8859737/75‑1 (MQ=255)
                                                                 gCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCAt       <  2:6696289/76‑1 (MQ=255)
                                                                 gCGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCAt       <  2:4343377/76‑1 (MQ=255)
                                                                  cGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCCTCACGATTTTCAGTACCCCATa      >  2:6364635/1‑76 (MQ=255)
                                                                  cGTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATa      <  2:2926897/76‑1 (MQ=255)
                                                                   gTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCCTCACGATTTTCCGTAACGCacac     >  2:4522403/1‑76 (MQ=255)
                                                                   gTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGGTTTTCCGGACCGCATAc     >  2:1461711/1‑76 (MQ=255)
                                                                   gTCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACCATTTTCAGTAACCCATaa     >  2:5680595/1‑75 (MQ=255)
                                                                    tCTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACa    <  1:1795908/76‑1 (MQ=255)
                                                                     cTGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAg   >  1:1199567/1‑76 (MQ=255)
                                                                      tGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGa  <  1:2291315/76‑1 (MQ=255)
                                                                      tGCTAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGGTTAGGCCGCATCACGATTTTCAGGACCGCATACAGa  >  2:7494495/1‑76 (MQ=255)
                                                                         |                                                                        
AAAGTCACCGTAACATTAGCCATTATGCTTTCCCCCAATCTTTACGTGGATCAAAAGACGACTGGGCGTCTGCCAGTTGCTGCCACAGCGCGGCAGTGTTTTCATCCGGTTTCGGCGGCATCACGATTTTCAGTACCGCATACAGA  >  REL606/1079123‑1079268

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: