Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,638,708 G→C 7.4% G235A (GGC→GCC)  PP_3205 → fumarylacetoacetate hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,638,7080GC 7.4% 170.7 / 3.2 54G235A (GGC→GCC) PP_3205fumarylacetoacetate hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (27/23);  new base C (2/2);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

GGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCG  >  NC_002947/3638563‑3638854
                                                                                                                                                 |                                                                                                                                                  
ggcttTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGGACTGGAACGGCGAGTGGTTCGGCCATCCCCATGgcggcg                                                                                                                                                  <  8:246615/148‑1 (MQ=255)
 gcttTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCgccgcc                                                                                                                                                <  7:12766/149‑1 (MQ=255)
      ctcCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCa                                                                                                                                           >  8:184575/1‑149 (MQ=255)
       tcCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCAc                                                                                                                                          >  7:93542/1‑149 (MQ=255)
       tcCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGCCGCCATgg                                                                                                                                            >  2:390348/1‑146 (MQ=255)
       tcCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGCCGCCATgg                                                                                                                                            <  1:390348/147‑2 (MQ=255)
             ccAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGt                                                                                                                                    >  1:350812/1‑149 (MQ=255)
                  cccGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTcca                                                                                                                               <  7:184575/149‑1 (MQ=255)
                    cGCCAGCAGTTTCGCGCCCGTTGCGGTCACACCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTccac                                                                                                                              >  8:254781/1‑148 (MQ=255)
                        agcagTTTCGCGGCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGtt                                                                                                                         <  8:138051/149‑1 (MQ=255)
                             tttCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTg                                                                                                                    <  3:326442/149‑1 (MQ=255)
                             tttCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTg                                                                                                                    >  1:249990/1‑149 (MQ=255)
                                  cgcCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCAt                                                                                                               >  7:348871/1‑149 (MQ=255)
                                              tCACCCCGGATGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgacg                                                                                                   <  8:243754/149‑1 (MQ=255)
                                                  cccGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGGACTTCGGTTTCCACCAGTTGATTGCCCATGCTGCGCTgacgcgac                                                                                               <  2:350812/149‑1 (MQ=255)
                                                  cccGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgacgcgac                                                                                               >  6:74529/1‑149 (MQ=255)
                                                     ggACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCt                                                                                            >  6:264396/1‑149 (MQ=255)
                                                     ggACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCAGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCt                                                                                            >  1:235356/1‑149 (MQ=255)
                                                        cGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGa                                                                                         <  2:249990/149‑1 (MQ=255)
                                                        cGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGa                                                                                         <  4:274244/149‑1 (MQ=255)
                                                           aCTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCCCCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTg                                                                                      >  8:95683/1‑149 (MQ=255)
                                                               ggCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCAGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCgg                                                                                   <  2:235356/148‑1 (MQ=255)
                                                                gCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGc                                                                                 <  8:93542/149‑1 (MQ=255)
                                                                       cTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTAc                                                                          >  2:274998/1‑149 (MQ=255)
                                                                       cTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTAc                                                                          <  7:254781/149‑1 (MQ=255)
                                                                            gcgcCGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGc                                                                     <  1:349816/149‑1 (MQ=255)
                                                                                  ttccagTGCACGTGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGc                                                               >  2:345635/6‑149 (MQ=255)
                                                                                   ccGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCa                                                              >  1:300937/1‑149 (MQ=255)
                                                                                    cGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCAc                                                             >  6:364102/1‑149 (MQ=255)
                                                                                     gggTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCAcc                                                            <  5:364102/149‑1 (MQ=255)
                                                                                      agTGCACGTGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCg                                                           >  5:320329/2‑149 (MQ=255)
                                                                                          cACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTc                                                       <  3:193495/149‑1 (MQ=255)
                                                                                           aCTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGATTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTc                                                       >  6:191710/1‑148 (MQ=255)
                                                                                                ccGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACACGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGc                                                 >  8:128417/1‑149 (MQ=255)
                                                                                                 cGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCTACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGcc                                                <  2:374751/149‑1 (MQ=255)
                                                                                                   cTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACACGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGa                                              <  7:128417/149‑1 (MQ=255)
                                                                                                      cACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAgc                                            <  1:274998/148‑1 (MQ=255)
                                                                                                              cTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTgg                                   >  3:125370/1‑149 (MQ=255)
                                                                                                                 gAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCt                                >  1:118803/1‑149 (MQ=255)
                                                                                                                 gAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCt                                >  3:367204/1‑149 (MQ=255)
                                                                                                                 gAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCt                                <  5:74529/149‑1 (MQ=255)
                                                                                                                    cGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTCATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCgg                             >  5:307829/1‑149 (MQ=255)
                                                                                                                       cGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGt                          >  7:269995/1‑149 (MQ=255)
                                                                                                                       cGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGt                          <  8:348871/149‑1 (MQ=255)
                                                                                                                             gTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCg                    >  4:269502/1‑149 (MQ=255)
                                                                                                                              ttCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGc                   <  5:264396/149‑1 (MQ=255)
                                                                                                                              ttCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGc                   <  3:359519/149‑1 (MQ=255)
                                                                                                                                 ggCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGcc                  >  8:260732/1‑147 (MQ=255)
                                                                                                                                 ggCCATCCCCATGGCGGAGCCATGCACTTCGGTTGCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGa                <  6:293742/149‑1 (MQ=255)
                                                                                                                                  gCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAg               >  3:83696/1‑149 (MQ=255)
                                                                                                                                    cATCCCAATGGCGGCGCCATGCACTTCGGTTTCCACAAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGTGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGCGCGCAGCGTGGGCTCGGCGTGCATCGCcgcgcg             >  5:225741/1‑149 (MQ=255)
                                                                                                                                        tgcATGGCGCCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCgg                                                                >  2:97978/3‑94 (MQ=255)
                                                                                                                                        tgcATGGCGCCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCgg                                                                <  1:97978/92‑1 (MQ=255)
                                                                                                                                         ccATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGCGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAg               <  3:233807/142‑1 (MQ=255)
                                                                                                                                         ccATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGCGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAg               >  4:233807/1‑142 (MQ=255)
                                                                                                                                          cATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAgcgcacgct       <  7:113187/149‑6 (MQ=255)
                                                                                                                                                ggcgCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCg  >  5:122977/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCG  >  NC_002947/3638563‑3638854

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