Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,026,252 2 bp→GC coding (110‑111/1164 nt) fadA ← thiolase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,026,2520CG99.1% 294.7 / ‑4.6 113T37T (ACG→ACCfadAthiolase I
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (58/54);  total (58/55)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
*NC_0009134,026,2530GC99.1% 294.8 / ‑3.6 115coding (110/1164 nt)fadAthiolase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (59/55);  minor base . (0/1);  total (59/56)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 115 total aligned reads displayed.

CTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGCGTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCC  >  NC_000913/4026183‑4026320
                                                                      ||                                                                   
ctgctgCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGc                                                                >  1:3971911/1‑76 (MQ=255)
   ctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCAtt                                                             <  1:4154540/76‑1 (MQ=255)
    tgcACACACCCCCAGTAAATATCGGCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTa                                                            <  2:1384036/76‑1 (MQ=255)
    tgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCCGCAGGCTACGCCTTa                                                            >  2:3716955/1‑76 (MQ=255)
    tgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGGTACGCGCCAGCCGGCTACGCATTa                                                            >  1:619341/1‑76 (MQ=255)
     gcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGCCTACGCAGTaa                                                           >  2:2668995/1‑76 (MQ=255)
      cacacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTaaa                                                          >  1:4422792/1‑76 (MQ=255)
      cacacaACCCCAGTAAATATCGGCGGGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTaaa                                                          <  2:2265813/76‑1 (MQ=255)
       ccacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAAt                                                         <  2:3629749/75‑1 (MQ=255)
       acacaACCCCAGTAAATATCGTCGAGGGCCCCCGCTTCCCGCGCCGGGTTCCGCGCCCGCAGGCTCCGCATTAAAt                                                         >  2:4129988/1‑76 (MQ=255)
        cacaCCCCCAGTAAATAGCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATg                                                        <  1:1168438/76‑1 (MQ=255)
        cacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATg                                                        <  2:2008693/76‑1 (MQ=255)
        aacaACCCCAGTAAGTAGGTGGGAGGGCCGGCGCTTCCAGCGACGGGTTACGCGCCAGCAGGCTACGCATTAAATg                                                        <  2:1491667/75‑1 (MQ=255)
         acaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGa                                                       <  2:3971911/76‑1 (MQ=255)
          cacccccAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAg                                                      <  2:2417930/73‑1 (MQ=255)
          caACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAg                                                      >  1:1355167/1‑76 (MQ=255)
            cccccAGAAAATATCGTCGAGGCCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCg                                                    <  1:2704799/75‑1 (MQ=255)
            aCCCCAGTAAATATCGTCGGGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCg                                                    <  2:2621629/76‑1 (MQ=255)
            aCCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCg                                                    >  1:4647366/1‑76 (MQ=255)
            aCCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCg                                                    >  1:4353709/1‑76 (MQ=255)
            aCCCCAGAAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCg                                                    <  1:1592444/76‑1 (MQ=255)
             ccccAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCgg                                                   <  2:2771831/76‑1 (MQ=255)
               ccAGAAAATACCGTCGGGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgag                                                 <  2:4632381/76‑1 (MQ=255)
                  gTAAATATCGTTGAGGGCCGCCGGTTGCAGCGCCGGGGTACGCGCCAGCAGGCTACGCATTAAATGAGCGtagaga                                              >  1:542163/1‑76 (MQ=255)
                  gTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgagaga                                              <  1:4113809/76‑1 (MQ=255)
                   tAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAt                                             >  1:948676/1‑76 (MQ=255)
                      atatCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtctt                                          >  2:2941688/1‑76 (MQ=255)
                      atatCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtctt                                          >  1:26082/1‑76 (MQ=255)
                       tatCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCACCAGGCTACGCATTAAACGAGCGGAGCGAgcttc                                         >  2:4474786/1‑76 (MQ=255)
                       aatCGCCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtcttc                                         <  1:1663571/75‑1 (MQ=255)
                           gtcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCCGGCTACGCATTAAATGAGCGGAGAGATCCTCTgca                                     >  1:4195123/1‑76 (MQ=255)
                           gtcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgca                                     <  1:3230413/76‑1 (MQ=255)
                           gtcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgca                                     >  2:4302689/1‑76 (MQ=255)
                             cgggggCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacg                                   <  1:581912/73‑1 (MQ=255)
                             cgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacg                                   >  2:2488054/1‑76 (MQ=255)
                             cgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacg                                   <  2:2038616/76‑1 (MQ=255)
                             cgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacg                                   <  2:749526/76‑1 (MQ=255)
                             cgAGGGCCGCCGCTTCCAGCGCCGGGGTACGCGCCAGCAGGCGACGCATTAAATGAGCGGAGAGATCTTCTgcacg                                   >  1:338138/1‑76 (MQ=255)
                              gAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgc                                  <  1:3657017/76‑1 (MQ=255)
                              gAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgc                                  >  1:795580/1‑76 (MQ=255)
                              gAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgc                                  >  1:2714772/1‑76 (MQ=255)
                              gAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgc                                  >  2:1260271/1‑76 (MQ=255)
                              gAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgc                                  >  1:1271397/1‑76 (MQ=255)
                              gAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgc                                  >  2:888395/1‑76 (MQ=255)
                                 ggCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGcacgcacg                               <  1:2941688/76‑1 (MQ=255)
                                   ccgccgCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGtt                             <  2:3202628/76‑1 (MQ=255)
                                    cgccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTa                            >  2:2803136/1‑76 (MQ=255)
                                     gccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTAc                           >  1:3861075/1‑76 (MQ=255)
                                     gccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTAc                           <  1:2710316/76‑1 (MQ=255)
                                     gccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTAc                           >  1:1344330/1‑76 (MQ=255)
                                     gccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTAc                           <  1:4639305/76‑1 (MQ=255)
                                      ccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACg                          <  2:2482947/76‑1 (MQ=255)
                                      ccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACg                          >  1:2502272/1‑76 (MQ=255)
                                       cgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTCCGa                         >  2:3767167/1‑76 (MQ=255)
                                       cgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGa                         >  1:3429483/1‑76 (MQ=255)
                                       cgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGCTACGa                         >  1:3098703/1‑76 (MQ=255)
                                        gcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaa                        >  1:1841578/1‑76 (MQ=255)
                                        gcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaa                        >  1:3936427/1‑76 (MQ=255)
                                        gcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaa                        <  2:1327279/76‑1 (MQ=255)
                                        gcTTCCAGCGACGGGTTACGCTCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaa                        <  1:1455438/76‑1 (MQ=255)
                                        gcTTACAGCGCCGGGGTACGCGCCACCAGGCTACGCATGAAATGAGCGGAGAGATCTTCTACACGCACCATACgta                        >  2:110568/1‑74 (MQ=255)
                                         cTTGCAGCTCCGGGTTATGCGCCAGCATGCTGCGCATTACCTGAGTGGCGAGCTCTTCTGCGCGCACGTTACGaaa                       <  1:3625756/76‑1 (MQ=255)
                                          ttCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaaaa                      <  1:1824113/76‑1 (MQ=255)
                                          gtCCAGCGCAGGGATATGCGCCAGGAGGCTACGATTTAAATGAGCGGAGAGATAGTCTGCACGCATGCTGCGaaaa                      >  1:419564/2‑76 (MQ=255)
                                           tCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaaaac                     >  2:45183/1‑75 (MQ=255)
                                           tCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAg                     >  1:517175/1‑76 (MQ=255)
                                           tCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAg                     <  1:2213707/76‑1 (MQ=255)
                                           tCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGACCGGAAAGATCTCCTCCACCCACGTGACGaaaac                     >  2:272961/1‑75 (MQ=255)
                                           tCCAGCGCCGGGTTACGCGCAAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAg                     <  1:2219162/76‑1 (MQ=255)
                                            ccAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGc                    <  1:4595987/76‑1 (MQ=255)
                                             cAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCa                   <  1:2203512/76‑1 (MQ=255)
                                              aGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAc                  <  1:2149514/76‑1 (MQ=255)
                                              aGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAc                  >  2:2787378/1‑76 (MQ=255)
                                                cgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCg                >  1:3829295/1‑76 (MQ=255)
                                                cgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCg                >  1:3187121/1‑76 (MQ=255)
                                                  ccGGGTTACGCGCCAGCAGGCTACGCATTAAAGGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGcc              <  2:3058065/76‑1 (MQ=255)
                                                   cGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGccc             >  2:490863/1‑76 (MQ=255)
                                                   cGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGccc             >  1:1312474/1‑76 (MQ=255)
                                                    gggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCt            >  2:4808559/1‑76 (MQ=255)
                                                    gggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCt            <  2:4096617/76‑1 (MQ=255)
                                                    gggTTACGCGCCAGCAGGCTACGCATGAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCt            >  1:586614/1‑76 (MQ=255)
                                                     ggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCtt           <  1:66937/76‑1 (MQ=255)
                                                     ggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCtt           <  1:4105667/76‑1 (MQ=255)
                                                      gTTACTCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTc          <  1:1230274/76‑1 (MQ=255)
                                                      gTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTc          <  1:323165/76‑1 (MQ=255)
                                                       ttACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCg         >  2:4522225/1‑76 (MQ=255)
                                                       ttACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCg         >  1:1620535/1‑76 (MQ=255)
                                                       ttACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTGACGAAAAGCACCGCCCGTCg         >  2:829035/1‑76 (MQ=255)
                                                        tACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGa        <  2:1714110/76‑1 (MQ=255)
                                                        tACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGa        >  1:2254599/1‑76 (MQ=255)
                                                        tACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGa        <  1:472931/76‑1 (MQ=255)
                                                         gcgcgCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGaa       <  2:908515/75‑1 (MQ=255)
                                                          cgcgCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAAc      <  1:3496454/76‑1 (MQ=255)
                                                           gcgcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACg     >  2:3967269/1‑76 (MQ=255)
                                                           gcgcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACg     >  2:1300513/1‑76 (MQ=255)
                                                           gcgcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACg     >  1:3361108/1‑76 (MQ=255)
                                                           gcgcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACg     >  1:2005326/1‑76 (MQ=255)
                                                            cgcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACgg    <  1:4668628/76‑1 (MQ=255)
                                                             gcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGc   >  1:1887877/1‑76 (MQ=255)
                                                              ccAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGcc  <  1:2152718/76‑1 (MQ=255)
                                                                      ||                                                                   
CTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGCGTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCC  >  NC_000913/4026183‑4026320

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: