Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 468,461 C→T 10.9% A832A (GCG→GCA)  PP_0386 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947468,4610CT10.9% 144.1 / 7.3 56A832A (GCG→GCA) PP_0386sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base T (3/3);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

GCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTGGGT  >  NC_002947/468326‑468607
                                                                                                                                       |                                                                                                                                                  
gcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTa                                                                                                                                       >  1:228257/1‑149 (MQ=255)
gcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGACAAGGCTGTCGTTGATGTGCTTGTAGCGGTCGAGGTCCAGGAACAGCACTGCGGCCTGGCGGTTGGACACCTGCGCACAGGTGAGt                                                                                                                                        <  1:136825/148‑1 (MQ=255)
        gCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTg                                                                                                                                >  1:131794/1‑148 (MQ=255)
         ctctTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGTGCGCAGGTgtccaac                                                                                                                                     <  5:137659/142‑7 (MQ=255)
         ctctTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGTGCGCAGGTgtccaac                                                                                                                                     >  6:137659/1‑136 (MQ=255)
                   cgcgGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                    <  6:156901/149‑1 (MQ=255)
                    gcTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                    >  2:320265/1‑148 (MQ=255)
                        ggCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACa                                                                                                               <  2:131794/149‑1 (MQ=255)
                          cGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAg                                                                                                             <  2:228257/149‑1 (MQ=255)
                                ccTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGa                                                                                                       >  2:110989/1‑149 (MQ=255)
                                 cTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGAt                                                                                                      >  6:211778/1‑149 (MQ=255)
                                 cTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCACGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCAAACGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGAt                                                                                                      >  2:288958/1‑149 (MQ=255)
                                 cTTGAGCAGCAGGTCGCCGACAGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGAt                                                                                                      >  3:97164/1‑149 (MQ=255)
                                     agcagcAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTggg                                                                                                  <  5:194794/149‑1 (MQ=255)
                                      gcagcaGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACTGCACTGCGCCCTGGCGGTTGGACACCTGCGCACTGGTGAGTACGCCCTGGAGCCGGTTTTCTACCAAGGTGCGATTggg                                                                                                  <  1:313813/148‑1 (MQ=255)
                                       cagcagGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTTCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCa                                                                                                <  1:110989/149‑1 (MQ=255)
                                          cagGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGc                                                                                             >  4:80673/1‑149 (MQ=255)
                                           agGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGcc                                                                                            >  4:66165/1‑149 (MQ=255)
                                             gTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCgg                                                                                          >  3:56299/1‑149 (MQ=255)
                                                  cGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGc                                                                                     >  8:226725/1‑149 (MQ=255)
                                                      cGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACTTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGt                                                                                 >  1:7127/1‑149 (MQ=255)
                                                            gCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCAGTGCGCCCTGGCGGTTGGACACCTGGGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCgtggt                                                                           <  1:288958/149‑1 (MQ=255)
                                                             ccAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGt                                                                                                                                            <  8:164280/83‑1 (MQ=255)
                                                              cAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGtggtgg                                                                         >  5:205824/1‑149 (MQ=255)
                                                              cAAGGCTGTCGTTGATGTGCTTGAAGCGGTCAAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGACGtggtgg                                                                         >  1:305033/1‑149 (MQ=255)
                                                               aaGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCgtggt                                                                           >  7:169589/1‑146 (MQ=255)
                                                               aaGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCgtggt                                                                           <  8:169589/146‑1 (MQ=255)
                                                                 ggCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCt                                                                      >  1:22649/1‑149 (MQ=255)
                                                                      tCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATa                                                                 <  5:211778/149‑1 (MQ=255)
                                                                         ttGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGTGCGCAGGTgt                                                                                                                                          >  7:297100/1‑72 (MQ=255)
                                                                         ttGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGTGCGCAGGTgt                                                                                                                                          <  8:297100/73‑2 (MQ=255)
                                                                               tgCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGttt                                                        <  3:66165/149‑1 (MQ=255)
                                                                                 cTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTgg                                                      <  6:205824/149‑1 (MQ=255)
                                                                                    gAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGt                                                          <  4:292431/142‑1 (MQ=255)
                                                                                    gAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGt                                                          >  3:292431/1‑142 (MQ=255)
                                                                                      aGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTgg                                                 <  4:56299/149‑1 (MQ=255)
                                                                                         ggTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCg                                              >  1:42217/1‑149 (MQ=255)
                                                                                           tCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTg                                            >  1:258033/1‑149 (MQ=255)
                                                                                           tCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTg                                            <  1:320265/149‑1 (MQ=255)
                                                                                                gTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACTTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGa                                       <  2:7127/149‑1 (MQ=255)
                                                                                                     ggAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTg                                  <  2:305033/149‑1 (MQ=255)
                                                                                                       aaCAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGga                                <  2:22649/149‑1 (MQ=255)
                                                                                                          aGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGAt                             >  6:267639/1‑149 (MQ=255)
                                                                                                              cTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGTTTGGGCAGGCCGGTCAGCGGGGCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCt                         <  7:226725/149‑1 (MQ=255)
                                                                                                              cTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCt                         <  3:80673/149‑1 (MQ=255)
                                                                                                                 cgcCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGc                       >  2:326988/1‑148 (MQ=255)
                                                                                                                 cgcCCTGGCGGTTGGACACCTGCGCACAGGGGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTg                        <  4:141997/147‑1 (MQ=255)
                                                                                                                       ggCGGTCGGACACCTGCGCACCGGTGAGTACGCCCCGGAGCCGGTGTACGACCAACGTGCGATTGGGCAGGCCGGGCAGCGGGGCGTGGTGGGCTTGATACTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAAcc                <  4:97164/149‑1 (MQ=255)
                                                                                                                        gCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAAcc                <  1:326988/148‑1 (MQ=255)
                                                                                                                          ggTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGGCGGTCAGCGGGGCGTGGTGGGCTTGATATTCGAGTTTGGCCTTGGCGTGCTTGATGCTGGAGATGTGTGGAAAAACCGcc             <  6:280848/149‑1 (MQ=255)
                                                                                                                             tactcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGc                                               >  1:102110/5‑112 (MQ=255)
                                                                                                                             tactcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGc                                               <  2:102110/108‑1 (MQ=255)
                                                                                                                               gACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTCCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGaa        >  4:21802/1‑149 (MQ=255)
                                                                                                                                  acCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTg     >  5:100743/1‑149 (MQ=255)
                                                                                                                                   ccTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTTCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTgg    <  3:306628/149‑1 (MQ=255)
                                                                                                                                    cTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTggg   >  4:270772/1‑149 (MQ=255)
                                                                                                                                       cgcACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTGGGt  <  3:270772/147‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                  
GCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTGGGT  >  NC_002947/468326‑468607

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