Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,455,242 C→T 7.1% L117L (CTG→CTA PP_4792 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,455,2420CT 7.1% 165.8 / 2.0 56L117L (CTG→CTAPP_4792hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/27);  new base T (2/2);  total (27/29)
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.39e-01

TGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTC  >  NC_002947/5455103‑5455379
                                                                                                                                           |                                                                                                                                         
tGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCTGGCCCGAAGCCTTGGCCAGCcgcgcg                                                                                                                                  <  1:218128/149‑1 (MQ=255)
tGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCcgcgcg                                                                                                                                  <  8:156486/149‑1 (MQ=255)
 gCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCgg                                                                                                                                 <  7:279744/149‑1 (MQ=255)
  cGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCgg                                                                                                                                 >  5:167059/1‑148 (MQ=255)
        gCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAgcg                                                                                                                          >  7:113685/1‑149 (MQ=255)
        gCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAgcg                                                                                                                          >  5:325856/1‑149 (MQ=255)
               cgcgCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGt                                                                                                                   >  7:55524/1‑149 (MQ=255)
                gcgcCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGtt                                                                                                                  >  5:336539/1‑149 (MQ=255)
                  gcCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTc                                                                                                                 >  5:209354/1‑148 (MQ=255)
                   ccATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCa                                                                                                               <  4:59395/149‑1 (MQ=255)
                          aCCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTgccgt                                                                                                        >  1:329879/1‑149 (MQ=255)
                          aCCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTgccgt                                                                                                        >  7:7424/1‑149 (MQ=255)
                           ccTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTgccgtg                                                                                                       >  3:151009/1‑149 (MQ=255)
                                aCAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTgccgtggccgt                                                                                                  <  4:269810/149‑1 (MQ=255)
                                aCAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTgccgtggccgt                                                                                                  <  8:7424/149‑1 (MQ=255)
                                  aGCGAAATCTTCATGCACGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGccgtggccgtgg                                                                                                <  8:104165/149‑1 (MQ=255)
                                      aaaTCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTg                                                                                            <  6:336539/149‑1 (MQ=255)
                                             cATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTATTGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTc                                                                                                                 >  4:139314/1‑121 (MQ=255)
                                             cATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTc                                                                                                                 <  3:139314/121‑1 (MQ=255)
                                             cATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCAcc                                                                                      <  6:325856/148‑1 (MQ=255)
                                              aTGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCTAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGc                                                                                    >  5:146945/1‑149 (MQ=255)
                                              aTGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCTAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGc                                                                                    >  5:153234/1‑149 (MQ=255)
                                                 cAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCAc                                                                                 >  2:59020/1‑149 (MQ=255)
                                                      gTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAg                                                                            <  2:339392/149‑1 (MQ=255)
                                                          cGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCTAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCg                                                                        <  6:153234/149‑1 (MQ=255)
                                                          cGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCTAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCg                                                                        <  6:146945/149‑1 (MQ=255)
                                                           gCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGGCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCgg                                                                       <  3:113907/149‑1 (MQ=255)
                                                               tCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTcc                                                                   <  8:113685/149‑1 (MQ=255)
                                                                      cgcctcgcTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCgg                                                                       <  6:405757/138‑1 (MQ=255)
                                                                      cgcctcgcTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCgg                                                                       >  5:405757/1‑138 (MQ=255)
                                                                       gcctcgcTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGGCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCgcgc                                                             <  4:22347/147‑1 (MQ=255)
                                                                       gcctcgcTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCgcgc                                                             <  8:55524/147‑1 (MQ=255)
                                                                        cctcgcTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTgccgtggccgt                                                                                                  <  3:47148/109‑1 (MQ=255)
                                                                        cctcgcTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTgccgtggccgt                                                                                                  >  4:47148/1‑109 (MQ=255)
                                                                            gcTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCCCCGCCCCCCCAGGCCGGGTCCCCGCGCTCGCGcc                                                      >  1:401066/1‑149 (MQ=255)
                                                                            gcTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTcgcg                                                        <  2:401066/147‑1 (MQ=255)
                                                                                     cgTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCa                                             <  4:151009/149‑1 (MQ=255)
                                                                                      gTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAg                                            <  2:329879/149‑1 (MQ=255)
                                                                                            tAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGcacca                                       >  2:12692/1‑148 (MQ=255)
                                                                                                gCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGc                                           <  2:2296/140‑1 (MQ=255)
                                                                                                gCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGc                                           >  1:2296/1‑140 (MQ=255)
                                                                                                 cAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGca                                          >  8:366056/1‑140 (MQ=255)
                                                                                                 cAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGca                                          <  7:366056/140‑1 (MQ=255)
                                                                                                     ccaccaCTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATc                             >  4:186970/1‑149 (MQ=255)
                                                                                                         cacTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGTGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCgg                                                                       <  1:324010/103‑1 (MQ=255)
                                                                                                         cacTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCgg                                                                       >  2:324010/1‑103 (MQ=255)
                                                                                                              gTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGc                     <  6:167059/148‑1 (MQ=255)
                                                                                                              gTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCt                    >  7:348303/1‑149 (MQ=255)
                                                                                                                  acacTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGcc                >  6:165125/1‑149 (MQ=255)
                                                                                                                   cacaGGCAGGCCCGAAGCCTTGGCCCGCCGCGCCGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGACGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGc                 >  2:242835/1‑147 (MQ=255)
                                                                                                                   cacTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCg                                                    >  3:37260/1‑112 (MQ=255)
                                                                                                                   cacTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCg                                                    <  4:37260/112‑1 (MQ=255)
                                                                                                                   cacTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCa               >  7:331299/1‑149 (MQ=255)
                                                                                                                        gCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAggccg          <  3:186970/149‑1 (MQ=255)
                                                                                                                           ggCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCATCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAggccggcc       <  6:209354/149‑1 (MQ=255)
                                                                                                                                gAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCACGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTc  >  4:373832/1‑149 (MQ=255)
                                                                                                                                          ccAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGcacca                                       <  5:41980/102‑1 (MQ=255)
                                                                                                                                          ccAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCAGGGTCCCCGCGCTCGCGCCCGGCGCCCAGcacca                                       >  6:41980/1‑102 (MQ=255)
                                                                                                                                           |                                                                                                                                         
TGCGGCTGGCCTCTTCGCGCCATTTCACCTCGACAGCGAAATCTTCATGCAAGCGTTCCGCCATCAGCCTCGCCTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTC  >  NC_002947/5455103‑5455379

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: