Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,158,617 T→C L43P (CTG→CCG)  oxyR → oxidative and nitrosative stress transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,158,6170TC96.9% 218.3 / ‑2.8 65L43P (CTG→CCG) oxyRoxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (31/32);  total (32/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 = 7.03e-01

ATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGC  >  NC_000913/4158483‑4158747
                                                                                                                                      |                                                                                                                                  
atggatAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGc                                                                                                                                >  1:1576259/1‑139 (MQ=255)
   gatAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTg                                                                                                                             >  2:1628121/1‑139 (MQ=255)
   gatAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGACTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTg                                                                                                                             >  1:1356998/1‑139 (MQ=255)
     tAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGACTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGAt                                                                                                                           <  2:1356998/139‑1 (MQ=255)
      aaTGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATg                                                                                                                          >  1:584879/1‑139 (MQ=255)
      aaTGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATg                                                                                                                          <  2:743649/139‑1 (MQ=255)
                               cTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGAc                                                                                                            >  1:406360/1‑128 (MQ=255)
                               cTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGAc                                                                                                            <  2:406360/128‑1 (MQ=255)
                               cTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAg                                                                                                 >  1:1345001/1‑139 (MQ=255)
                               cTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAg                                                                                                 <  2:991117/139‑1 (MQ=255)
                                 ggtggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCg                                                                                                                               <  2:1662484/107‑1 (MQ=255)
                                 ggtggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCg                                                                                                                               >  1:1662484/1‑107 (MQ=255)
                                 ggtggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAACCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAgtg                                                                                               >  1:1678646/1‑139 (MQ=255)
                                   tggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGtgtt                                                                                             >  1:451533/1‑139 (MQ=255)
                                         tGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGc                                                                                                                      <  2:908503/108‑1 (MQ=255)
                                         tGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGc                                                                                                                      >  1:908503/1‑108 (MQ=255)
                                         tGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCAc                                                                                       <  2:1345001/139‑1 (MQ=255)
                                         tGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCAc                                                                                       >  1:955618/1‑139 (MQ=255)
                                          ggCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCAcc                                                                                      <  1:911204/139‑1 (MQ=255)
                                            cTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCa                                                                                    <  2:955618/139‑1 (MQ=255)
                                                 cacCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGTCCAGCCGTAAAGTGTTGTTCACCCAGGCgg                                                                               >  1:785598/1‑139 (MQ=255)
                                                    cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGc                                                                                                                      >  1:412175/1‑97 (MQ=255)
                                                    cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGc                                                                                                                      <  2:412175/97‑1 (MQ=255)
                                                    cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGc                                                                                                                >  2:1288373/1‑103 (MQ=255)
                                                    cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGc                                                                                                                <  1:1288373/103‑1 (MQ=255)
                                                    cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCAcc                                                                                      >  1:1109638/1‑129 (MQ=255)
                                                    cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCAcc                                                                                      <  2:1109638/129‑1 (MQ=255)
                                                    cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGaa                                                                            >  1:798753/1‑139 (MQ=255)
                                                    cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGaa                                                                            >  1:564102/1‑139 (MQ=255)
                                                    cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGaa                                                                            >  1:301935/1‑139 (MQ=255)
                                                    cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGaa                                                                            >  2:1356735/1‑139 (MQ=255)
                                                      ccATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAATTTCGTAAGCTGGAAGATGAGCCGGGCGTGa                                                                                                                            >  2:786195/1‑89 (MQ=255)
                                                      ccATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGa                                                                                                                            <  1:786195/89‑1 (MQ=255)
                                                      ccATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGTCCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAAtg                                                                          <  2:785598/139‑1 (MQ=255)
                                                      ccATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAAtg                                                                          <  1:1356735/139‑1 (MQ=255)
                                                       cATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAAtgc                                                                         <  1:1658994/139‑1 (MQ=255)
                                                                          aTTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGcgcg                                                      >  2:10629/1‑139 (MQ=255)
                                                                            tCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGAc                                                                                                            <  1:1166692/83‑1 (MQ=255)
                                                                            tCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGAc                                                                                                            >  2:1166692/1‑83 (MQ=255)
                                                                            tCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTa                                                    >  2:571497/1‑139 (MQ=255)
                                                                            tCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTa                                                    >  2:1547108/1‑139 (MQ=255)
                                                                                    cGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGtgctg                                            >  1:849103/1‑139 (MQ=255)
                                                                                      ttAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTa                                                    >  1:1657789/1‑129 (MQ=255)
                                                                                      ttAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTa                                                    <  2:1657789/129‑1 (MQ=255)
                                                                                               cGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGTTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGaa                                 <  2:1576259/139‑1 (MQ=255)
                                                                                               cGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGaa                                 <  1:1547108/139‑1 (MQ=255)
                                                                                               cGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAACCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGaa                                 <  2:1678646/139‑1 (MQ=255)
                                                                                                 aCGCATAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAg                               <  2:390672/139‑1 (MQ=255)
                                                                                                  cGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGt                              <  2:301935/139‑1 (MQ=255)
                                                                                                         cGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCtgg                                                                   <  2:921265/95‑1 (MQ=255)
                                                                                                         cGGGCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCtgg                                                                   >  1:921265/1‑95 (MQ=255)
                                                                                                           ggCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGAt                                                              <  1:1264285/98‑1 (MQ=255)
                                                                                                           ggCAAATTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGAt                                                              >  2:1264285/1‑98 (MQ=255)
                                                                                                              aaaTTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGAt                  <  2:451533/139‑1 (MQ=255)
                                                                                                              aaaTTCGTAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGAt                  <  2:564102/139‑1 (MQ=255)
                                                                                                                     tAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGaa                                 >  1:1376330/1‑117 (MQ=255)
                                                                                                                     tAAGCTGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGaa                                 <  2:1376330/117‑1 (MQ=255)
                                                                                                                          tGGAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCcagca      <  1:571497/139‑1 (MQ=255)
                                                                                                                            gAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGAc                                                                                                            >  1:326873/1‑35 (MQ=255)
                                                                                                                            gAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGAc                                                                                                            <  2:326873/35‑1 (MQ=255)
                                                                                                                            gAAGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAgg    >  2:730093/1‑139 (MQ=255)
                                                                                                                              aGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGc  <  1:4507/139‑1 (MQ=255)
                                                                                                                              aGATGAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGc  <  1:1628121/139‑1 (MQ=255)
                                                                                                                                  gAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTgttgtt                                                                                          >  1:1188839/1‑47 (MQ=255)
                                                                                                                                  gAGCCGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTgttgtt                                                                                          <  2:1188839/47‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
ATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGC  >  NC_000913/4158483‑4158747

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