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% 289.5 / NA 71K90K (AAG→AAAfkpApeptidyl‑prolyl cis‑trans isomerase FkpA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/51);  total (20/51)

CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTA  >  NC_000913/3477104‑3477195
                                              |                                             
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                             <  1:3741341/49‑1 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                             <  1:577721/49‑1 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                             <  1:800775/49‑1 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                            >  1:3978486/1‑50 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                            >  1:4429059/1‑50 (MQ=255)
 gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           <  1:3244829/50‑1 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATc                                          <  1:3711384/50‑1 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           >  1:2703635/1‑48 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           >  1:2237195/1‑48 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:2378704/50‑1 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:810856/50‑1 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:4176275/50‑1 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:1296659/50‑1 (MQ=255)
    tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                        >  1:707588/1‑50 (MQ=255)
    tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                        <  1:821455/50‑1 (MQ=255)
      cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                       <  1:359469/49‑1 (MQ=255)
      cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                      <  1:2159936/50‑1 (MQ=255)
       ttGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                      >  1:4263675/1‑49 (MQ=255)
        tGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                                    >  1:3692849/1‑50 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAtt                                  <  1:1755641/50‑2 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                   <  1:2381994/49‑1 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                   <  1:2410510/49‑1 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                  <  1:1106264/50‑1 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                  <  1:3729671/50‑1 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                  <  1:539879/50‑1 (MQ=255)
            gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGc                                 <  1:3322262/49‑1 (MQ=255)
            gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCa                                >  1:885874/1‑50 (MQ=255)
             agTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                               >  1:2153672/1‑50 (MQ=255)
             agTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                               >  1:3411902/1‑50 (MQ=255)
               tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcc                             <  1:719688/50‑1 (MQ=255)
               tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcc                             <  1:3239704/50‑1 (MQ=255)
                  gTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                           >  1:2768996/1‑49 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                           <  1:3217053/48‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                           <  1:4297687/48‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                           <  1:1839733/48‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         <  1:3041499/50‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         <  1:2696792/50‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGa                          <  1:2951516/49‑1 (MQ=255)
                     cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAac                        <  1:2330586/49‑1 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                      >  1:2542601/1‑50 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                      >  1:517473/1‑50 (MQ=255)
                        tctcTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                      >  1:1481965/1‑48 (MQ=255)
                          tctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCa                    <  1:861427/48‑1 (MQ=255)
                          tctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  >  1:777708/1‑50 (MQ=255)
                          tctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  >  1:3353518/1‑50 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCa                    <  1:2423158/46‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                   <  1:508417/47‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                   <  1:3410062/47‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                   <  1:1297003/47‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  <  1:2452785/48‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  <  1:2063722/48‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  <  1:3303386/48‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  <  1:2912245/48‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCg                 <  1:747565/47‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt               <  1:2457879/49‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt               >  1:2384959/1‑49 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATc              <  1:2386562/49‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATc              <  1:1867096/49‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATc              <  1:2311986/49‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa             <  1:439442/50‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa             <  1:3359046/50‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa             <  1:2706031/50‑1 (MQ=255)
                                tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg            <  1:606137/50‑1 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          >  1:1071188/1‑50 (MQ=255)
                                     gagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGAt       <  1:728234/50‑1 (MQ=255)
                                      agTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc      <  1:2121123/50‑1 (MQ=255)
                                      agTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc      >  1:988313/1‑50 (MQ=255)
                                       gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt     >  1:2519402/1‑50 (MQ=255)
                                        tttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt    <  1:2278897/50‑1 (MQ=255)
                                        tttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt    <  1:37807/50‑1 (MQ=255)
                                        tttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt    <  1:3206777/50‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCttt   <  1:2896893/49‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCttt   <  1:3324358/49‑1 (MQ=255)
                                           gCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTa  <  1:2129903/49‑1 (MQ=255)
                                              |                                             
CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTA  >  NC_000913/3477104‑3477195

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