Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,711,461 T→C 7.0% E77G (GAG→GGG)  paaC ← ring 1,2‑phenylacetyl‑CoA epoxidase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,711,4610TC 7.0% 162.6 / 2.2 57E77G (GAG→GGG) paaCring 1,2‑phenylacetyl‑CoA epoxidase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base C (2/2);  total (30/27)
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.53e-01

GGCAATCCCCGCGATTCGCGCATCGCTGGACTCGACCAACCCCTGCAGCACGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGCTCTTCCTCGATGGCCGGGGCACGGCCA  >  NC_002947/3711313‑3711607
                                                                                                                                                    |                                                                                                                                                  
ggcAATCCCCGCGATTCGCGCATCGCTGGACTCGACCAACCCCTGCAGCACGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCt                                                                                                                                                    >  4:229942/1‑149 (MQ=255)
 gcAATCCCCGCGATTCGCGCATCGCTGGACTCGACCAACCCCTGCAGCACGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTc                                                                                                                                                   <  6:286547/149‑1 (MQ=255)
       cccGCGATTCGCGCATCGCTGGACTCGACCAACCCCTGCAGCACGGCGTAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCag                                                                                                                                              <  1:373126/148‑2 (MQ=255)
       cccGCGATTCGCGCATCGCTGGACTCGACCAACCCCTGCAGCACGGCGAAGAGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCAcg                                                                                                                                             >  2:238207/1‑149 (MQ=255)
           cgATTCGCGCATCGCTGGACTCGACCAACCCCTGCAGCACGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCAcggcgg                                                                                                                                         <  1:177696/149‑1 (MQ=255)
              ttCGCGCATCGCTGGACTCGACCAACCCCTGCAGCACGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAAc                                                                                                                                      <  3:229942/149‑1 (MQ=255)
              ttCGCGCATCGCTGGACTCGACCAACCCCTGCAGCACGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGCTGCGGTAGGCGCGCTCGTCACGGCGGAAc                                                                                                                                      <  1:238207/149‑1 (MQ=255)
                 gcgcATCGCTGGACTCGACCAACCCCTGCAGCACGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGc                                                                                                                                    >  4:93427/1‑148 (MQ=255)
                  cgcATCGCTGGACTCGACCAACCCCTGCAGCACGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCa                                                                                                                                  <  3:93427/149‑1 (MQ=255)
                            gACTCGACCAACCCCTGCAGCACGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGggcatcgg                                                                                                                         <  5:329444/148‑1 (MQ=255)
                              cTCGACCAACCCCTGCAGCACGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGggcatcggcat                                                                                                                      >  8:179888/1‑149 (MQ=255)
                               tCGACCAACCCCTGCAGCACGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATAACCGTTGGGTTGCTCGACCAGCAGCAGGTAGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGggcatcggcatc                                                                                                                     >  6:62746/1‑149 (MQ=255)
                                 gACCAACCGCTGCAGCACGGCGAAGTGCCAGGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGTTTGCTCGACCAGCAGCAGGTTGCTGTAGGCGCGCTCGGCACGGCGGAACGCCAGGGCATCGGCAtctc                                                                                                                   <  4:180296/149‑1 (MQ=255)
                                             cagcaCGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGa                                                                                                       <  7:179888/149‑1 (MQ=255)
                                                  cGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGagcagc                                                                                                  <  7:373887/149‑1 (MQ=255)
                                                            ccACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCa                                                                                         >  3:13732/1‑148 (MQ=255)
                                                            ccACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCCCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCggcggc                                                                                          >  4:364839/1‑147 (MQ=255)
                                                            ccACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCCCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCggcggc                                                                                          >  4:370979/1‑147 (MQ=255)
                                                            ccACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCCCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCggcggc                                                                                          <  3:370979/147‑1 (MQ=255)
                                                            ccACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCCCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCggcggc                                                                                          <  3:364839/147‑1 (MQ=255)
                                                                  gtcgtAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCAtc                                                                                                            >  1:10886/1‑123 (MQ=255)
                                                                  gtcgtAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCAtc                                                                                                            <  2:10886/123‑1 (MQ=255)
                                                                       agagGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCg                                                                                                                  >  2:140740/1‑112 (MQ=255)
                                                                         agGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTccagccag                                                                             <  4:70691/147‑1 (MQ=255)
                                                                         agGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTccagccag                                                                             >  3:70691/1‑147 (MQ=255)
                                                                                 cTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCa                                                                                         <  6:223063/127‑1 (MQ=255)
                                                                                 cTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCa                                                                                         >  5:223063/1‑127 (MQ=255)
                                                                                      tcatggtcaCGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGcc                                                               >  6:417928/1‑148 (MQ=255)
                                                                                          ggtcaCGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCggcggc                                                                                          >  6:144059/1‑117 (MQ=255)
                                                                                          ggtcaCGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCggcggc                                                                                          <  5:144059/117‑1 (MQ=255)
                                                                                                ggCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGggcatcggc                                                                                                                        >  4:290313/1‑81 (MQ=255)
                                                                                                ggCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGggcatcggc                                                                                                                        <  3:290313/81‑1 (MQ=255)
                                                                                                  cAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAgg                                                  <  5:417928/149‑1 (MQ=255)
                                                                                                   aaaaTCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGc                                                 <  3:32008/149‑1 (MQ=255)
                                                                                                       tCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGAc                                             >  6:157742/1‑149 (MQ=255)
                                                                                                             ttGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATGCTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAg                                                   >  6:121656/1‑137 (MQ=255)
                                                                                                             ttGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATGCTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAg                                                   <  5:121656/137‑1 (MQ=255)
                                                                                                                ggTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGTCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAg                                    <  5:62746/149‑1 (MQ=255)
                                                                                                                ggTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTtcat                                       >  4:93031/1‑146 (MQ=255)
                                                                                                                ggTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTtcat                                       <  3:93031/146‑1 (MQ=255)
                                                                                                                     cTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTAGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGc                                 >  5:344435/1‑147 (MQ=255)
                                                                                                                      tCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAg                              >  5:74586/1‑149 (MQ=255)
                                                                                                                          ccAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGctct                          >  4:391537/1‑149 (MQ=255)
                                                                                                                          ccAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGctct                          >  4:296768/1‑149 (MQ=255)
                                                                                                                                 cagGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGCTCTTCCTCGa                   >  3:376344/1‑149 (MQ=255)
                                                                                                                                 cagGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGA‑CAGATCGAGGCCGACGTTCATCAGGGCCAGCTCTTCCTCGAt                  <  4:13732/149‑1 (MQ=255)
                                                                                                                                    gTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGCTCTTCCTCGATgg                >  5:159256/1‑149 (MQ=255)
                                                                                                                                     ttGCGGAAGGCGCGCTCGTCAAGGCGGAACGCCAGGGCATAGGCATCTCGGACATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTAGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGCTCTTCATCGATGGc               >  1:9328/1‑149 (MQ=255)
                                                                                                                                        cGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGCTCTTCCTCGATGGcgaa            <  4:376344/149‑4 (MQ=255)
                                                                                                                                        cGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGCTCTTCCTCGATGGCCgg            <  5:157742/149‑1 (MQ=255)
                                                                                                                                          gTAGGCGCGCTCGTCACGGTGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGCTCTTCCTCGATGGCCgggg          >  8:358330/1‑149 (MQ=255)
                                                                                                                                          gTAGGCGAGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGCTCTTCCTCGATGGCCgggg          >  2:396795/1‑149 (MQ=255)
                                                                                                                                              gcgcgcTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGCTCTTCCTCGATGGCCGGGGCACg      >  1:38538/1‑149 (MQ=255)
                                                                                                                                               cgcgcTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCTAGGCCGACGTTCATCAGGGCCAGCTCTTCCTCGATGGCCGGGGCACgg     <  2:9328/149‑1 (MQ=255)
                                                                                                                                                  gcTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGCTCTTCCTCGATGGCCGTGCCCCGGCca  >  2:20641/1‑149 (MQ=255)
                                                                                                                                                  gcTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGCTCTTCCTCGATGGCCGGGGCACGGCca  >  8:166112/1‑149 (MQ=255)
                                                                                                                                                   cTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGc                                 >  2:355955/1‑117 (MQ=255)
                                                                                                                                                   cTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGc                                 <  1:355955/117‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGCAATCCCCGCGATTCGCGCATCGCTGGACTCGACCAACCCCTGCAGCACGGCGAAGTGCCACGCGTCGTAGAGGAACTGCTTGGTCATGGTCACGGCAAAATCACCGTTGGGTTGCTCGACCAGCAGCAGGTTGCGGTAGGCGCGCTCGTCACGGCGGAACGCCAGGGCATCGGCATCTCGGCCATCGTCGAGCAGCTCGGCGGCATACTCCAGCCAGTTGCGGGCCTGGCCGACCAGATCGAGGCCGACGTTCATCAGGGCCAGCTCTTCCTCGATGGCCGGGGCACGGCCA  >  NC_002947/3711313‑3711607

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