Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,445,402 G→C 19.0% A51A (GCC→GCG lexA‑I ← transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,445,4020GC19.0% 132.9 / 29.1 58A51A (GCC→GCGlexA‑Itranscriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (20/27);  new base C (5/6);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-01

TGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTTGGCGTGGCGTC  >  NC_002947/2445265‑2445543
                                                                                                                                         |                                                                                                                                             
tgCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCCCGCAAGGGcgcga                                                                                                                                   <  2:351517/150‑5 (MQ=255)
 gCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGa                                                                                                                                  <  2:706800/150‑1 (MQ=255)
 gCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGa                                                                                                                                  <  2:256753/150‑1 (MQ=255)
  cTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAg                                                                                                                                 >  2:1058915/1‑150 (MQ=255)
  cTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAg                                                                                                                                 <  1:700391/150‑1 (MQ=255)
   tGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTg                                                                                                                                   >  1:640465/1‑147 (MQ=255)
   tGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTg                                                                                                                                   <  2:640465/147‑1 (MQ=255)
   tGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAgg                                                                                                                                <  1:562656/150‑1 (MQ=255)
        ggCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGctc                                                                                                                           >  2:881897/1‑150 (MQ=255)
        ggCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGctc                                                                                                                           <  1:1058915/150‑1 (MQ=255)
         gCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGctcc                                                                                                                          >  1:184768/1‑150 (MQ=255)
              aaTCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGTATGCCGCGGTAGCCGGCTGGCTGGATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTcggc                                                                                                                     >  2:467706/1‑150 (MQ=255)
              aaTCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTcggc                                                                                                                     <  2:1124343/150‑1 (MQ=255)
                     gcACCGGCAGCGTCTCGGCCGATGATGTGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTg                                                                                                              >  2:258184/1‑150 (MQ=255)
                        ccGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCATGCCAGGTACGCGCTGGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGc                                                                                                           <  1:664604/150‑1 (MQ=255)
                                   ccGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGcc                                                                                                <  2:360735/150‑1 (MQ=255)
                                   ccGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCCCGCAAGGGc                                                                                                                                       >  2:1129278/1‑111 (MQ=255)
                                   ccGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCCCGCAAGGGc                                                                                                                                       <  1:1129278/111‑1 (MQ=255)
                                             atgGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCCCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGa                                                                                      <  2:950108/150‑1 (MQ=255)
                                              tgGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAg                                                                                     <  2:184768/150‑1 (MQ=255)
                                                ggCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAggg                                                                                                                                        >  1:80042/1‑97 (MQ=255)
                                                ggCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAggg                                                                                                                                        <  2:80042/97‑1 (MQ=255)
                                                 gCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCATTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCaa                                                                                  >  1:22476/1‑150 (MQ=255)
                                                  cAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAAt                                                                                 <  1:467706/150‑1 (MQ=255)
                                                                cAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCgg                                                                       >  1:57902/1‑146 (MQ=255)
                                                                cAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCgg                                                                       <  2:57902/146‑1 (MQ=255)
                                                                cAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAg                                                                   <  1:1037369/150‑1 (MQ=255)
                                                                     ttGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCggg                                                              >  2:632501/1‑150 (MQ=255)
                                                                          ttCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCGTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGcc                                                         >  1:473819/1‑150 (MQ=255)
                                                                           tCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCg                                                        >  1:233990/1‑150 (MQ=255)
                                                                             caggccaggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTGGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGTAATCTCGGCGCGGGTAGTCGGGAAGCCGtt                                                      <  1:881897/150‑1 (MQ=255)
                                                                             caggccaggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTTGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGtt                                                      >  1:210888/1‑150 (MQ=255)
                                                                             caggccaggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGtt                                                      <  2:22476/150‑1 (MQ=255)
                                                                             caggccaggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCATTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGtt                                                      <  1:258184/150‑1 (MQ=255)
                                                                                  caggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCgg                                                               >  1:953564/1‑136 (MQ=255)
                                                                                  caggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCgg                                                               <  2:953564/136‑1 (MQ=255)
                                                                                           ggatgcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCACACCGGGGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAgc                                        <  1:80024/150‑1 (MQ=255)
                                                                                            gatgcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAAt                                                                                 <  2:144143/108‑1 (MQ=255)
                                                                                            gatgcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAAt                                                                                 >  1:144143/1‑108 (MQ=255)
                                                                                               gcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCgg                                                                       <  1:604017/115‑1 (MQ=255)
                                                                                               gcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCgg                                                                       >  2:604017/1‑115 (MQ=255)
                                                                                                 cgcgGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGCCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTtga                                  <  1:632501/150‑1 (MQ=255)
                                                                                                          cgcCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCga                         <  2:1274330/150‑1 (MQ=255)
                                                                                                                     aTTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCt              >  1:1080950/1‑150 (MQ=255)
                                                                                                                      tttCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCtt             <  2:233990/150‑1 (MQ=255)
                                                                                                                           aTCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAgc                                        >  1:564719/1‑118 (MQ=255)
                                                                                                                           aTCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAgc                                        <  2:564719/118‑1 (MQ=255)
                                                                                                                           aTCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTtggcgt        >  2:962033/1‑150 (MQ=255)
                                                                                                                             aaggcCCTTGCCCGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGcc                                                         <  2:1012539/96‑1 (MQ=255)
                                                                                                                             aaggcCCTTGCCCGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGcc                                                         >  1:1012539/4‑99 (MQ=255)
                                                                                                                             aaggcCCTTGCCCGCAAGGGCCTTGAGGCGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGtt                                                      >  2:486948/4‑102 (MQ=255)
                                                                                                                             aaggcCCTTGCCCGCAAGGGCCTTGAGGCGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGtt                                                      <  1:486948/99‑1 (MQ=255)
                                                                                                                               ggcCCTTGCCCGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCa                                              >  2:755117/2‑108 (MQ=255)
                                                                                                                               ggcCCTTGCCCGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCa                                              <  1:755117/107‑1 (MQ=255)
                                                                                                                                gcCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTTggcgtggcgt   >  1:806415/1‑150 (MQ=255)
                                                                                                                                gcCCTTGCCCGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCtt             >  2:673159/1‑140 (MQ=255)
                                                                                                                                gcCCTTGCCCGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCtt             <  1:673159/140‑1 (MQ=255)
                                                                                                                                 cccTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTTGGCGTGGCGTc  <  2:473819/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                             
TGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTTGGCGTGGCGTC  >  NC_002947/2445265‑2445543

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: