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,1500CT92.2% 438.0 / 21.4 128K90K (AAG→AAAfkpApeptidyl‑prolyl cis‑trans isomerase FkpA
Reads supporting (aligned to +/- strand):  ref base C (6/4);  new base T (58/60);  total (64/64)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

Only 100 of 133 total aligned reads displayed.

CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTAT  >  NC_000913/3477104‑3477196
                                              |                                              
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                              <  1:1252954/49‑1 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                              <  2:3985325/49‑1 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                              <  1:2799167/49‑1 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                             >  1:1477401/1‑50 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                            >  2:3799351/1‑49 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                            >  2:1764187/1‑49 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                            >  2:144181/1‑49 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                            >  2:3128078/1‑49 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATc                                           <  1:3621187/50‑1 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATc                                           <  1:721926/50‑1 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                          <  1:3851285/50‑1 (MQ=255)
    tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                         >  1:2209501/1‑50 (MQ=255)
     gCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                         >  1:4496691/1‑49 (MQ=255)
      cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                        <  1:4242923/49‑1 (MQ=255)
       ttGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                       >  1:5342234/1‑49 (MQ=255)
       ttGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                       >  1:818661/1‑49 (MQ=255)
        tGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                                     >  1:4772498/1‑50 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGGGTTTGCTCTTATCAGCaaagg                                   <  2:2165488/50‑3 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                    <  1:987577/49‑1 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                    <  2:1173602/49‑1 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                    <  1:3203426/49‑1 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                   >  2:1383889/1‑50 (MQ=255)
           agagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                   >  1:1235617/1‑49 (MQ=255)
            gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCa                                 <  1:1825680/50‑1 (MQ=255)
            gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCa                                 >  1:3647852/1‑50 (MQ=255)
            gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCa                                 >  2:465104/1‑50 (MQ=255)
             agTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                                <  1:395190/50‑1 (MQ=255)
             agTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                                <  1:2953494/50‑1 (MQ=255)
             agTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCAt                                >  1:2117009/1‑50 (MQ=255)
              gTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                               <  2:1324136/50‑1 (MQ=255)
              gTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                               >  2:3576417/1‑50 (MQ=255)
               tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcc                              <  1:4894736/50‑1 (MQ=255)
                cTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                             >  2:4923862/1‑50 (MQ=255)
                cTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                             >  2:3050658/1‑50 (MQ=255)
                 tGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                             >  1:2466060/1‑49 (MQ=255)
                  gTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                            >  1:939773/1‑49 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                             <  2:4869963/47‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                             <  1:4217427/47‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                             <  1:1452041/47‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                             <  2:2849634/47‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                            <  1:3860338/48‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                            <  1:3346062/48‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                            <  1:1437067/48‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                            <  1:4629749/48‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                          <  1:328044/50‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGa                           <  1:373481/49‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGa                           <  1:1503664/49‑1 (MQ=255)
                    tCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                          >  2:4871255/1‑49 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       >  1:284706/1‑50 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       >  2:2239106/1‑50 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAaca                        <  1:492107/49‑1 (MQ=255)
                       aTCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       >  2:4138833/1‑49 (MQ=255)
                       aTCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAacac                       >  1:4589508/1‑49 (MQ=255)
                        tctcTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       >  1:1458004/1‑48 (MQ=255)
                        tctcTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       >  2:209427/1‑48 (MQ=255)
                        tctcTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       >  1:3813798/1‑48 (MQ=255)
                        tctcTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       >  2:3678759/1‑48 (MQ=255)
                        tctcTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       >  1:688198/1‑48 (MQ=255)
                        tctcTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       >  1:3449869/1‑48 (MQ=255)
                           ctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                   >  2:2797443/1‑49 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       <  2:1152134/44‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACAcc                      <  1:5296393/45‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCa                     <  1:218049/46‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                    <  1:209358/47‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGa                 >  1:1504714/1‑50 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGa                 >  2:3276421/1‑50 (MQ=255)
                            ttGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGc                   <  1:4651084/48‑1 (MQ=255)
                             tGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt                >  2:1237891/1‑50 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCg                  <  1:977115/47‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCg                  <  2:1466729/47‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGa                 <  1:4725618/48‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt                >  1:61714/1‑49 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt                >  2:4698209/1‑49 (MQ=255)
                              ggTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGa                 <  1:1216369/48‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATc               <  2:3303/49‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa              <  2:5196302/50‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa              <  1:1662162/50‑1 (MQ=255)
                                tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg             >  1:4172906/1‑50 (MQ=255)
                                 cGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc            >  2:486825/1‑50 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc            <  1:2915118/49‑1 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt           >  1:3231808/1‑50 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt           <  2:1786805/50‑1 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt           <  2:1307768/50‑1 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt           <  2:1570473/50‑1 (MQ=255)
                                   gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt           >  1:2782672/1‑49 (MQ=255)
                                   gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt           >  1:1170251/1‑49 (MQ=255)
                                    agagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGa         >  2:2941685/1‑50 (MQ=255)
                                    agagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGa         >  2:4057603/1‑50 (MQ=255)
                                    agagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGa         <  1:2226591/50‑1 (MQ=255)
                                     gagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGAt        <  1:1791077/50‑1 (MQ=255)
                                        tttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt     <  2:4567149/50‑1 (MQ=255)
                                        tttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt     <  1:2064730/50‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCttt    <  1:4346094/49‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCttt    <  2:4994086/49‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCttt    <  2:3382959/49‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCttt    <  2:2448597/49‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTa   >  1:4054204/1‑50 (MQ=255)
                                           gCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTa   <  1:3218233/49‑1 (MQ=255)
                                           gCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTAt  >  1:44617/1‑50 (MQ=255)
                                            ctctTATCAGCAAATGCATCCTGAACGCCAGCGATCAGCTGATCTTTAt  <  2:2241732/49‑1 (MQ=255)
                                              |                                              
CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTAT  >  NC_000913/3477104‑3477196

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