Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,960,578 2 bp→CG coding (1266‑1267/1485 nt) gppA ← guanosine pentaphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,960,5780GC100.0% 354.3 / NA 132L423V (CTC→GTC) gppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (60/72);  total (60/72)
*NC_0009133,960,5790CG99.2% 363.3 / ‑5.0 132R422R (CGG→CGCgppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (59/72);  minor base A (1/0);  total (60/72)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 132 total aligned reads displayed.

CGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAGCCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTA  >  NC_000913/3960509‑3960646
                                                                     ||                                                                   
cGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcgg                                                                >  1:1539593/1‑76 (MQ=255)
 gagaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcggg                                                               >  1:5365395/1‑75 (MQ=255)
 gagaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcggc                                                               <  1:307086/76‑1 (MQ=255)
 gagaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcggc                                                               <  1:3783056/76‑1 (MQ=255)
 gagaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcggc                                                               <  1:4702144/76‑1 (MQ=255)
  agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcggcg                                                              <  2:2761647/76‑1 (MQ=255)
    aTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGt                                                            >  2:3947458/1‑76 (MQ=255)
    aTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGt                                                            <  2:5000075/76‑1 (MQ=255)
    aTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGt                                                            <  1:671074/76‑1 (MQ=255)
     tCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTa                                                           >  1:513224/1‑76 (MQ=255)
      cGTCACGGCGACGGCTGTCAAAAATGATGTCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTAc                                                          <  1:5583435/76‑1 (MQ=255)
       gTCGCGTCGTCGGCTGGCAAAAATGCTGGCGAGGCGTAGTAAAGTGCCGAGTTGTTCTGCGACGCGCGGCGGTTCg                                                         <  2:2245476/76‑1 (MQ=255)
       gTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACg                                                         >  2:3867731/1‑76 (MQ=255)
       gTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACg                                                         >  2:1710550/1‑76 (MQ=255)
        tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACgg                                                        >  1:1578354/1‑76 (MQ=255)
         cACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGc                                                       >  2:805927/1‑76 (MQ=255)
         cACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGc                                                       <  1:319276/76‑1 (MQ=255)
          aCGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCa                                                      <  2:3013182/76‑1 (MQ=255)
           cGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGcct                                                     >  2:4503362/1‑74 (MQ=255)
           cGGCGACGGCGGGCAAAAATGATGGCCAGGCGTAGTAAACGGCACAGGTGTTCAGAGACGCGCGGCGGCACgggct                                                     >  2:4372951/1‑73 (MQ=255)
            ggCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtt                                                    >  1:2800957/1‑76 (MQ=255)
              cGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtttt                                                  <  1:4601176/76‑1 (MQ=255)
               gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTg                                                 >  1:3610882/1‑76 (MQ=255)
               gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTg                                                 >  2:937449/1‑76 (MQ=255)
               gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTg                                                 <  1:3305082/76‑1 (MQ=255)
               gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTg                                                 >  1:2317126/1‑76 (MQ=255)
                aCGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGc                                                >  1:90022/1‑76 (MQ=255)
                aCGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGc                                                <  1:5077018/76‑1 (MQ=255)
                aCGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGc                                                <  1:4412451/76‑1 (MQ=255)
                  ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTg                                              >  1:5649437/1‑76 (MQ=255)
                  ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTg                                              >  1:4903114/1‑76 (MQ=255)
                  ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTg                                              >  1:5821789/1‑76 (MQ=255)
                   gCTGGCAAAAATGAGGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                             <  2:5267732/76‑1 (MQ=255)
                    gtGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGAt                                            <  2:5036343/75‑1 (MQ=255)
                    cTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGAt                                            <  2:4475729/76‑1 (MQ=255)
                    cTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGAt                                            >  1:2072488/1‑76 (MQ=255)
                     gggCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATg                                           <  1:347259/75‑1 (MQ=255)
                      ggCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGc                                          >  1:7654/1‑76 (MQ=255)
                       gCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCa                                         >  1:2330830/1‑76 (MQ=255)
                       gCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCa                                         >  2:1907435/1‑76 (MQ=255)
                       gCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCa                                         <  2:513224/76‑1 (MQ=255)
                         aaaagTGATGTGCAGGCATAGTAAACGGCAGAGGTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGc                                       <  1:4372952/76‑1 (MQ=255)
                           aaaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGa                                     <  2:4353805/76‑1 (MQ=255)
                           aaaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGa                                     >  2:2339938/1‑76 (MQ=255)
                           aaaTGAGGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGa                                     <  2:5748751/76‑1 (MQ=255)
                            aaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCCTTTTGCCGCTGCAGCgct                                    >  2:588078/1‑74 (MQ=255)
                            aaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGAt                                    <  1:5466917/76‑1 (MQ=255)
                            aaTGATGGCCAGGAGAAGAAAACGGCAGAGGTGTTGCGCGACGCGCGGCCGCCGGGGATGTTGGTGATGCAGCgct                                    >  1:3222460/1‑74 (MQ=255)
                             atgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATg                                   >  2:2234423/1‑76 (MQ=255)
                             atgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATg                                   <  2:2152991/76‑1 (MQ=255)
                               gatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgag                                 <  2:5063302/76‑1 (MQ=255)
                               gatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgag                                 <  1:461051/76‑1 (MQ=255)
                                 tgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagag                               >  1:2888159/1‑76 (MQ=255)
                                 tgGCCAGGCGTAGGAAAGGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagag                               >  1:2594137/1‑76 (MQ=255)
                                  ggCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagaga                              >  1:481656/1‑76 (MQ=255)
                                  ggCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagaga                              <  2:880839/76‑1 (MQ=255)
                                   gcgaGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTTCGGCATTTTGCTGATGCAGCGGTGAGAGAt                             <  2:3718830/73‑1 (MQ=255)
                                    ccAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATc                            <  1:3155168/76‑1 (MQ=255)
                                    ccAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATc                            <  1:1967425/76‑1 (MQ=255)
                                    ccAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATc                            <  2:1421098/76‑1 (MQ=255)
                                     cAGGCGTGGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCg                           <  2:2929369/76‑1 (MQ=255)
                                     cAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGCGAGATCg                           >  1:4171072/1‑76 (MQ=255)
                                     cAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCg                           <  2:3556064/76‑1 (MQ=255)
                                     cAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGGACGGCATTTTGCTGATGCAGGGGTGAGAGATCg                           >  2:5633442/1‑76 (MQ=255)
                                      aGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGa                          >  1:563759/1‑76 (MQ=255)
                                      aGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGa                          >  1:1741278/1‑76 (MQ=255)
                                      aGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGa                          <  1:4338733/76‑1 (MQ=255)
                                      aGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCCATGATAGATCGa                          >  1:1781389/1‑76 (MQ=255)
                                       ggCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGAc                         <  1:2558986/76‑1 (MQ=255)
                                        gCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGAcc                        <  1:4184222/76‑1 (MQ=255)
                                         cGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCCTTTTGCTGATGCAGCGATGAGAGATCGACCg                       >  2:1477573/1‑76 (MQ=255)
                                          gtagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCgg                      <  1:1536238/76‑1 (MQ=255)
                                          gtagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCgg                      >  1:2679742/1‑76 (MQ=255)
                                          gtagtaAACGGCAGAGGTGTTCTGCGACGCGCGGAGGAACGGGAGTTTGCCGATGCAGCCGTGCGAGATGGATCgg                      >  1:3491763/1‑76 (MQ=255)
                                           tagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGa                     <  1:5259213/76‑1 (MQ=255)
                                             gtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGAtt                   >  1:1007872/1‑76 (MQ=255)
                                             gtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGAtt                   <  2:4227331/76‑1 (MQ=255)
                                               aaaCGGCAGAGTTGTTCTGCGGCGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAg                 <  2:549137/76‑1 (MQ=255)
                                                aaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGt                <  2:1007872/76‑1 (MQ=255)
                                                aaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGt                >  2:4717175/1‑76 (MQ=255)
                                                 aCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTc               >  1:491899/1‑76 (MQ=255)
                                                 aCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTc               <  1:3947459/76‑1 (MQ=255)
                                                 aCGGCAGAGTTGTTCTGCGACGCGCCGCGGCAAGGCACCTTGCTGATGCAGCGATGGGAGGTCGATCGGGATAGTc               >  1:209492/1‑76 (MQ=255)
                                                  cGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCt              <  2:2305579/76‑1 (MQ=255)
                                                   ggCAGAGTTGTTCTGCGACACGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTg             >  1:4095266/1‑76 (MQ=255)
                                                    gCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTgg            <  2:2399727/76‑1 (MQ=255)
                                                     cAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGt           <  2:3491763/76‑1 (MQ=255)
                                                      agagTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGtt          <  2:2631798/76‑1 (MQ=255)
                                                       gagTTGTTCTGCGACGCGCGGCGTTACGGCATTTTGCTGATGCAGCGATGAGAGATCGGCCCGCTTAGTCTgggcg         >  2:482071/1‑73 (MQ=255)
                                                       gagTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCCGCGATGAGAGATCGGTCGGATTAGTCTGGTTg         <  1:589075/76‑1 (MQ=255)
                                                       gagTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGCGATCGCCCGGATTAGTCTGGTTg         >  2:110424/1‑76 (MQ=255)
                                                       gagTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTg         <  2:2038296/76‑1 (MQ=255)
                                                        agTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTgg        >  2:4755612/1‑75 (MQ=255)
                                                         gttgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAg       >  2:4313861/1‑76 (MQ=255)
                                                          ttgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGc      <  1:5106050/76‑1 (MQ=255)
                                                           tgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCa     <  1:2241529/76‑1 (MQ=255)
                                                            gttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAg    <  1:2489041/76‑1 (MQ=255)
                                                             ttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGTGCAGt   >  2:3835414/1‑76 (MQ=255)
                                                              tCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTa  <  2:1848203/76‑1 (MQ=255)
                                                              tCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTa  <  1:4503363/76‑1 (MQ=255)
                                                                     ||                                                                   
CGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAGCCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTA  >  NC_000913/3960509‑3960646

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