Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,049,443 A→T 10.3% K102M (AAG→ATG) ‡ spoT → bifunctional (p)ppGpp synthase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,049,4430AT10.3% 162.4 / 7.1 58K102M (AAG→ATG) ‡spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base T (2/4);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAG  >  NC_002947/6049296‑6049583
                                                                                                                                                   |                                                                                                                                            
ccAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAg                                                                                                                                             >  6:154430/1‑149 (MQ=255)
 cAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCaa                                                                                                                                              <  1:24604/147‑1 (MQ=255)
 cAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCaa                                                                                                                                              >  2:24604/1‑147 (MQ=255)
 cAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGc                                                                                                                                            >  2:69151/1‑149 (MQ=255)
  aGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCt                                                                                                                                           >  5:338913/1‑149 (MQ=255)
        cTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCTTGCTGAccc                                                                                                                                      <  5:108961/148‑1 (MQ=255)
           gCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCTTGCTGAccccatc                                                                                                                                  <  3:364982/149‑5 (MQ=255)
           gCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGa                                                                                                                                   <  4:111914/148‑1 (MQ=255)
           gCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGTCAGCTTGCTGAccccatc                                                                                                                                  <  6:191489/148‑5 (MQ=255)
                    cacaTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                         <  4:340195/149‑1 (MQ=255)
                    cacaTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                         <  8:277292/149‑1 (MQ=255)
                    cacaTGGACCATCAGAGCCTGATGGCGGCCATGCTCCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                         <  6:324845/149‑1 (MQ=255)
                     acaTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTc                                                                                                                           >  2:299952/1‑146 (MQ=255)
                     acaTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                         <  1:299952/148‑1 (MQ=255)
                       aTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGAc                                                                                                                      >  2:231541/1‑149 (MQ=255)
                       aTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGAc                                                                                                                      <  1:69151/149‑1 (MQ=255)
                           aCCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAg                                                                                                                  <  6:374560/149‑1 (MQ=255)
                                                 cATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCGGGCGGAAAACt                                                                                            >  7:66690/1‑149 (MQ=255)
                                                    gctgcACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGgcgc                                                                                                        <  6:148871/134‑1 (MQ=255)
                                                    gctgcACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACAGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACAAAGGCCGAGgcgc                                                                                                        >  5:148871/1‑134 (MQ=255)
                                                     ctgcACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTcc                                                                                         <  2:155161/148‑1 (MQ=255)
                                                        cACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCagaa                                                                                     <  6:33702/149‑1 (MQ=255)
                                                           gttgtgATCGAAGACACCGGCATCGCCAAGGAAGCCCTAAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGat                                                                                  <  8:319038/147‑1 (MQ=255)
                                                           gATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGat                                                                                  >  3:181040/1‑149 (MQ=255)
                                                               tgATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggcc                                                                              <  1:231541/149‑1 (MQ=255)
                                                                 aTCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccat                                                                            >  6:22918/1‑149 (MQ=255)
                                                                    gAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggc                                                                         >  4:113536/1‑149 (MQ=255)
                                                                      aGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggcca                                                                       >  7:96769/1‑149 (MQ=255)
                                                                        acacCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatg                                                                     <  6:279859/149‑1 (MQ=255)
                                                                          acCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                         <  2:164740/95‑1 (MQ=255)
                                                                          acCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                         >  1:164740/1‑95 (MQ=255)
                                                                          acCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAg                                                                                                      <  7:95784/114‑1 (MQ=255)
                                                                          acCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAg                                                                                                      >  8:95784/1‑114 (MQ=255)
                                                                                 tCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGat                                                            <  8:96769/149‑1 (MQ=255)
                                                                                 tCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGa                                                             >  6:322694/1‑148 (MQ=255)
                                                                                    ccAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGTGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCATGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATc                                                         >  7:172430/1‑149 (MQ=255)
                                                                                     cAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATcc                                                        <  3:113536/149‑1 (MQ=255)
                                                                                      aaGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCcg                                                       >  7:8236/1‑149 (MQ=255)
                                                                                       aGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCcgc                                                      <  6:338913/149‑1 (MQ=255)
                                                                                        ggAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCcgcg                                                     >  4:256145/1‑149 (MQ=255)
                                                                                         gAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGt                                                    <  5:322694/149‑1 (MQ=255)
                                                                                            gCCCTCAGCAAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCTATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGAt                                                 <  7:240276/149‑1 (MQ=255)
                                                                                                  agccagcAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTgg                                            <  8:8236/148‑1 (MQ=255)
                                                                                                  agccagcAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGt                                           <  3:256145/149‑1 (MQ=255)
                                                                                                          gTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCGAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTCGTCAAGCTgg                                   >  6:317055/1‑149 (MQ=255)
                                                                                                             tGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGcc                                 <  5:154430/148‑1 (MQ=255)
                                                                                                               gCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGAc                              >  4:280854/1‑149 (MQ=255)
                                                                                                                    aCCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCagaa                                                                                     >  4:219515/1‑89 (MQ=255)
                                                                                                                    aCCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCagaa                                                                                     <  3:219515/89‑1 (MQ=255)
                                                                                                                      cGTGGCCGAATTGGTCGATGGGGTCGGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTg                        >  6:208709/1‑148 (MQ=255)
                                                                                                                           ccGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCacaa                   >  1:210926/1‑148 (MQ=255)
                                                                                                                                   gTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCa           <  2:210926/148‑1 (MQ=255)
                                                                                                                                       tctgggTCAGCTTGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAg                                       <  7:308042/113‑1 (MQ=255)
                                                                                                                                       tctgggTCAGCTTGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAg                                       >  8:308042/4‑116 (MQ=255)
                                                                                                                                       tctgggTCAGCTTGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTg      >  7:337985/4‑149 (MQ=255)
                                                                                                                                          gggTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGaa   <  5:22918/149‑1 (MQ=255)
                                                                                                                                           ggTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTg                                             >  4:322413/1‑106 (MQ=255)
                                                                                                                                           ggTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTg                                             <  3:322413/106‑1 (MQ=255)
                                                                                                                                           ggTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGAGATATCCGCGTGATCCTGGTCAAGCTGGACGACCGCCTGCACAACATGCGCACACTGGAAg  >  1:380495/1‑149 (MQ=255)
                                                                                                                                            gTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatgg                                                                          >  3:215979/1‑76 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
CCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAG  >  NC_002947/6049296‑6049583

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: