Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,192,526 T→G 10.4% Q312P (CAG→CCG)  prr ← gamma‑aminobutyraldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,192,5260TG10.4% 157.8 / 10.1 58Q312P (CAG→CCG) prrgamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base G (3/3);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

AAGAACCCGGCACCTGCGGCGCGCTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGTTGTAAAAGCCGAAGGCGCGAATGCCCTCGACGGCGGC  >  NC_002947/3192385‑3192674
                                                                                                                                             |                                                                                                                                                     
aagaaCCCGGCACCTGCGGCGCGCTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCAGGATGACCCGGCAAcc                                                                                                                                                 >  3:302234/1‑148 (MQ=255)
   aaCCCGGCACCTGCGGCGCGCTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCAGGATGACCCGGCAAccgaa                                                                                                                                              >  4:200244/1‑145 (MQ=255)
   aaCCCGGCACCTGCGGCGCGCTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCAGGATGACCCGGCAAccgaa                                                                                                                                              <  3:200244/148‑4 (MQ=255)
    aCCCGGCACCTGCGGCGCGCTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACAt                                                                                                                                            >  6:191561/1‑149 (MQ=255)
      ccGGCACCTGCGGCGCGCTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACAttt                                                                                                                                          <  3:53893/149‑1 (MQ=255)
      ccGGCACCTGCGGCGCGCTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCAGGATGACCCGGCAAccgaaatgg                                                                                                                                          <  4:302234/149‑8 (MQ=255)
             cTGCGGCGCGCTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGc                                                                                                                                    <  2:173228/147‑1 (MQ=255)
             cTGCGGCGCGCTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGc                                                                                                                                    >  1:173228/1‑147 (MQ=255)
               gcggcgCGCTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCAGGATGACCCGGCAAc                                                                                                                                                  >  6:22586/1‑132 (MQ=255)
               gcggcgCGCTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCAGGATGACCCGGCAAc                                                                                                                                                  <  5:22586/132‑1 (MQ=255)
                   cgcgcTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACg                                                                                                                            >  6:110700/1‑149 (MQ=255)
                    gcgcTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACgg                                                                                                                           <  5:191561/149‑1 (MQ=255)
                          gccgccCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCTGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCg                                                                                                                      <  5:110700/148‑1 (MQ=255)
                                 gTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCtt                                                                                                              <  5:276647/149‑1 (MQ=255)
                                       acCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGcgaac                                                                                                        <  1:161766/149‑1 (MQ=255)
                                             atatGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATTCTGGACACGGCTTCGCCAAGCATGCGAACGAACTc                                                                                                  >  2:253665/1‑149 (MQ=255)
                                                       gctCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATg                                                                                                                                     <  3:14017/104‑1 (MQ=255)
                                                       gctCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATg                                                                                                                                     >  4:14017/1‑104 (MQ=255)
                                                       gctCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGAcc                                                                                        <  7:188982/149‑1 (MQ=255)
                                                        ctCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTg                                                                                                                                                      <  2:255674/87‑1 (MQ=255)
                                                        ctCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTg                                                                                                                                                      >  1:255674/1‑87 (MQ=255)
                                                        ctCCCTGGCCCGGTAGACGAACCCCTCGACACGCTGCAGGTGCTCTTCGGTGATCACCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGAc                                                                                         <  6:194695/147‑1 (MQ=255)
                                                          cccTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGCCACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGcc                                                                                     >  6:206555/1‑149 (MQ=255)
                                                           ccTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTc                                                                                                                       <  2:82324/114‑1 (MQ=255)
                                                           ccTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTc                                                                                                                       >  1:82324/1‑114 (MQ=255)
                                                           ccTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCt                                                                                    >  6:108344/1‑149 (MQ=255)
                                                               gCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCggg                                                                                              >  4:73918/1‑135 (MQ=255)
                                                               gCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCggg                                                                                              <  3:73918/135‑1 (MQ=255)
                                                                 ccGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGc                                                                                                            <  6:154645/119‑1 (MQ=255)
                                                                 ccGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGc                                                                                                            >  5:154645/1‑119 (MQ=255)
                                                                   ggTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGcaca                                                                            <  5:206555/149‑1 (MQ=255)
                                                                    gTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTTATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATAGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGcaca                                                                            <  1:253665/149‑1 (MQ=255)
                                                                          gAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGGTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGc                                                                     >  2:149483/1‑149 (MQ=255)
                                                                          gAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGc                                                                     >  2:40210/1‑149 (MQ=255)
                                                                                 tCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCGAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTg                                                                                <  2:182580/131‑1 (MQ=255)
                                                                                 tCGACACGCTGCAGGTGCTCTTGGGTGAACAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCGAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTg                                                                                >  1:182580/1‑131 (MQ=255)
                                                                                   gACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCgg                                                            >  6:130596/1‑149 (MQ=255)
                                                                                        gCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGGCATCCTGCAACCCACATTTTATGCTGGACACCGCTTCGCCAAGCTTGCGAACGAACACCTGCTAGACCGCCTTCTGCACATATAGGCGACAAGCGGCCGTa                                                       >  7:123696/1‑149 (MQ=255)
                                                                                                tGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTGATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGcc                                                          >  1:12528/1‑138 (MQ=255)
                                                                                                tGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTGTGCACATAGAGGCGACACGCGGcc                                                          <  2:12528/138‑1 (MQ=255)
                                                                                                 gCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTgg                                              >  5:331561/1‑149 (MQ=255)
                                                                                                    ctTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTc                                                  <  4:229466/142‑1 (MQ=255)
                                                                                                    ctTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTc                                                  >  3:229466/1‑142 (MQ=255)
                                                                                                     ttGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGTGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCgg                                          >  2:130878/1‑149 (MQ=255)
                                                                                                       gggTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGTGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGc                                         <  1:130878/148‑1 (MQ=255)
                                                                                                        ggTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGt                                       >  1:22845/1‑149 (MQ=255)
                                                                                                         gTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGtt                                      <  1:40210/149‑1 (MQ=255)
                                                                                                             tCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGt                                       <  7:13203/144‑1 (MQ=255)
                                                                                                             tCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGt                                       >  8:13203/1‑144 (MQ=255)
                                                                                                                    gCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGTTGTAAAAGCCg                            >  4:275577/1‑148 (MQ=255)
                                                                                                                    gCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGTTGTAAAAGCCg                            >  4:233295/1‑148 (MQ=255)
                                                                                                                     cccATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGc                                                                                                                >  2:303220/1‑63 (MQ=255)
                                                                                                                     cccATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGc                                                                                                                <  1:303220/63‑1 (MQ=255)
                                                                                                                      ccATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGTTGTAAAAGCCg                            >  4:180802/1‑146 (MQ=255)
                                                                                                                      ccATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGTTGTAAAAGCCGaa                          <  3:180802/148‑1 (MQ=255)
                                                                                                                       cATGTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGAACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGTTGTAAAAGCCGAAgg                        <  8:123696/149‑1 (MQ=255)
                                                                                                                             ggTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGTTGTAAAAGCCGAAGGCGCGa                   <  5:130596/148‑1 (MQ=255)
                                                                                                                              gTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGTTGTAAAAGCCGAAGGCGCGAAt                 >  7:299143/1‑149 (MQ=255)
                                                                                                                                             tGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGTTGTAAAAGCCGAAGGCGCGAATGCCCTCGACggcggc  >  1:50341/1‑149 (MQ=255)
                                                                                                                                             tGCAACACACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGTTGTAAAAGCCGAAGGCGCGAATGCCCTCGACggcggc  >  5:101778/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
AAGAACCCGGCACCTGCGGCGCGCTTGCCGCCCGTGACGACCTCTATATGGGGCTGCTCCCTGGCCCGGTTGACGAAGCCCTCGACACGCTGCAGGTGCTCTTGGGTGATCAGCGGGCCCATTTCGGTTGCCGGGTCATCCTGCAACCCACATTTGATGCTGGACACGGCTTCGCCAAGCTTGCGAACGAACTCCGGGTAGACCGCCTTCTGCACATAGAGGCGACACGCGGCCGTACAGTCCTGGCCGGCGTTGTAAAAGCCGAAGGCGCGAATGCCCTCGACGGCGGC  >  NC_002947/3192385‑3192674

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: