Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,490,375 C→T A255A (GCG→GCA idnR ← transcriptional repressor, 5‑gluconate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,490,3750CT100.0% 214.0 / NA 66A255A (GCG→GCAidnRtranscriptional repressor, 5‑gluconate‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/42);  total (24/42)

CCGACCAATATCAAAACGGGGAGTAATTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGCGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATGGATGAGATGGCGCGTGGATTCATGCGTAATGGGGAAAGA  >  NC_000913/4490240‑4490510
                                                                                                                                       |                                                                                                                                       
ccGACCAATATCAATACGGGGAGTAATTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTACAGCCCATGAAAACCCGCAATTGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCa                                                                                                                                      <  1:1015142/139‑1 (MQ=255)
ccGACCAATATCAAAACGGGGAGTAATTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCa                                                                                                                                      <  2:682923/139‑1 (MQ=255)
ccGACCAATATCAAAACGGGGAGTAATTACGCTGGCAAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCa                                                                                                                                      <  2:1199246/139‑1 (MQ=255)
  gACCAATATCAAAACGGGGAGTAATTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCAcc                                                                                                                                    <  1:362946/139‑1 (MQ=255)
         atCAAAACGGGGAGTAATTACGCTGGCGAGGCTTGGGATCATTCGCCTGCCAATCTCCAGCCCATGAAAACCCGCACTGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCa                                                                                                                             <  2:1075258/139‑1 (MQ=255)
                  gggAGTAATTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCa                                                                                                                    >  2:576121/1‑139 (MQ=255)
                  gggAGTAATTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCa                                                                                                                    >  2:725320/1‑139 (MQ=255)
                   ggAGTAATTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCAt                                                                                                                   <  1:725320/139‑1 (MQ=255)
                   ggAGTAATTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCAt                                                                                                                   <  1:219687/139‑1 (MQ=255)
                   ggAGTAATTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCAt                                                                                                                   >  1:135916/1‑139 (MQ=255)
                      gTAATTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTgg                                                                                                                >  1:317638/1‑139 (MQ=255)
                       tAATTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGt                                                                                                               <  1:1406337/139‑1 (MQ=255)
                       tAATTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGt                                                                                                               <  1:441661/139‑1 (MQ=255)
                         aTTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGc                                                                                                             <  1:3955/139‑1 (MQ=255)
                               cTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATa                                                                                                       >  1:515072/1‑139 (MQ=255)
                                 ggCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATaca                                                                                                     <  2:135916/139‑1 (MQ=255)
                                 ggCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATaca                                                                                                     <  1:18113/139‑1 (MQ=255)
                                  gCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATacac                                                                                                    >  2:1008464/1‑139 (MQ=255)
                                  gCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATacac                                                                                                    >  1:133795/1‑139 (MQ=255)
                                       gCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCt                                                                                               <  1:159589/139‑1 (MQ=255)
                                                 aTTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCAt                                                                                                                                 >  2:1008141/1‑95 (MQ=255)
                                                 aTTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCAt                                                                                                                                 <  1:1008141/95‑1 (MQ=255)
                                                 aTTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGt                                                                                                 <  1:949403/127‑1 (MQ=255)
                                                 aTTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGt                                                                                                 >  2:949403/1‑127 (MQ=255)
                                                 aTTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGt                                                                                     <  2:711706/139‑1 (MQ=255)
                                                          ccAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCAcc                                                                                                                                    <  2:1192008/83‑1 (MQ=255)
                                                          ccAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCAcc                                                                                                                                    >  1:1192008/1‑83 (MQ=255)
                                                          ccAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTaa                                                                            >  2:445112/1‑139 (MQ=255)
                                                          ccAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTaa                                                                            >  2:1544574/1‑139 (MQ=255)
                                                                tcCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCAt                                                                      <  1:1544574/139‑1 (MQ=255)
                                                                tcCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTACAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCAt                                                                      <  1:1008464/139‑1 (MQ=255)
                                                                 ccAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATc                                                                     <  1:1653740/139‑1 (MQ=255)
                                                                  cAGCCCATGAAAGCCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCa                                                                    <  2:1576928/139‑1 (MQ=255)
                                                                   aGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCAc                                                                   <  2:632865/139‑1 (MQ=255)
                                                                     cccATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGc                                                                 >  1:1366306/1‑139 (MQ=255)
                                                                        aTGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCAtt                                                              <  2:1266831/139‑1 (MQ=255)
                                                                         tGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTa                                                             >  2:1496690/1‑139 (MQ=255)
                                                                         tGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTa                                                             >  1:975744/1‑139 (MQ=255)
                                                                         tGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTa                                                             >  2:849287/1‑139 (MQ=255)
                                                                          gAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTaa                                                            >  1:1491469/1‑139 (MQ=255)
                                                                             aaCCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTg                                                         <  2:1491469/139‑1 (MQ=255)
                                                                               cccGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCa                                                       <  1:879357/139‑1 (MQ=255)
                                                                                ccGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCAt                                                      <  1:473858/139‑1 (MQ=255)
                                                                                   cAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCa                                                                                                                             >  1:582855/1‑65 (MQ=255)
                                                                                   cAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCa                                                                                                                             <  2:582855/65‑1 (MQ=255)
                                                                                   cAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTcc                                                   >  1:374957/1‑139 (MQ=255)
                                                                                     aTGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCAcc                                                                                                                                    >  2:1596114/1‑56 (MQ=255)
                                                                                     aTGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCAcc                                                                                                                                    <  1:1596114/56‑1 (MQ=255)
                                                                                            tttGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAaatgaat                                          <  1:576121/139‑1 (MQ=255)
                                                                                              tGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATgg                                        <  2:975744/139‑1 (MQ=255)
                                                                                              tGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATgg                                        <  1:445112/139‑1 (MQ=255)
                                                                                               gTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATGGa                                       <  1:1469675/139‑1 (MQ=255)
                                                                                                ttCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTg                                                                                                                 <  2:250658/64‑1 (MQ=255)
                                                                                                ttCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCAACCATCGCAATATCGTCATTg                                                                                                                 >  1:250658/1‑64 (MQ=255)
                                                                                                          gCGAGGTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATGGATGAGATGgcgc                            <  2:317638/139‑1 (MQ=255)
                                                                                                               gTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCt                                                                              <  2:264173/84‑1 (MQ=255)
                                                                                                               gTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCt                                                                              >  1:264173/1‑84 (MQ=255)
                                                                                                               gTTTCGCTCGCGGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATGGATGAGATGGCGCGTGGa                       >  2:432906/1‑139 (MQ=255)
                                                                                                                         cgGCACAGCAGAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATGGATGAGATGGCGCGTGGATTCATGCGTa             <  1:432906/139‑1 (MQ=255)
                                                                                                                              cagcagAAGTGCACCCATCGCAATATCGTCATTGGTGCAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATGGATGAGATGGCGCGTGGATTCATGCGTAATgggg       <  2:178916/139‑1 (MQ=255)
                                                                                                                               agcagAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATGGATGAGATGGCGCGTGGATTCATGCGTAATgggg       <  1:1567806/139‑1 (MQ=255)
                                                                                                                                gcagAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATGGATGAGATGGCGCGTGGATTCATGCGTAATGGGGa      <  1:849287/139‑1 (MQ=255)
                                                                                                                                gcagAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATGGATGAGATGGCGCGTGGATTCATGCGTAATGGGGa      <  1:1496690/139‑1 (MQ=255)
                                                                                                                                gcagAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATGGATGAGATGGCACGTGGATTCATGCGTAATGGGGa      >  2:436600/1‑139 (MQ=255)
                                                                                                                                 cagAAGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATGGATGAGATGGCGCGTGGATTCATGCGTAATGGGGaa     <  2:1366306/139‑1 (MQ=255)
                                                                                                                                    aaGTGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATGGATGAGATGGCGCGTGGATTCATGCGTAATGGGGAAAGa  <  2:446763/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
CCGACCAATATCAAAACGGGGAGTAATTACGCTGGCGAGGCTTGGGATCATTTGCCTGCCAATCTCCAGCCCATGAAAACCCGCAATGGAAATTTGTTCTGGCACGGCGAGGTTTCGCTCGCGGCACAGCAGAAGCGCACCCATCGCAATATCGTCATTGGTGCAAAATACACCGTCTAAATCCGGGTTTGCGCTAAGGGCATCACGCATTAACTGCATTCCTAAATGAATGGATGAGATGGCGCGTGGATTCATGCGTAATGGGGAAAGA  >  NC_000913/4490240‑4490510

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: