Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,572,250 G→A 8.9% L124L (TTG→TTA) ‡ pcaB → 3‑carboxy‑cis,cis‑muconate cycloisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,572,2500GA 8.9% 156.8 / 4.8 56L124L (TTG→TTA) ‡pcaB3‑carboxy‑cis,cis‑muconate cycloisomerase
Reads supporting (aligned to +/- strand):  ref base G (27/24);  new base A (3/2);  total (30/26)
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.46e-01

GCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGTGC  >  NC_002947/1572102‑1572393
                                                                                                                                                    |                                                                                                                                               
gCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTg                                                                                                                                                 <  8:231222/149‑1 (MQ=255)
 cAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGa                                                                                                                                                <  8:55316/149‑1 (MQ=255)
  aaCTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGAt                                                                                                                                               >  4:108655/1‑149 (MQ=255)
       ggCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGa                                                                                                                                            <  3:41046/147‑1 (MQ=255)
          cATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGcc                                                                                                                                        >  8:147836/1‑148 (MQ=255)
            ttCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGAc                                                                                                                                     >  7:202785/1‑149 (MQ=255)
               cGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGTACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTc                                                                                                                                  >  2:162491/1‑149 (MQ=255)
               cGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTc                                                                                                                                  >  4:65254/1‑149 (MQ=255)
                    gTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGc                                                                                                                               >  1:100707/1‑147 (MQ=255)
                       aaaGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCt                                                                                                                          >  8:274670/1‑149 (MQ=255)
                       aaaGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATCAAATCGAgg                                                                                                                                          <  3:127322/133‑1 (MQ=255)
                       aaaGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATCAAATCGAgg                                                                                                                                          >  4:127322/1‑133 (MQ=255)
                         aGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATCAAATCGAggg                                                                                                                                         >  2:89501/1‑131 (MQ=255)
                         aGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATCAAATCGAggg                                                                                                                                         <  1:89501/132‑2 (MQ=255)
                           cGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATTGAGGCCGACCTCGGCAAGCTGGcc                                                                                                                      <  3:226989/149‑1 (MQ=255)
                           cGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGcc                                                                                                                      >  8:33155/1‑149 (MQ=255)
                                ggCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCCAGCTTGCCGagg                                                                                                                  <  7:147836/148‑3 (MQ=255)
                                 gCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATAcc                                                                                                                <  7:33155/149‑1 (MQ=255)
                                  cAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATAccc                                                                                                               >  2:163188/1‑149 (MQ=255)
                                         aTCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCg                                                                                                          >  7:248637/1‑147 (MQ=255)
                                          tCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCgcag                                                                                                       <  7:274670/149‑1 (MQ=255)
                                              cACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAgg                                                                                                   >  3:177042/1‑149 (MQ=255)
                                                     gTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGaa                                                                                            <  5:217570/149‑1 (MQ=255)
                                                     gTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGaa                                                                                            <  5:251675/149‑1 (MQ=255)
                                                       gCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGc                                                                                          >  4:69664/1‑149 (MQ=255)
                                                                tGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATa                                                                                 >  3:141295/1‑149 (MQ=255)
                                                                  aGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATAc                                                                                <  8:128529/148‑1 (MQ=255)
                                                                           tgtaCCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTtggt                                                                       <  1:162491/145‑1 (MQ=255)
                                                                                cTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCg                                                                 >  8:2603/1‑149 (MQ=255)
                                                                                cTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCg                                                                 >  4:178601/1‑149 (MQ=255)
                                                                                       ccaccaGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTgg                                                          >  6:60883/1‑149 (MQ=255)
                                                                                       acaccacCCCGCATGAGATCGATACCGGTCGGCGACTGCAGCTGCGCGATGCCCACGATTTGAGCGAGGCCGACCTCGGCAAGCTCGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTgg                                                          <  8:288532/148‑1 (MQ=255)
                                                                                                       gATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCAccc                                          <  8:248637/149‑1 (MQ=255)
                                                                                                       gATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCAccc                                          >  2:269028/1‑149 (MQ=255)
                                                                                                        aTGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTc                                                                  <  8:289467/124‑1 (MQ=255)
                                                                                                        aTGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTc                                                                  >  7:289467/1‑124 (MQ=255)
                                                                                                        aTGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCAcccc                                         <  2:100707/149‑1 (MQ=255)
                                                                                                             tACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGa                                    >  6:281306/1‑149 (MQ=255)
                                                                                                                cGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGAccc                                 >  5:21549/1‑149 (MQ=255)
                                                                                                                       gTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCAt                          >  5:40097/1‑149 (MQ=255)
                                                                                                                        ttCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCAt                          >  1:112327/1‑148 (MQ=255)
                                                                                                                           tGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGaaa                      >  3:21302/1‑149 (MQ=255)
                                                                                                                             cAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTggg                             >  4:291805/1‑140 (MQ=255)
                                                                                                                             cAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTggg                             <  3:291805/140‑1 (MQ=255)
                                                                                                                             cAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGaaa                      <  3:65254/147‑1 (MQ=255)
                                                                                                                              aGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGATGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTg                   <  1:163188/149‑1 (MQ=255)
                                                                                                                               gCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTgg                  <  3:178601/149‑1 (MQ=255)
                                                                                                                                   cgcgATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCg               <  3:108655/148‑1 (MQ=255)
                                                                                                                                    cgaatctccTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCgg              <  1:110131/141‑1 (MQ=255)
                                                                                                                                          gccTCGATCAAATCGAGGCCGCCCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCATGGGCATGAAACTGGCCGGTGTACTg       >  3:248724/2‑149 (MQ=255)
                                                                                                                                          cccTCGAGTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGGCGTACCTGGTTGCACCACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTg       <  4:3883/149‑1 (MQ=255)
                                                                                                                                           ccTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTg                                                       >  2:275238/1‑100 (MQ=255)
                                                                                                                                           ccTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTgg      <  4:141295/149‑1 (MQ=255)
                                                                                                                                            cTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACAGGCCgg              >  2:110131/1‑140 (MQ=255)
                                                                                                                                               gATTTGATCGAGTCCGACCTCGGCAAGCTGGCCGATACACTGTCGCAGCAGGCCTTGAAGCACGCAGATACGCACTTGGTGGGTCGTACCTGGTTGCAACACGCCCCCCCGGTGACCCTGGGCAGGAAACTGGCCGGTGGACTggggcc  >  8:220869/1‑146 (MQ=255)
                                                                                                                                               gATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGt    <  5:60883/147‑1 (MQ=255)
                                                                                                                                               gATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCCGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGTg   <  5:281306/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGTGC  >  NC_002947/1572102‑1572393

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑