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,1500CT99.3% 542.0 / ‑3.8 139K90K (AAG→AAAfkpApeptidyl‑prolyl cis‑trans isomerase FkpA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (59/79);  total (59/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

Only 100 of 143 total aligned reads displayed.

CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTA  >  NC_000913/3477104‑3477195
                                              |                                             
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                             <  1:5297504/49‑1 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                             <  1:1389118/49‑1 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                             <  1:8136577/49‑1 (MQ=255)
 gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                            <  1:5626690/49‑1 (MQ=255)
 gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                            <  1:127334/49‑1 (MQ=255)
 gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           <  1:1724992/50‑1 (MQ=255)
 gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           <  1:7721373/50‑1 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           >  1:4458265/1‑49 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           >  1:3288577/1‑49 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           >  1:7066061/1‑49 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:6360234/50‑1 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:3513516/50‑1 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:5013358/50‑1 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:429014/50‑1 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:4206436/50‑1 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTAGGAGAGTTTGCTTTTATCa                                         <  1:4532621/50‑1 (MQ=255)
    tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                        <  1:4774563/50‑1 (MQ=255)
    tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                        >  1:5096460/1‑50 (MQ=255)
      cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                       >  1:353069/1‑49 (MQ=255)
      cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                       >  1:6681774/1‑49 (MQ=255)
      cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                       <  1:2299292/49‑1 (MQ=255)
        tttAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                                    >  1:6551357/3‑50 (MQ=255)
        tGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                                    >  1:6780324/1‑50 (MQ=255)
        tGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                                    >  1:2545707/1‑50 (MQ=255)
         gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                                    <  1:3160648/49‑1 (MQ=255)
         gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                   >  1:2146259/1‑50 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                  <  1:159756/50‑1 (MQ=255)
           agagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                  >  1:2471009/1‑49 (MQ=255)
           agagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                  >  1:3340430/1‑49 (MQ=255)
            gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCa                                <  1:4666233/50‑1 (MQ=255)
            gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATCCa                                <  1:4666143/50‑1 (MQ=255)
             agTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                               <  1:4696430/50‑1 (MQ=255)
               tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcc                             <  1:6606464/50‑1 (MQ=255)
               tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcc                             <  1:7235368/50‑1 (MQ=255)
               tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcc                             <  1:1413384/50‑1 (MQ=255)
                cTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                            >  1:6205002/1‑50 (MQ=255)
                 tGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                            >  1:4045629/1‑49 (MQ=255)
                  gTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                           >  1:2314495/1‑49 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                            <  1:6476131/47‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                           <  1:2849786/48‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                           <  1:5222249/48‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         <  1:1909914/50‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         <  1:6613860/50‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         <  1:7939189/50‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGa                          <  1:2077527/49‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGa                          <  1:4870118/49‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGa                          <  1:6847891/49‑1 (MQ=255)
                   ttCGATCTCTTAGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                            <  1:815043/47‑1 (MQ=39)
                    tCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         >  1:1229299/1‑49 (MQ=255)
                    tCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         >  1:326503/1‑49 (MQ=255)
                    tCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         >  1:477264/1‑49 (MQ=255)
                     cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAaca                       >  1:5496444/1‑50 (MQ=255)
                     cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAaca                       >  1:2626096/1‑50 (MQ=255)
                     cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAac                        <  1:5204754/49‑1 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                      >  1:7151160/1‑50 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                      >  1:7701078/1‑50 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                      >  1:3147964/1‑50 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAaca                       <  1:208812/49‑1 (MQ=255)
                         ctctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                   >  1:3963862/1‑50 (MQ=255)
                           ctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  >  1:2357909/1‑49 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCa                    <  1:7624941/46‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCa                    <  1:4372818/46‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCa                    <  1:2277934/46‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCa                    <  1:2316891/46‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                   <  1:638685/47‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  <  1:3461566/48‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  <  1:3573814/48‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGa                >  1:2361635/1‑50 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGa                >  1:4215043/1‑50 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt               <  1:3613766/49‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt               >  1:1645717/1‑49 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt               <  1:1125718/49‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa             <  1:3616443/50‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa             <  1:7734406/50‑1 (MQ=255)
                               gTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATc              <  1:6376910/49‑1 (MQ=255)
                                tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg            <  1:1724813/50‑1 (MQ=255)
                                tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg            <  1:2973041/50‑1 (MQ=255)
                                 cGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc           >  1:1105051/1‑50 (MQ=255)
                                 cGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc           >  1:1818815/1‑50 (MQ=255)
                                 cGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc           >  1:5053029/1‑50 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc           <  1:6486604/49‑1 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc           <  1:7906781/49‑1 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc           <  1:7569286/49‑1 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc           <  1:3106313/49‑1 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          >  1:8317108/1‑50 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          >  1:5892531/1‑50 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          <  1:3127817/50‑1 (MQ=255)
                                   gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          >  1:4388695/1‑49 (MQ=255)
                                   gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTg         <  1:2716533/50‑1 (MQ=255)
                                    agagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGa        >  1:7383859/1‑50 (MQ=255)
                                      agTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc      >  1:955561/1‑50 (MQ=255)
                                      agTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc      <  1:6569313/50‑1 (MQ=255)
                                      agTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc      >  1:551571/1‑50 (MQ=255)
                                       gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt     >  1:1027518/1‑50 (MQ=255)
                                        tttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt    <  1:3473227/50‑1 (MQ=255)
                                        tttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt    >  1:6591322/1‑50 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCttt   <  1:6317898/49‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTa  >  1:332917/1‑50 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTa  >  1:1635553/1‑50 (MQ=255)
                                           gCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTa  <  1:3261188/49‑1 (MQ=255)
                                              |                                             
CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTA  >  NC_000913/3477104‑3477195

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