Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,477,150 C→T K90K (AAG→AAA fkpA ← peptidyl‑prolyl cis‑trans isomerase FkpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,477,1500CT100.0% 321.1 / NA 79K90K (AAG→AAAfkpApeptidyl‑prolyl cis‑trans isomerase FkpA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/51);  total (28/51)

CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTAT  >  NC_000913/3477104‑3477196
                                              |                                              
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                              <  1:3540544/49‑1 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                              <  1:6359990/49‑1 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                             >  1:2207630/1‑50 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                             >  1:1867039/1‑50 (MQ=255)
 gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                            <  1:3404869/50‑1 (MQ=255)
 gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                            >  1:6057201/1‑50 (MQ=255)
 gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                            >  1:2164847/1‑50 (MQ=255)
    tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                         <  1:6400718/50‑1 (MQ=255)
    tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                         <  1:3566998/50‑1 (MQ=255)
      cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                        >  1:2090472/1‑49 (MQ=255)
      cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                        <  1:3413036/49‑1 (MQ=255)
       ttGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                       >  1:4601084/1‑49 (MQ=255)
        tGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                                     >  1:5678470/1‑50 (MQ=255)
         gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                                     <  1:935344/49‑1 (MQ=255)
         gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaa                                      <  1:285829/48‑1 (MQ=255)
         gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaa                                      <  1:1053761/48‑1 (MQ=255)
         gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                    >  1:548119/1‑50 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                    <  1:3963433/49‑1 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                    <  1:6244231/49‑1 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                    <  1:5295861/49‑1 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                   >  1:2270697/1‑50 (MQ=255)
            gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGc                                  <  1:1116924/49‑1 (MQ=255)
            gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCa                                 >  1:3194565/1‑50 (MQ=255)
             agTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                                <  1:6388540/50‑1 (MQ=255)
             agTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                                <  1:2805013/50‑1 (MQ=255)
              gTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                               <  1:2189149/50‑1 (MQ=255)
               tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcc                              <  1:706017/50‑1 (MQ=255)
                  gTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                            >  1:3436957/1‑49 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                             <  1:309926/47‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                            <  1:5838686/48‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                            <  1:5685796/48‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                            <  1:2015722/48‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                          <  1:2049195/50‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGa                           <  1:3312256/49‑1 (MQ=255)
                     cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAac                         <  1:5447223/49‑1 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAaca                        <  1:6382508/49‑1 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAaca                        <  1:1857473/49‑1 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAaca                        <  1:2194717/49‑1 (MQ=255)
                       aTCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       >  1:3712973/1‑49 (MQ=255)
                       aTCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       >  1:4204254/1‑49 (MQ=255)
                          tctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                   >  1:4638395/1‑50 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       <  1:3700549/44‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       <  1:4851244/44‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                    <  1:4157321/47‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                    <  1:756555/47‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                    <  1:4656362/47‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                   <  1:5976534/48‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                   <  1:5886782/48‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                   <  1:1753900/48‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                   <  1:4756524/48‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAACGCATCCTGCACACCAGCGa                 >  1:3379665/1‑50 (MQ=255)
                             tGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt                >  1:4835086/1‑50 (MQ=255)
                             tGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt                >  1:3708296/1‑50 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCg                  <  1:1780329/47‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCg                  <  1:2209673/47‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCg                  <  1:5791072/47‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa              <  1:4514367/50‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa              <  1:6034640/50‑1 (MQ=255)
                                tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg             >  1:1672413/1‑50 (MQ=255)
                                tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg             >  1:155230/1‑50 (MQ=255)
                                tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg             <  1:2380936/50‑1 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc            <  1:4500078/49‑1 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc            <  1:3824435/49‑1 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt           >  1:94735/1‑50 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt           <  1:3361760/50‑1 (MQ=255)
                                   gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt           >  1:2520311/1‑49 (MQ=255)
                                   gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTg          <  1:5031459/50‑1 (MQ=255)
                                    agagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGa         <  1:1125250/50‑1 (MQ=255)
                                    agagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGa         >  1:868268/1‑50 (MQ=255)
                                      agTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc       >  1:5488210/1‑50 (MQ=255)
                                      agTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc       >  1:6265394/1‑50 (MQ=255)
                                      agTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc       >  1:3050958/1‑50 (MQ=255)
                                       gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt      >  1:3948113/1‑50 (MQ=255)
                                       gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt      <  1:3840580/50‑1 (MQ=255)
                                       gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt      <  1:6404888/50‑1 (MQ=255)
                                       gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt      <  1:2862042/50‑1 (MQ=255)
                                       gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt      <  1:798077/50‑1 (MQ=255)
                                        tttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt     <  1:1033492/50‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCttt    <  1:4934012/49‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTa   >  1:4215987/1‑50 (MQ=255)
                                           gCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTAt  >  1:5402415/1‑50 (MQ=255)
                                              |                                              
CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTAT  >  NC_000913/3477104‑3477196

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