Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,120,767 T→C 6.9% D263G (GAC→GGC)  wcaJ ← colanic biosynthesis UDP‑glucose lipid carrier transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,120,7670TC 6.9% 174.4 / 2.9 59D263G (GAC→GGC) wcaJcolanic biosynthesis UDP‑glucose lipid carrier transferase
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base C (2/2);  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 = 9.38e-01

AAAATAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCG  >  NC_000913/2120626‑2120904
                                                                                                                                             |                                                                                                                                         
aaaaTAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAg                                                                                                                                    >  4:244518/1‑149 (MQ=255)
  aaTAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCg                                                                                                                                  <  8:132367/149‑1 (MQ=255)
  aaTAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCg                                                                                                                                  <  3:216812/149‑1 (MQ=255)
         ggCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGcaccac                                                                                                                           <  6:377098/149‑1 (MQ=255)
          gCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGGGTCGTAAAGCGTCacgacc                                                                                                                          >  4:70173/1‑149 (MQ=255)
           cccTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCg                                                                                                                         <  3:129263/149‑1 (MQ=255)
              tGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCg                                                                                                                                  >  4:246823/1‑137 (MQ=255)
              tGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCg                                                                                                                                  <  3:246823/137‑1 (MQ=255)
                   aaCTGAGTTTCACCGCCATCGCAATACAGCCCAGCACCGGGGAGATCCGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGTATGCGGTTAACCCCGGAAAGCGGCGTTGCGTAAATCGGCACCACCGGGACGCCg                                                                                                                 <  1:100971/149‑1 (MQ=255)
                      tGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGCCGTAAAGCGGCACCACCGGTACGCCGtt                                                                                                               >  7:229570/1‑148 (MQ=255)
                             cACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTACCACCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCCTCTCtt                                                                                                       >  2:218728/1‑149 (MQ=255)
                                    aGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGGTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCg                                                                                                >  6:324846/1‑149 (MQ=255)
                                     gcgcAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGt                                                                                               <  5:349464/149‑1 (MQ=255)
                                     gcgcAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGt                                                                                               <  3:190847/149‑1 (MQ=255)
                                        cAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGCCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGaa                                                                                            <  8:229570/149‑1 (MQ=255)
                                          aTACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATg                                                                                          <  8:159178/149‑1 (MQ=255)
                                          aTACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATg                                                                                          >  4:80175/1‑149 (MQ=255)
                                              agcacagcacCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGgag                                                                                       >  8:294302/1‑148 (MQ=255)
                                                cacagcacCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAAt                                                                                    <  5:333007/149‑1 (MQ=255)
                                                      acCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGaa                                                                              <  7:222018/149‑1 (MQ=255)
                                                            gagaTCAGCTGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTGTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCCGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTaaa                                                                        <  4:353399/149‑1 (MQ=255)
                                                               aTCAGAAGCAGAATAAGCGTCGCCAGCACAATCTCTTCCGCACGGTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGa                                                                      <  2:297602/148‑1 (MQ=255)
                                                                  agcagcagAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTc                                                                  >  3:63119/1‑149 (MQ=255)
                                                                   gcagcagAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCg                                                                 >  1:309882/1‑149 (MQ=255)
                                                                    cagcagAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCgg                                                                >  4:260241/1‑149 (MQ=255)
                                                                      gcagAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGCCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCggg                                                               >  6:304529/1‑148 (MQ=255)
                                                                        agAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGAt                                                            <  4:63119/149‑1 (MQ=255)
                                                                            tAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATcagc                                                        >  2:244335/1‑149 (MQ=255)
                                                                            tAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATcag                                                         >  3:51535/1‑148 (MQ=255)
                                                                               gCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATcagcag                                                      >  8:133560/1‑148 (MQ=255)
                                                                               gCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATcagcag                                                      >  3:93861/1‑148 (MQ=255)
                                                                                    gCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGcc                                                                                                                  <  7:199735/83‑1 (MQ=255)
                                                                                    gCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGcc                                                                                                                  >  8:199735/1‑83 (MQ=255)
                                                                                           caATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAggtg                                         <  7:133560/149‑1 (MQ=255)
                                                                                           caATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAggtg                                         <  6:98640/149‑1 (MQ=255)
                                                                                            aaTGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAggtgg                                        >  1:174763/1‑149 (MQ=255)
                                                                                                  ttCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGaa                                                                              <  8:35159/105‑1 (MQ=255)
                                                                                                  ttCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGaa                                                                              >  7:35159/1‑105 (MQ=255)
                                                                                                  ttCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCg                                  >  5:374631/1‑149 (MQ=255)
                                                                                                  ttCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCg                                  >  3:259678/1‑149 (MQ=255)
                                                                                                      gCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAg                              >  1:36557/1‑149 (MQ=255)
                                                                                                      gCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAg                              <  3:244518/149‑1 (MQ=255)
                                                                                                          gTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTaaa                                                                        >  8:80973/1‑103 (MQ=255)
                                                                                                          gTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTaaa                                                                        <  7:80973/103‑1 (MQ=255)
                                                                                                          gTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGa                          <  5:154102/149‑1 (MQ=255)
                                                                                                           tttGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGCCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGAt                         <  5:304529/149‑1 (MQ=255)
                                                                                                             tGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATg                        <  4:93861/148‑1 (MQ=255)
                                                                                                              gAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATgg                       <  2:36557/148‑1 (MQ=255)
                                                                                                                gCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAggtg                                         <  5:152790/128‑1 (MQ=255)
                                                                                                                gCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAggtg                                         >  6:152790/1‑128 (MQ=255)
                                                                                                                 cAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGCACAGGTGGGGTCCGCCAGGTGATGGAc                     >  8:149152/1‑147 (MQ=255)
                                                                                                                  aGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAg                  >  6:190437/1‑149 (MQ=255)
                                                                                                                       gtgaaCCCCGGAAAGCGGCGGGTCGTACAGAGGCACCACCGGTACGCCGTTCATCGCTTCGACGCGTGAATGGAGAAGGTTGAAGGTAAAGACGGCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGCCCAGttttt             <  2:163486/146‑1 (MQ=255)
                                                                                                                        ttAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGtttttt            >  3:358778/1‑149 (MQ=255)
                                                                                                                           aCCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCa          <  7:149152/148‑1 (MQ=255)
                                                                                                                                 gAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTcgcgc   <  3:80175/149‑1 (MQ=255)
                                                                                                                                  aaaGCGGCGTGTCGTTAAGCTGCTCGACCGGTACGACGTTCATCTCTTCGAGGCGTGAATGGAGAATTTTGATGGTAAAGACTTAGGTGATCAGCAGCACCGTACTGGTGGTGTCCGACAGTTGATGGAGCAGTTTTTTctctatagat  <  2:309882/149‑7 (MQ=255)
                                                                                                                                  aaaGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCAc                                                   <  4:225990/100‑1 (MQ=255)
                                                                                                                                  aaaGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCAc                                                   >  3:225990/1‑100 (MQ=255)
                                                                                                                                             |                                                                                                                                         
AAAATAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCG  >  NC_000913/2120626‑2120904

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: