Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,422,430 A→G D82G (GAC→GGC)  ulaD → 3‑keto‑L‑gulonate‑6‑phosphate decarboxylase UlaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,422,4300AG100.0% 222.2 / NA 60D82G (GAC→GGC) ulaD3‑keto‑L‑gulonate‑6‑phosphate decarboxylase UlaD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/34);  total (26/34)

ATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGGCGC  >  NC_000913/4422294‑4422563
                                                                                                                                        |                                                                                                                                     
aTTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACgctggctg                                                                                                                                    <  1:122214/140‑1 (MQ=255)
  tCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTggg                                                                                                                                  <  2:31099/140‑1 (MQ=255)
    tgtgCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTg                                                                                                                                >  2:480904/1‑140 (MQ=255)
         gTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGt                                                                                                                           >  2:906442/1‑140 (MQ=255)
           gggCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTaa                                                                                                                         >  2:651949/1‑140 (MQ=255)
                aaGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGc                                                                                                                    >  2:188958/1‑140 (MQ=255)
                aaGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGc                                                                                                                    >  2:564270/1‑140 (MQ=255)
                 aGGCGTGCGTGCGGTTCGCGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCt                                                                                                                   <  2:537222/140‑1 (MQ=255)
                    cgtgcgtgCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCtgtg                                                                                                                <  1:612059/140‑1 (MQ=255)
                    cgtgcgtgCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCtgtg                                                                                                                <  2:781436/140‑1 (MQ=255)
                      tgcgtgcgGTTCGTGACCTGAACGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCg                                                                                                              <  1:651949/140‑1 (MQ=255)
                      tgcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCg                                                                                                              >  2:485612/1‑140 (MQ=255)
                      tgcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCg                                                                                                              <  2:464856/140‑1 (MQ=255)
                      tgcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCg                                                                                                              <  2:561486/140‑1 (MQ=255)
                      tgcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCg                                                                                                              <  2:191123/140‑1 (MQ=255)
                      tgcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCg                                                                                                              <  2:749980/140‑1 (MQ=255)
                      ggcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCg                                                                                                              <  2:166939/139‑1 (MQ=255)
                       gcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCtgtg                                                                                                                >  2:829720/1‑137 (MQ=255)
                       gcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCtgtg                                                                                                                <  1:829720/137‑1 (MQ=255)
                       gcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCgg                                                                                                             <  1:610093/140‑1 (MQ=255)
                       gcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCgg                                                                                                             >  2:153759/1‑140 (MQ=255)
                         gtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGat                                                                                                           <  1:383375/140‑1 (MQ=255)
                          tgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGata                                                                                                          >  1:598119/1‑140 (MQ=255)
                               ttCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAac                                                                                                     <  2:664244/140‑1 (MQ=255)
                                       cTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCaa                                                                                             <  1:872017/140‑1 (MQ=255)
                                                ctctACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAg                                                                                           <  2:67888/133‑1 (MQ=255)
                                                ctctACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAg                                                                                           >  1:67888/1‑133 (MQ=255)
                                                ctctACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCt                                                                                    >  1:452826/1‑140 (MQ=255)
                                                        gCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTgg                                                                            <  1:480904/140‑1 (MQ=255)
                                                          acaAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCa                                                                          >  1:491953/1‑140 (MQ=255)
                                                             aaaTCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCaaaa                                                                       <  2:452826/140‑1 (MQ=255)
                                                                tCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAagag                                                                    >  2:430066/1‑140 (MQ=255)
                                                                   tACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGttt                                                                 >  2:24581/1‑140 (MQ=255)
                                                                      tGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAAc                                                              >  2:521343/1‑140 (MQ=255)
                                                                        gCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACgg                                                            >  2:403617/1‑140 (MQ=255)
                                                                         cAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGc                                                           <  1:557691/140‑1 (MQ=255)
                                                                          aGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCg                                                          >  2:346345/1‑140 (MQ=255)
                                                                            aCGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGAc                                                        >  1:659569/1‑140 (MQ=255)
                                                                                         cGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAAc                                           <  1:504097/140‑1 (MQ=255)
                                                                                           aTGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTg                                         >  2:69791/1‑140 (MQ=255)
                                                                                            tGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGa                                        <  1:153759/140‑1 (MQ=255)
                                                                                             gCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGAc                                       <  2:659569/140‑1 (MQ=255)
                                                                                              cAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACt                                          <  1:310155/136‑1 (MQ=255)
                                                                                              cAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACt                                          >  2:310155/1‑136 (MQ=255)
                                                                                              cAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGAcc                                      <  2:599602/140‑1 (MQ=255)
                                                                                               aGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCg                                     <  1:346345/140‑1 (MQ=255)
                                                                                                      aTCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTg                              <  1:430066/140‑1 (MQ=255)
                                                                                                         cTTTCGCGTATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGAc                           <  1:485526/140‑1 (MQ=255)
                                                                                                                 tATGTGCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAAc                   <  1:470479/140‑1 (MQ=255)
                                                                                                                      gCTTCGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGgcg              <  1:521343/140‑1 (MQ=255)
                                                                                                                          cGAAGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAAc          >  2:65351/1‑140 (MQ=255)
                                                                                                                            aaGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAg        >  2:511315/1‑140 (MQ=255)
                                                                                                                            aaGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAg        >  2:250950/1‑140 (MQ=255)
                                                                                                                            aaGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAg        >  1:249685/1‑140 (MQ=255)
                                                                                                                            aaGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAg        >  1:452901/1‑140 (MQ=255)
                                                                                                                            aaGCCAACGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAg        >  1:707719/1‑140 (MQ=255)
                                                                                                                                 aaCGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGgcg   <  1:511315/140‑1 (MQ=255)
                                                                                                                                 aaCGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGgcg   <  1:564270/140‑1 (MQ=255)
                                                                                                                                 aaCGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGgcg   <  2:725324/140‑1 (MQ=255)
                                                                                                                                 aaCGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGgcg   <  1:188958/140‑1 (MQ=255)
                                                                                                                                  aCGCTGGCTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGgcgc  <  1:146405/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
ATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTGCAGATCGAACTGACCGGTTACTGGACCTGGGAACAGGCGCAACAGTGGCGC  >  NC_000913/4422294‑4422563

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