Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,917,650 G→C 7.8% L236L (CTG→CTC) ‡ ygiF → thiamine triphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,917,6500GC 7.8% 154.9 / 4.3 51L236L (CTG→CTC) ‡ygiFthiamine triphosphatase
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base C (2/2);  total (24/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.39e-01

ACTGATGCCGTGCGACATCAGCAAGGCTGAGCGCGGTTACCGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTGGGCCAGGCAGCACCGCGCGCCACCACCCGCGACCTGCGCG  >  NC_002947/5917506‑5917798
                                                                                                                                                |                                                                                                                                                      
aCTGATGCCGTGCGACATCAGCAAGGCTGAGCGCGGTTACCGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGCCGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTggg                                                                                                                                                      >  4:97258/1‑147 (MQ=255)
aCTGATGCCGTGCGACATCAGCAAGGCTGAGCGCGGTTACCGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGgca                                                                                                                                                    >  3:13721/1‑149 (MQ=255)
     tGCCGTGCGACATCAGCAAGGCTGAGCGCGGTTACCGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCagcagc                                                                                                                                               <  2:84244/149‑1 (MQ=255)
         gTGCGACATCAGCAAGGCTGAGCGCGGTTACCGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCTGCTGCCCAGCAGctgccagg                                                                                                                                           <  7:61470/149‑2 (MQ=255)
         gTGCGACATCAGCAAGGCTGAGCGCGGTTACCGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGcagccagc                                                                                                                                           >  8:211123/1‑149 (MQ=255)
                 tCAGCAAGGCTGAGCGCGGTTACCGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCAGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGcc                                                                                                                                   >  8:261104/1‑149 (MQ=255)
                          cTGAGCGCGGTTACCGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTAc                                                                                                                          <  1:47131/149‑1 (MQ=255)
                            gAGCGCGGTTACCGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCg                                                                                                                        <  3:127070/149‑1 (MQ=255)
                                   gTTACCGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTgg                                                                                                                                                       <  4:253090/111‑1 (MQ=255)
                                   gTTACCGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTgg                                                                                                                                                       >  3:253090/1‑111 (MQ=255)
                                        cGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCa                                                                                                            <  1:21676/149‑1 (MQ=255)
                                        cGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCa                                                                                                            >  8:279421/1‑149 (MQ=255)
                                         gCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCAc                                                                                                           >  1:245494/1‑149 (MQ=255)
                                                    ccGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCt                                                                                              >  4:9119/1‑149 (MQ=255)
                                                     cgtaCAGCTACGAGCTTGACCTGCCGCTCTCCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGGCAGCGCCTGGGCGATCAGTACCGCCACAACGGCCTCTGG‑‑CGCCTGCTc                                                                                             <  2:229832/146‑1 (MQ=255)
                                                     cGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTc                                                                                             <  7:261104/149‑1 (MQ=255)
                                                        aCAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCCGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAg                                                                                          <  1:6265/149‑1 (MQ=255)
                                                            cTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGGCTCGACTGCTCCAGGa                                                                                        >  1:114442/1‑149 (MQ=255)
                                                             tACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTg                                                                                     >  6:164411/1‑149 (MQ=255)
                                                             tACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTg                                                                                     >  6:106440/1‑149 (MQ=255)
                                                               cGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGTCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTg                                                                                     <  7:211123/147‑1 (MQ=255)
                                                               cGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTggg                                                                                   <  2:54002/149‑1 (MQ=255)
                                                                  gCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCAGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTc                                                                                 >  8:101118/1‑148 (MQ=255)
                                                                   cTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTc                                                                                 >  5:251872/1‑147 (MQ=255)
                                                                   cTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGa                                                                               <  8:74290/149‑1 (MQ=255)
                                                                      gACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAg                                                                              <  2:245494/147‑1 (MQ=255)
                                                                       aCCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGc                                                                           >  1:66990/1‑149 (MQ=255)
                                                                        ccTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGa                                                                               <  6:215172/144‑1 (MQ=255)
                                                                        ccTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGa                                                                               >  5:215172/1‑144 (MQ=255)
                                                                        ccTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGcc                                                                          <  5:106440/149‑1 (MQ=255)
                                                                            ccGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTg                                                                        <  7:303678/147‑1 (MQ=255)
                                                                            ccGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTg                                                                        >  8:303678/1‑147 (MQ=255)
                                                                               cacaCCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTgagcga                                                                   <  2:66990/149‑1 (MQ=255)
                                                                                acacCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTgagccag                                                                  >  8:217884/1‑149 (MQ=255)
                                                                                  acCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCt                                                                >  5:224958/1‑149 (MQ=255)
                                                                                  acCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCt                                                                >  7:13144/1‑149 (MQ=255)
                                                                                             cGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCa                                                     >  4:173810/1‑149 (MQ=255)
                                                                                              gCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGCGCGCAGCTCGCTCag                                                    <  3:13717/149‑2 (MQ=255)
                                                                                                cGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCg                                                  <  3:9119/149‑1 (MQ=255)
                                                                                                        cccTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTg                                          <  7:101118/149‑1 (MQ=255)
                                                                                                         ccTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTgg                                         >  1:16832/1‑149 (MQ=255)
                                                                                                          cTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTggg                                        <  7:280867/149‑1 (MQ=255)
                                                                                                          cTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTggg                                        <  5:164411/149‑1 (MQ=255)
                                                                                                          cTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTggg                                        >  8:280867/1‑149 (MQ=255)
                                                                                                                gacgCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTGGGCCAGgc                                  <  6:251872/149‑1 (MQ=255)
                                                                                                                  cgCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCGCTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTGGGCCAGgcag                                <  4:13721/149‑1 (MQ=255)
                                                                                                                  cgCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTGGGCCAGgcag                                <  6:224958/149‑1 (MQ=255)
                                                                                                                                    cTGGCTGCTGCCCAGCAGCTGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTGGGCCAGGCAGCACCGCGCGCCACCAccc              >  5:302000/1‑149 (MQ=255)
                                                                                                                                    cTGGCTGCTGCCCAGCAGCTGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTGGGCCAGGCAGCACCGCGCGCCACCAccc              <  6:302000/149‑1 (MQ=255)
                                                                                                                                    cTGGCTGCTGCCCAGCAGCTGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTGGGCCAGGCAGCACCGCGCGCCACCAccc              >  6:320170/1‑149 (MQ=255)
                                                                                                                                       gCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTGGGCCAGGCAGCACCGCGCGCCACCACCcgcg           <  8:13144/149‑1 (MQ=255)
                                                                                                                                       gCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTGGGCCAGGCAGCACCGCGCGCCACCACCcgcg           >  8:6093/1‑149 (MQ=255)
                                                                                                                                         aGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTGGGCCAGGCAGCACCGCGCGCCACCACCCGCGAc         <  3:173810/149‑1 (MQ=255)
                                                                                                                                          gctgctGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTGGGCCAGGCAGCACCGCGCGCCACCACCCGCGAcc        >  8:27842/1‑149 (MQ=255)
                                                                                                                                                gggCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTGGGCCAGGCAGCACCGCGCGCCACCACCCGCGACCTgcgcg  >  4:142880/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
ACTGATGCCGTGCGACATCAGCAAGGCTGAGCGCGGTTACCGCCTGCTGGAACCGGACAGCTACGAGCTGGACCTGCCGCACACCGCGCTGGACGCCGAAATGGCCCTGGACGACGCCTATGCCGCCCTCGCCTGGCAGCTGCTGGGCAGCAGCCAGCGCCTGGCCGAGCAGTACCGCCACAACGGCCACTGG‑‑CGCCTGCTCCAGGACTGGGTCGAGTGCCTGAGCGAGCTGCGCGCGCTCACCGCCAGCCTGGGCCAGGCAGCACCGCGCGCCACCACCCGCGACCTGCGCG  >  NC_002947/5917506‑5917798

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: