Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,809,176 G→T 10.5% R76L (CGT→CTT)  eno → enolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,809,1760GT10.5% 170.8 / 9.5 57R76L (CGT→CTT) enoenolase
Reads supporting (aligned to +/- strand):  ref base G (23/28);  new base T (3/3);  total (26/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 = 9.02e-01

GACAACGGCATCATCGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCC  >  NC_002947/1809031‑1809319
                                                                                                                                                 |                                                                                                                                               
gACAACGGCATCATCGGCAGCGCTTGTGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGa                                                                                                                                              <  7:36032/149‑1 (MQ=255)
gACAACGGCATCATCGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGt                                                                                                                                                <  5:94324/147‑1 (MQ=255)
         atcatcGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACctgctgct                                                                                                                                     >  6:171535/1‑149 (MQ=255)
         atcatcGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACctgctgct                                                                                                                                     >  1:312817/1‑149 (MQ=255)
            atcGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGa                                                                                                                                              <  3:225253/137‑1 (MQ=255)
            atcGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGa                                                                                                                                              >  4:225253/1‑137 (MQ=255)
            atcGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCAGTGACCTGCTGCTgg                                                                                                                                   <  1:366836/148‑1 (MQ=255)
            atcGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGACTCGCGCGAAGCGGTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCATCATCAACGGCCCGATCCGTGACCTGCTGCTggg                                                                                                                                  <  1:169212/149‑1 (MQ=255)
               ggCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCaa                                                                                                                               <  4:170062/149‑1 (MQ=255)
               ggCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCaa                                                                                                                               <  5:186350/149‑1 (MQ=255)
               ggCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCaa                                                                                                                               <  5:223024/149‑1 (MQ=255)
               ggCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCaa                                                                                                                               <  5:215357/149‑1 (MQ=255)
               ggCAGCGCTTGCGCGCCGTCCGGTGCTTCTACAGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGCCAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGc                                                                                                                                 >  5:285281/1‑147 (MQ=255)
                gCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTCGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGCAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAg                                                                                                                              >  1:272574/1‑149 (MQ=255)
                 cAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCacg                                                                                                                              >  4:151667/1‑146 (MQ=255)
                     gcTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATcg                                                                                                                         <  3:151667/149‑2 (MQ=255)
                      cTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATcgg                                                                                                                        >  2:269606/1‑147 (MQ=255)
                       ttGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCtt                                                                                                                       >  6:77220/1‑149 (MQ=255)
                         gcgcgcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTc                                                                                                                      <  7:189553/148‑1 (MQ=255)
                         gcgcgcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGACAAGGATCCTTc                                                                                                                      <  2:312817/148‑1 (MQ=255)
                          cgcgcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATcgggcc                                                                                                                     <  1:269606/148‑6 (MQ=255)
                          cgcgcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTcc                                                                                                                     <  5:218870/148‑1 (MQ=255)
                            cgcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGAt                                                                                                                  <  5:77220/149‑1 (MQ=255)
                             gcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTc                                                                                                                      >  8:39391/1‑144 (MQ=255)
                             gcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTc                                                                                                                      <  7:39391/144‑1 (MQ=255)
                             gcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATc                                                                                                                 <  2:286624/149‑1 (MQ=255)
                                   cGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAgg                                                                                                           >  1:326410/1‑149 (MQ=255)
                                     gTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTTCGATCAGAAGGct                                                                                                         >  2:139477/1‑149 (MQ=255)
                                     gTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGct                                                                                                         >  7:6363/1‑149 (MQ=255)
                                         ttCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTg                                                                                                      >  4:264126/1‑148 (MQ=255)
                                           cTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGa                                                                                                    <  8:322205/148‑1 (MQ=255)
                                              ccGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACcgc                                                                                                >  3:206323/1‑149 (MQ=255)
                                                 gcgcgcgcgAAGCGCGGGAGCGGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAAAATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGcc                                                                                             <  6:285281/146‑1 (MQ=255)
                                                     gcgcgAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGa                                                                                         <  3:264126/149‑1 (MQ=255)
                                                           aGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAAc                                                                                   >  6:205996/1‑149 (MQ=255)
                                                             cgcTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTg                                                                                 <  5:270813/149‑1 (MQ=255)
                                                              gcTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTgg                                                                                <  2:326410/149‑1 (MQ=255)
                                                                tGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGAc                                                                              >  1:339626/1‑149 (MQ=255)
                                                                    gCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTa                                                                          >  7:87550/1‑149 (MQ=255)
                                                                          tGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATcgg                                                                                                                        >  8:335040/1‑95 (MQ=25)
                                                                          tGATGGCGACAAGAGCCGTTACCTGGGAAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATcgg                                                                                                                        <  7:335040/97‑3 (MQ=25)
                                                                                 gACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGc                                                             >  7:221317/1‑149 (MQ=255)
                                                                                           gTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGgc                                                   <  8:221317/149‑1 (MQ=255)
                                                                                               ccTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGcc                                                >  4:173655/1‑148 (MQ=255)
                                                                                                   ggCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGc                                           >  7:76523/1‑149 (MQ=255)
                                                                                                    gCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGcc                                          <  8:6363/149‑1 (MQ=255)
                                                                                                        gggCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATcc                                      >  5:41588/1‑149 (MQ=255)
                                                                                                           cGTGCTGAAGGTCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTgg                                   <  1:128275/149‑1 (MQ=255)
                                                                                                            gTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGc                                  <  5:205996/149‑1 (MQ=255)
                                                                                                                  aaGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCAATCCTGGCTGTTTc                            >  8:52794/1‑149 (MQ=255)
                                                                                                                           accAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGc                                           >  4:103595/2‑125 (MQ=255)
                                                                                                                           accAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGc                                           <  3:103595/124‑1 (MQ=255)
                                                                                                                                      cGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAgctg        <  5:171535/149‑1 (MQ=255)
                                                                                                                                       ggCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACAGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAgctgc       >  4:241032/1‑149 (MQ=255)
                                                                                                                                        gCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTTCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAgctgct      <  1:139477/149‑1 (MQ=255)
                                                                                                                                            gATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGccc  <  6:41588/149‑1 (MQ=255)
                                                                                                                                            gATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGcc   >  6:346227/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GACAACGGCATCATCGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCC  >  NC_002947/1809031‑1809319

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: