Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 475,296 Δ1 bp coding (10/1677 nt) ybaL ← predicted transporter with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606475,2960C.98.9% 859.4 / 2.0 179coding (10/1677 nt)ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base C (2/0);  major base . (62/113);  minor base G (2/0);  total (66/113)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 186 total aligned reads displayed.

GTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGT  >  REL606/475224‑475366
                                                                        |                                                                      
  ttATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATg                                                                  <  2:6350816/76‑1 (MQ=255)
    aTTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGGGGCGGGATGc                                                                 >  1:3547869/1‑76 (MQ=255)
      tGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCAtt                                                              <  2:5640523/76‑1 (MQ=255)
      tGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGGGGGTGATGCAtt                                                              >  1:3896500/1‑76 (MQ=255)
       ggCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTc                                                             <  2:7531663/76‑1 (MQ=255)
       ggCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTc                                                             <  2:4544686/76‑1 (MQ=255)
        gCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTcc                                                            >  1:7355179/1‑76 (MQ=255)
        gCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTcc                                                            >  1:6354378/1‑76 (MQ=255)
         ccAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCg                                                           <  2:8146193/76‑1 (MQ=255)
         ccAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCg                                                           <  1:665421/76‑1 (MQ=255)
         ccAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCg                                                           <  1:854978/76‑1 (MQ=255)
         ccAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCg                                                           <  2:1509483/76‑1 (MQ=255)
          cAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGt                                                          <  1:4997444/76‑1 (MQ=255)
           aGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGtc                                                         <  2:3214997/76‑1 (MQ=255)
            gCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGGGATGCATTCCGtcc                                                        >  1:7632098/1‑75 (MQ=255)
             cATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGtctc                                                       <  2:1442569/76‑1 (MQ=255)
              aTGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTcc                                                      >  2:2503016/1‑76 (MQ=255)
               tGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCt                                                     <  2:4607278/76‑1 (MQ=255)
                gCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCtt                                                    <  1:589226/76‑1 (MQ=255)
                gCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCtt                                                    >  1:4557688/1‑76 (MQ=255)
                gCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCtt                                                    <  2:4870091/76‑1 (MQ=255)
                gCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCtt                                                    >  2:290264/1‑76 (MQ=255)
                gCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCtt                                                    <  2:2818549/76‑1 (MQ=255)
                 ccGAGGATAAAGGCGAGCAGAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCttt                                                   <  2:3783505/76‑1 (MQ=255)
                 ccGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCttt                                                   >  2:1095964/1‑76 (MQ=255)
                  cGAGGATAAAGGCGAGCACAAGGCTGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTGCtttt                                                  <  2:5547845/76‑1 (MQ=255)
                  cGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCtttt                                                  <  1:3810297/76‑1 (MQ=255)
                  cGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCtttt                                                  >  1:1558708/1‑76 (MQ=255)
                   gAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTc                                                 <  1:5238634/76‑1 (MQ=255)
                   gAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTc                                                 >  2:5087249/1‑76 (MQ=255)
                   gAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTc                                                 >  2:5365867/1‑76 (MQ=255)
                   gAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTc                                                 >  1:6266667/1‑76 (MQ=255)
                     ggATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCt                                               <  1:2547997/76‑1 (MQ=255)
                       aTAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCtgg                                             <  1:4318029/76‑1 (MQ=255)
                        tAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCtggt                                            >  1:4760062/1‑76 (MQ=255)
                        tAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCtggt                                            >  1:53684/1‑76 (MQ=255)
                        tAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCtggt                                            <  2:7147515/76‑1 (MQ=255)
                         aaaGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCtggtg                                           >  2:5394956/1‑76 (MQ=255)
                         aaaGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCtggtg                                           >  2:7269907/1‑76 (MQ=255)
                          aaGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCtggtgg                                          <  2:1814281/76‑1 (MQ=255)
                          aaGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCtggtgg                                          >  1:1221234/1‑76 (MQ=255)
                           aGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTggtggt                                         <  2:5834887/76‑1 (MQ=255)
                           aGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTggtggt                                         <  2:6224457/76‑1 (MQ=255)
                            ggCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGtt                                        <  1:6021835/76‑1 (MQ=255)
                             gCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTa                                       >  2:5934094/1‑76 (MQ=255)
                             gCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTa                                       <  2:6271499/76‑1 (MQ=255)
                              cGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTAt                                      <  1:1917886/76‑1 (MQ=255)
                               gAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTAtt                                     <  2:8077952/76‑1 (MQ=255)
                                aGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTg                                    <  1:107574/76‑1 (MQ=255)
                                  cacaAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTc                                  <  2:3187922/76‑1 (MQ=255)
                                  cacaAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTc                                  >  1:6048566/1‑76 (MQ=255)
                                    caAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCa                                <  2:3252372/76‑1 (MQ=255)
                                    caAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCa                                <  1:5643496/76‑1 (MQ=255)
                                    caAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCa                                <  1:8025658/76‑1 (MQ=255)
                                    caAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCAATCCGTCTCCTTTTCCTGGTGGTTATTGTCCa                                <  2:6466971/76‑1 (MQ=255)
                                     aaGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTccct                               >  2:539515/1‑74 (MQ=255)
                                        gccgccAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAtttt                            >  2:5138618/1‑76 (MQ=255)
                                        gccgccAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAtttt                            >  1:2670900/1‑76 (MQ=255)
                                        gccgccAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAtttt                            <  1:3088429/76‑1 (MQ=255)
                                         ccgccAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAttttt                           <  2:5031916/76‑1 (MQ=255)
                                          cgccAACAATGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTg                          <  1:8125214/76‑1 (MQ=255)
                                             caacaaTGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGcc                       <  1:3749121/76‑1 (MQ=255)
                                              aacaaTGGGGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCg                      >  2:4166755/1‑76 (MQ=255)
                                               acaaTGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCgg                     >  1:7535328/1‑76 (MQ=255)
                                               acaaTGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCgg                     >  1:5177431/1‑76 (MQ=255)
                                                 aaTGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGa                   <  2:3599934/76‑1 (MQ=255)
                                                 aaTGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGa                   >  2:7699197/1‑76 (MQ=255)
                                                 aaTGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGa                   <  2:2963584/76‑1 (MQ=255)
                                                 aaTGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGa                   >  2:3596207/1‑76 (MQ=255)
                                                  aTGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaa                  >  1:5272012/1‑76 (MQ=255)
                                                  aTGGTGGTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaa                  >  1:4605191/1‑76 (MQ=255)
                                                   tggtggTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaaa                 >  2:4245857/1‑76 (MQ=255)
                                                   tggtggTGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaaa                 <  1:4649345/76‑1 (MQ=255)
                                                    tgtggtgATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaaaa                <  2:3515886/75‑1 (MQ=255)
                                                    ggtggtGATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaaaa                <  1:5028995/76‑1 (MQ=255)
                                                     gtggtgATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaaaa                >  2:6055307/1‑76 (MQ=255)
                                                     gtggtgATAAGCGGGGTGGGGGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAAc               >  1:5263060/1‑76 (MQ=255)
                                                      tggtgATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAAcc              <  1:5488595/76‑1 (MQ=255)
                                                       ggtgATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCa             <  2:3428380/76‑1 (MQ=255)
                                                       ggtgATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCa             <  2:4603377/76‑1 (MQ=255)
                                                       ggtgATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCa             <  1:7208390/76‑1 (MQ=255)
                                                        gtgATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaa            <  1:4338537/76‑1 (MQ=255)
                                                        gtgATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaa            <  1:4832295/76‑1 (MQ=255)
                                                        gtgATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaa            >  2:7910953/1‑76 (MQ=255)
                                                        gtgATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaa            <  1:5320269/76‑1 (MQ=255)
                                                         tgATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaaa           >  2:7573207/1‑76 (MQ=255)
                                                          gATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaaaa          <  2:3670651/76‑1 (MQ=255)
                                                          gATAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaaaa          <  2:2328147/76‑1 (MQ=255)
                                                           aTAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAAt         <  2:7862942/76‑1 (MQ=255)
                                                           aTAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAAt         <  2:2138583/76‑1 (MQ=255)
                                                            tAAGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAAtt        <  1:2985490/76‑1 (MQ=255)
                                                              aGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTAc      <  1:4952994/76‑1 (MQ=255)
                                                              aGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTAc      <  1:2754416/76‑1 (MQ=255)
                                                              aGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTAc      <  2:6784658/76‑1 (MQ=255)
                                                              aGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTAc      <  2:182805/76‑1 (MQ=255)
                                                              aGCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTAc      <  1:580466/76‑1 (MQ=255)
                                                               gCGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACa     <  1:2959897/76‑1 (MQ=255)
                                                                cGGGGTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAg    <  2:247656/76‑1 (MQ=255)
                                                                 ggggTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAgg   <  1:5067223/76‑1 (MQ=255)
                                                                 ggggTGGGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAgg   >  2:844617/1‑76 (MQ=255)
                                                                        |                                                                      
GTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGT  >  REL606/475224‑475366

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