Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 819,703 T→A 8.1% A61A (GCT→GCA)  ada → fused O(6)‑alkylguanine DNA alkyltransferase/transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947819,7030TA 8.1% 172.8 / 3.1 59A61A (GCT→GCA) adafused O(6)‑alkylguanine DNA alkyltransferase/transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (30/24);  new base A (3/2);  total (33/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 = 1.00e+00

GAGAGCCGCGACACTGCCGCTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAGTGCACGGCTGGCCGTCAG  >  NC_002947/819563‑819846
                                                                                                                                            |                                                                                                                                               
gagagCCGCGACACTGCCGCTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTa                                                                                                                                         >  2:155159/1‑149 (MQ=255)
      cgcgACACTGCCGCTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGc                                                                                                                                   >  8:76585/1‑149 (MQ=255)
                 cGCTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCAGCCTCGGCAGCCgcggca                                                                                                                                        >  4:116547/1‑128 (MQ=255)
                 cGCTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCAGCCTCGGCAGCCgcggca                                                                                                                                        <  3:116547/133‑6 (MQ=255)
                   cTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGt                                                                                                                      <  3:194566/149‑1 (MQ=255)
                           cACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCaa                                                                                                              <  5:163852/149‑1 (MQ=255)
                             cTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACCACACTGCCGCGGCTGCCGAGGCAGCCGGCTACCGTGCTTGCCAACGCTGTATCAGCAAgg                                                                                                            >  5:172032/1‑149 (MQ=255)
                                tGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGc                                                                                                           <  3:243002/147‑1 (MQ=255)
                                                accaccGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCAcag                                                                                         >  1:248397/1‑149 (MQ=255)
                                                  caccGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCaa                                                                                                              <  1:66232/126‑1 (MQ=255)
                                                  caccGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCaa                                                                                                              >  2:66232/1‑126 (MQ=255)
                                                  caccGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCAcagc                                                                                        >  2:75865/1‑148 (MQ=255)
                                                    ccGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCAcagccag                                                                                     <  2:248397/149‑1 (MQ=255)
                                                       gCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTg                                                                                  <  8:82017/149‑1 (MQ=255)
                                                            tACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGGTAc                                                                             >  3:47409/1‑149 (MQ=255)
                                                              cTGCGACCCCGGCTGCAATTCACGGCAGGGCAAGCGTATGATCGTCGAGTTGTACGACACTGGCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACtg                                                                           <  4:47409/149‑2 (MQ=255)
                                                               tGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCa                                                                          >  2:12909/1‑149 (MQ=255)
                                                                   aCCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCCCCCGCCCCAGCCCGCTGGTTACTCACGcc                                                                      >  4:28374/1‑149 (MQ=255)
                                                                   aCCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGcc                                                                      >  4:171094/1‑149 (MQ=255)
                                                                      ccGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCAc                                                                         <  5:250975/143‑1 (MQ=255)
                                                                      ccGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCAc                                                                         >  6:250975/1‑143 (MQ=255)
                                                                         gCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGc                                                                <  3:171094/149‑1 (MQ=255)
                                                                         gCTGCAAATAACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGc                                                                >  8:207405/1‑149 (MQ=255)
                                                                             cAACTCACGCCAGGCCACGCGTATGAAGGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTg                                                             <  6:172032/148‑1 (MQ=255)
                                                                             cAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCa                                                                          <  7:49903/135‑1 (MQ=255)
                                                                             cAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCa                                                                          >  8:49903/1‑135 (MQ=255)
                                                                              aaaTCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGAt                                                           <  1:155159/149‑1 (MQ=255)
                                                                              aaaTCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGAt                                                           <  5:132179/149‑1 (MQ=255)
                                                                               aaTCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATc                                                          >  6:84390/1‑149 (MQ=255)
                                                                               aaTCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATc                                                          <  3:163949/149‑1 (MQ=255)
                                                                                 tCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCg                                                         >  5:174420/1‑148 (MQ=255)
                                                                                     gCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAg                                                                                                        >  3:75299/1‑97 (MQ=255)
                                                                                     gCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAg                                                                                                        <  4:75299/97‑1 (MQ=255)
                                                                                           ccAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGAc                                              >  2:267033/1‑149 (MQ=255)
                                                                                            cAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGAcc                                             >  3:85601/1‑149 (MQ=255)
                                                                                            cAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGAcc                                             <  5:84390/149‑1 (MQ=255)
                                                                                                 gTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGtt                                        <  4:52352/149‑1 (MQ=255)
                                                                                                         gTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCt                                >  1:172078/1‑149 (MQ=255)
                                                                                                          tCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAgcg                                                                                                      >  4:125417/1‑78 (MQ=255)
                                                                                                          tCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAgcg                                                                                                      <  3:125417/78‑1 (MQ=255)
                                                                                                          tCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGcc                                 <  2:898/147‑1 (MQ=255)
                                                                                                          tCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGcc                                 >  1:898/1‑147 (MQ=255)
                                                                                                          tCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTc                               >  5:135724/1‑149 (MQ=255)
                                                                                                          tCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCAGCCACAGCCCGCTGGTTACTAACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTGCCAAGCCTc                               >  8:240048/1‑149 (MQ=255)
                                                                                                           cGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGc                                  >  1:123019/1‑145 (MQ=255)
                                                                                                           cGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGc                                  <  2:123019/145‑1 (MQ=255)
                                                                                                            gAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCgg                             <  1:75865/149‑1 (MQ=255)
                                                                                                                 cTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCaa                         <  6:174420/148‑1 (MQ=255)
                                                                                                                 cTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCaa                         <  2:172078/148‑1 (MQ=255)
                                                                                                                 cTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCaa                         >  8:233930/1‑148 (MQ=255)
                                                                                                                     gACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAg                    >  6:31636/1‑149 (MQ=255)
                                                                                                                      acacTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAGt                   >  3:94051/1‑149 (MQ=255)
                                                                                                                             cgcgGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAGTGCACGGc            >  7:199198/1‑149 (MQ=255)
                                                                                                                              ccgGCAGCCTCGGCAGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGAt                                                           <  3:112129/100‑1 (MQ=255)
                                                                                                                              ccgGCAGCCTCGGCAGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGAt                                                           >  4:112129/2‑101 (MQ=255)
                                                                                                                                    gCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAGTGCACGGCTGGCCGt     >  6:63685/1‑149 (MQ=255)
                                                                                                                                    gCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAGTGCACGGCTGGCCGt     >  6:237511/1‑149 (MQ=255)
                                                                                                                                       gcggCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAGTGCACGGCTGGCCGTCAg  <  5:237511/147‑1 (MQ=255)
                                                                                                                                       gAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAGTGCACGGCTGGCCGTCAg  >  7:203249/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                               
GAGAGCCGCGACACTGCCGCTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAGTGCACGGCTGGCCGTCAG  >  NC_002947/819563‑819846

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: ‑