Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,586,533 G→A 10.1% G251G (GGC→GGT PP_4914 ← phytase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,586,5330GA10.1% 163.1 / 6.7 59G251G (GGC→GGTPP_4914phytase‑like protein
Reads supporting (aligned to +/- strand):  ref base G (25/28);  new base A (3/3);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

TGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGCCAC  >  NC_002947/5586389‑5586679
                                                                                                                                                |                                                                                                                                                  
tGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCt                                                                                                                                                >  8:244996/1‑149 (MQ=255)
  cTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGCAAGCTGTTTACCCTCg                                                                                                                                              <  1:17336/149‑1 (MQ=255)
   ttCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCg                                                                                                                                              >  5:179781/1‑148 (MQ=255)
   ttCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCg                                                                                                                                              >  1:245161/1‑148 (MQ=255)
     cGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGCAAGCTGTTTACCCTCGa                                                                                                                                             <  6:35985/147‑1 (MQ=255)
      ggccaggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAAGGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAc                                                                                                                                           <  3:302533/148‑1 (MQ=255)
       gccaggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAcca                                                                                                                                         <  7:244996/149‑1 (MQ=255)
        ccaggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccag                                                                                                                                        <  6:206741/149‑1 (MQ=255)
         caggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccagg                                                                                                                                       <  3:130588/149‑1 (MQ=255)
            gccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggcca                                                                                                                                    <  6:179781/149‑1 (MQ=255)
             ccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggccag                                                                                                                                   <  7:185994/149‑1 (MQ=255)
                 aCGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGATAAACAGCTTGCCCTCGAACCAGGCCAGGt                                                                                                                                 >  1:13025/1‑147 (MQ=255)
                         gcgTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTc                                                                                                                       <  6:205493/149‑1 (MQ=255)
                           gTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTcgc                                                                                                                     <  8:20791/149‑1 (MQ=255)
                             aGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCAATGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGCAAGCTGTTTACCCTCGAAcgtaac                                                                                                                                      >  2:253145/1‑125 (MQ=38)
                             aGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCAATGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGCAAGCTGTTTACCCTCGAAcgtaac                                                                                                                                      <  1:253145/130‑6 (MQ=38)
                              gcgcCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTcgcgcg                                                                                                                  >  7:261795/1‑149 (MQ=255)
                                   gTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAgg                                                                                                             <  1:145696/149‑1 (MQ=255)
                                    tGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAggg                                                                                                            <  8:120736/149‑1 (MQ=255)
                                        tCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCt                                                                                                        >  1:7495/1‑149 (MQ=255)
                                        tCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCt                                                                                                        >  5:121137/1‑149 (MQ=255)
                                                cgcagccgcagcGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCTAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTgg                                                                                                <  7:124547/149‑1 (MQ=255)
                                                 gcATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTggg                                                                                               >  4:147855/1‑149 (MQ=255)
                                                 gcATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTggg                                                                                               <  5:182297/149‑1 (MQ=255)
                                                     cagcagCGAACGACCCGCACCGCTCGACCTTGCCACTGTCGATATCACGGCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGTACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATc                                                                                           <  3:245319/149‑1 (MQ=255)
                                                      agcagcGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCt                                                                                          <  3:147855/149‑1 (MQ=255)
                                                                gACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCa                                                                                                                               <  6:276573/102‑1 (MQ=255)
                                                                gACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCa                                                                                                                               >  5:276573/1‑102 (MQ=255)
                                                                gACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTg                                                                                <  2:85871/149‑1 (MQ=255)
                                                                gACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTg                                                                                >  1:249931/1‑149 (MQ=255)
                                                                 aCCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGc                                                                               <  8:261795/149‑1 (MQ=255)
                                                                     gCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGACGTCGCGCGACAGGGCCCTGGCATTGGGCAGCTGCGCGGGGTGCAttt                                                                           >  4:170229/1‑149 (MQ=255)
                                                                     gCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCAttt                                                                           <  6:121137/149‑1 (MQ=255)
                                                                                  ttGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTc                                                              >  5:54332/1‑149 (MQ=255)
                                                                                      cACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTc                                                          <  6:54332/149‑1 (MQ=255)
                                                                                           tCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagca                                                      <  8:189764/148‑1 (MQ=255)
                                                                                             gATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCagcagc                                                    <  5:132555/148‑1 (MQ=255)
                                                                                              atatCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCAc                                                  >  6:50831/1‑149 (MQ=255)
                                                                                                     gTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCa                                           <  4:301803/149‑1 (MQ=255)
                                                                                                     gTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCa                                           >  1:222879/1‑149 (MQ=255)
                                                                                                     gTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCa                                           >  3:148783/1‑149 (MQ=255)
                                                                                                             acCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCg                                     >  1:172613/1‑147 (MQ=255)
                                                                                                               ccGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccag                                                                                                                                        <  6:259839/46‑1 (MQ=255)
                                                                                                               ccGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccag                                                                                                                                        >  5:259839/1‑46 (MQ=255)
                                                                                                                 ggTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAgcagca                                                   <  4:280662/129‑1 (MQ=255)
                                                                                                                 ggTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAgcagca                                                   >  3:280662/1‑129 (MQ=255)
                                                                                                                 ggTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTc                               <  2:172613/149‑1 (MQ=255)
                                                                                                                         ttACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGa                       <  2:7495/149‑1 (MQ=255)
                                                                                                                             gTTCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGActgc                   >  7:165749/1‑149 (MQ=255)
                                                                                                                                cGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGActgctg                 <  4:148783/148‑1 (MQ=255)
                                                                                                                                cGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACACCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTgg                >  7:3990/1‑149 (MQ=255)
                                                                                                                                 gAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTgg                <  2:249931/148‑1 (MQ=255)
                                                                                                                                 gAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGgc               >  2:19057/1‑149 (MQ=255)
                                                                                                                                     gTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGCCTGCTGGCGCTc           >  3:232020/1‑149 (MQ=255)
                                                                                                                                      tAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGCCCGACTGCTGGCGCTc           >  5:296867/1‑148 (MQ=255)
                                                                                                                                             cTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCt                                                           >  6:1147/1‑93 (MQ=255)
                                                                                                                                             cTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGcca   >  8:83114/1‑149 (MQ=255)
                                                                                                                                             cTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGcca   >  6:290912/1‑149 (MQ=255)
                                                                                                                                              ttGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGccac  >  6:196185/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGCCAC  >  NC_002947/5586389‑5586679

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: