Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 773,987 C→T P316S (CCG→TCG)  cydB → cytochrome d terminal oxidase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913773,9870CT100.0% 183.9 / NA 52P316S (CCG→TCG) cydBcytochrome d terminal oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/27);  total (25/27)

GGGTGTGGTTCTGCCGCTGCTGACCATCCTGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTT  >  NC_000913/773854‑774116
                                                                                                                                     |                                                                                                                                 
gggTGTGGTTCTGCCGCTGCTGACCATCCTGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTcc                                                                                                                              >  2:919251/1‑139 (MQ=255)
      ggTTCTGCCGCTGCTGACCATCCTGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCAcc                                                                                                                        <  2:416444/139‑1 (MQ=255)
           tGCCGCTGCTGACCATCCTGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCatgat                                                                                                                   >  2:1173093/1‑139 (MQ=255)
               gctgctGACCATCCGGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAAc                                                                                                               <  1:782841/139‑1 (MQ=255)
                ctgctgACCATCCTGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACg                                                                                                              <  1:1002248/139‑1 (MQ=255)
                           ccTGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACa                                                                                                   <  1:486069/139‑1 (MQ=255)
                           ccTGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACa                                                                                                   <  2:276583/139‑1 (MQ=255)
                             tGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGa                                                                                                     <  1:127552/135‑1 (MQ=255)
                             tGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGa                                                                                                     >  2:127552/1‑135 (MQ=255)
                               aCTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAAtgt                                                                                               <  1:424994/139‑1 (MQ=255)
                                             tAAAGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTcc                                                                                 >  1:431915/1‑139 (MQ=255)
                                               aaGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTccag                                                                               >  1:181743/1‑139 (MQ=255)
                                                  ccGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTccagcca                                                                            >  1:979610/1‑139 (MQ=255)
                                                      gTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTg                                                                        <  2:181743/139‑1 (MQ=255)
                                                       tGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGa                                                                       <  2:532465/139‑1 (MQ=255)
                                                                       cctccCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGAc                                                       <  1:584565/139‑1 (MQ=255)
                                                                          cccTGACGCTGGCCTGCGTCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTg                                                    >  2:436622/1‑139 (MQ=255)
                                                                             tGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGt                                                 >  2:47076/1‑139 (MQ=255)
                                                                                cGCTGGCCCGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATgtcgtc                                                                                                                               >  2:486890/1‑58 (MQ=255)
                                                                                cGCTGGCCCGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATgtcgtc                                                                                                                               <  1:486890/58‑1 (MQ=255)
                                                                                  cTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAACCAGCTGACGCTTAACGTCATGACCTGGGTTGCgg                                            >  2:778270/1‑139 (MQ=255)
                                                                                   tGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCggt                                           <  2:576474/139‑1 (MQ=255)
                                                                                        tGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTc                                      <  1:1014402/139‑1 (MQ=255)
                                                                                          catcatCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTg                                    >  2:694311/1‑139 (MQ=255)
                                                                                               tcCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACg                                                                                                              >  2:538995/1‑60 (MQ=255)
                                                                                               tcCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACg                                                                                                              <  1:538995/60‑1 (MQ=255)
                                                                                               tcCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACa                                                                                                   >  2:928581/1‑71 (MQ=255)
                                                                                               tcCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACa                                                                                                   <  1:928581/71‑1 (MQ=255)
                                                                                               tcCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAAtgt                                                                                               >  1:936055/1‑75 (MQ=255)
                                                                                               tcCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAAtgt                                                                                               <  2:936055/75‑1 (MQ=255)
                                                                                               tcCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTAcc                               >  2:63617/1‑139 (MQ=255)
                                                                                                 cTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGa                             <  1:47076/139‑1 (MQ=255)
                                                                                                    aCAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAg                                                                                                          <  1:615915/59‑1 (MQ=255)
                                                                                                    aCAGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAg                                                                                                          >  2:615915/1‑59 (MQ=255)
                                                                                                      aGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCAtt                        <  2:286325/139‑1 (MQ=255)
                                                                                                      aGCCGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCAtt                        <  1:1173093/139‑1 (MQ=255)
                                                                                                        ccGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCaa                                                                                                           <  1:1027879/54‑1 (MQ=255)
                                                                                                        ccGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCaa                                                                                                           >  2:1027879/1‑54 (MQ=255)
                                                                                                        ccGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCt                      >  2:967885/1‑139 (MQ=255)
                                                                                                         cGGTATCGCAATGTTCCCGTTTGTGATGTCGTTCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGt                                  >  2:718130/1‑126 (MQ=255)
                                                                                                         cGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCa                                                                                                                          >  1:708854/1‑38 (MQ=255)
                                                                                                         cGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCa                                                                                                                          <  2:708854/38‑1 (MQ=255)
                                                                                                         cGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTg                     >  1:200051/1‑139 (MQ=255)
                                                                                                         cGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTg                     >  2:281181/1‑139 (MQ=255)
                                                                                                         cGGTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATATGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGt                                  <  1:718130/126‑1 (MQ=255)
                                                                                                           gTATCGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGct                   <  2:979610/139‑1 (MQ=255)
                                                                                                               cGCAATGTTCCCGTTTGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTac               <  1:694311/139‑1 (MQ=255)
                                                                                                                            tttGTGATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGtt  >  2:999394/1‑139 (MQ=255)
                                                                                                                              tgtgATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGAt                                                                                          <  2:571672/49‑1 (MQ=255)
                                                                                                                              tgtgATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGAt                                                                                          >  1:571672/1‑49 (MQ=255)
                                                                                                                                tgATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGAcc                                                      <  2:958501/83‑1 (MQ=255)
                                                                                                                                tgATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGAcc                                                      >  1:958501/1‑83 (MQ=255)
                                                                                                                                tgATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCAtt                        <  2:832608/113‑1 (MQ=255)
                                                                                                                                tgATGTCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCAtt                        >  1:832608/1‑113 (MQ=255)
                                                                                                                                     |                                                                                                                                 
GGGTGTGGTTCTGCCGCTGCTGACCATCCTGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTT  >  NC_000913/773854‑774116

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