Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,751 C→T A451V (GCT→GTT)  spoT → bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,823,7510CT100.0% 210.9 / NA 55A451V (GCT→GTT) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/26);  total (29/26)

CGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGC  >  NC_000913/3823613‑3823897
                                                                                                                                          |                                                                                                                                                  
cGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGc                                                                                                                                        >  1:225854/1‑151 (MQ=255)
cGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGc                                                                                                                                        <  2:396519/151‑1 (MQ=255)
    cGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCt                                                                                                                                    >  2:246705/1‑151 (MQ=255)
       cGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTgg                                                                                                                                 <  1:188028/151‑1 (MQ=255)
       cGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTgg                                                                                                                                 >  2:594561/1‑151 (MQ=255)
       cGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTgg                                                                                                                                 <  1:467778/151‑1 (MQ=255)
        gACTTCGATTATGCAGTGCATACCGATATCGGTCATGTCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGATTGGc                                                                                                                                >  1:230870/1‑151 (MQ=255)
            tCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGGATGCCGCTTGGCTGa                                                                                                                             >  1:582970/1‑150 (MQ=255)
                 tATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACtt                                                                                                                         <  2:273853/149‑1 (MQ=255)
                      aGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGtt                                                                                                                  <  1:228483/151‑1 (MQ=255)
                       gTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGtt                                                                                                                  >  1:59138/1‑150 (MQ=255)
                          cATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGc                                                                                                               >  1:520441/1‑150 (MQ=255)
                           aTACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTc                                                                                                             >  1:65170/1‑151 (MQ=255)
                             aCCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCg                                                                                                            <  1:351574/150‑1 (MQ=255)
                               cGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGaa                                                                                                          >  1:412954/1‑150 (MQ=255)
                                gATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGaa                                                                                                          >  2:334821/1‑149 (MQ=255)
                                    tCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcg                                                                                                    >  1:192286/1‑151 (MQ=255)
                                         cATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCa                                                                                               <  2:658538/151‑1 (MQ=255)
                                            gCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAGACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCaaa                                                                                             >  1:580964/1‑150 (MQ=255)
                                            gCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCaaa                                                                                             <  1:521644/150‑1 (MQ=255)
                                             ccTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAAt                                                                                           <  2:652276/151‑1 (MQ=255)
                                             ccTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAAt                                                                                           <  2:166435/151‑1 (MQ=255)
                                                gCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTc                                                                                         <  2:527636/150‑1 (MQ=255)
                                                         cACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCt                                                                               >  1:314025/1‑151 (MQ=255)
                                                            gcgTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCCAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGACTGCCGCTTGGCTGAACTTTGTCGTTAGCTCCAAAGCGCGCGCCAAAATTCGTCAGTTGCTGa                                                                             >  1:162176/1‑150 (MQ=255)
                                                              gTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGaaaa                                                                          <  2:143026/151‑1 (MQ=255)
                                                                       cAGCCTTCCCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGc                                                                 <  1:10323/151‑1 (MQ=255)
                                                                        aGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCg                                                                >  1:582406/1‑151 (MQ=255)
                                                                               cccGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTgatgat                                                         >  1:124259/1‑151 (MQ=255)
                                                                                ccGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGAtt                                                        >  2:515970/1‑151 (MQ=255)
                                                                                  gCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCt                                                      <  2:700502/151‑1 (MQ=255)
                                                                                    tGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGt                                                    >  1:493391/1‑151 (MQ=255)
                                                                                    tGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGt                                                    >  2:319341/1‑151 (MQ=255)
                                                                                     gTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGt                                                    <  1:612347/150‑1 (MQ=255)
                                                                                     gTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGt                                                    >  1:656922/1‑150 (MQ=255)
                                                                                      tCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTa                                                   >  1:310213/1‑150 (MQ=255)
                                                                                       cGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAg                                                 <  1:477787/151‑1 (MQ=255)
                                                                                              gCTTACCAGCGGTCAAACCGTTGTAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATGCTGTAAGTCt                                              <  1:478828/147‑1 (MQ=255)
                                                                                              gCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGc                                          >  1:355803/1‑151 (MQ=255)
                                                                                              gCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCACCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGc                                          >  1:173068/1‑151 (MQ=255)
                                                                                                 tACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGcc                                         <  1:616364/149‑1 (MQ=255)
                                                                                                     aGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCgtcgtc                                   >  2:228723/1‑151 (MQ=255)
                                                                                                       cGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTg                                 <  2:329714/151‑1 (MQ=255)
                                                                                                           cAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTc                              >  1:523377/1‑150 (MQ=255)
                                                                                                           cAAACCGTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTc                              >  1:272971/1‑150 (MQ=255)
                                                                                                           cAAACCGTTGAAAGCATTACCGCTCCGGGCGTTCGCCCGAATTCCGCTTGGCTGATCTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATGCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTATGTAAGCCTGGGACGGCGTCt                                  <  1:476430/146‑1 (MQ=255)
                                                                                                                 gTTGAAATCATTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATg                       <  2:497695/151‑1 (MQ=255)
                                                                                                                          aTTACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGggtg              <  2:638872/151‑1 (MQ=255)
                                                                                                                            tACCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGggtggt            >  1:206755/1‑151 (MQ=255)
                                                                                                                             aCCGCTCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTa           <  2:650899/151‑1 (MQ=255)
                                                                                                                                  tCCGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCTAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATGCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTCGggtggt            <  1:637195/145‑1 (MQ=255)
                                                                                                                                    cGGGCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTa     >  1:80272/1‑150 (MQ=255)
                                                                                                                                     gggCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAg   <  2:597571/151‑1 (MQ=255)
                                                                                                                                      ggCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGc  >  2:410726/1‑151 (MQ=255)
                                                                                                                                      ggCGTTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGc  <  2:225901/151‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
CGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGC  >  NC_000913/3823613‑3823897

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