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

CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTA  >  NC_000913/3477104‑3477195
                                              |                                             
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                             <  1:5133265/49‑1 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                            >  1:2245985/1‑50 (MQ=255)
 gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           <  1:3242399/50‑1 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           >  1:2662138/1‑49 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           >  1:3538179/1‑49 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           >  1:2191503/1‑49 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           >  1:3790147/1‑49 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATc                                          <  1:3070390/50‑1 (MQ=255)
  aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATc                                          <  1:3780793/50‑1 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:1375253/50‑1 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:3933185/50‑1 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:4968162/50‑1 (MQ=255)
    tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                        <  1:842631/50‑1 (MQ=255)
    tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                        >  1:1740606/1‑50 (MQ=255)
      cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                       >  1:2577185/1‑49 (MQ=255)
      cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                       <  1:2045677/49‑1 (MQ=255)
      cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                      <  1:1758238/50‑1 (MQ=255)
       ttGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                      >  1:1692415/1‑49 (MQ=255)
       ttGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                      >  1:4155502/1‑49 (MQ=255)
         gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaa                                     <  1:4741549/48‑1 (MQ=255)
         gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                   >  1:4133599/1‑50 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                   <  1:5051659/49‑1 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                   <  1:1570836/49‑1 (MQ=255)
          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                  >  1:572992/1‑50 (MQ=255)
           agagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                  >  1:1282400/1‑49 (MQ=255)
            gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGc                                 <  1:568437/49‑1 (MQ=255)
            gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCa                                <  1:1596903/50‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                            <  1:1405012/47‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                            <  1:4593825/47‑1 (MQ=255)
                    tCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         >  1:2480745/1‑49 (MQ=255)
                     cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAac                        <  1:576249/49‑1 (MQ=255)
                     cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAac                        <  1:4744847/49‑1 (MQ=255)
                     cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAac                        <  1:1161773/49‑1 (MQ=255)
                     cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAac                        <  1:2115913/49‑1 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                      >  1:3069795/1‑50 (MQ=255)
                      gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAaca                       <  1:3491391/49‑1 (MQ=255)
                         ctctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                   >  1:3746901/1‑50 (MQ=255)
                          tctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  >  1:1274919/1‑50 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACAcc                     <  1:2514376/45‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                   <  1:3833487/47‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  <  1:4969305/48‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  <  1:213824/48‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  <  1:500201/48‑1 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  <  1:2328432/48‑1 (MQ=255)
                             tGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt               >  1:1927114/1‑50 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCg                 <  1:1956783/47‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGa                <  1:3511637/48‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt               <  1:4047638/49‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt               >  1:2445486/1‑49 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATc              <  1:747086/49‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATc              <  1:4773904/49‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa             <  1:740697/50‑1 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa             <  1:3167082/50‑1 (MQ=255)
                                tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg            >  1:1136395/1‑50 (MQ=255)
                                tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg            <  1:1604850/50‑1 (MQ=255)
                                tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg            >  1:3471015/1‑50 (MQ=255)
                                 cGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc           >  1:2153930/1‑50 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          >  1:2032462/1‑50 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          >  1:840881/1‑50 (MQ=255)
                                   gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          >  1:210144/1‑49 (MQ=255)
                                   gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          >  1:4978911/1‑49 (MQ=255)
                                   gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTg         <  1:2976331/50‑1 (MQ=255)
                                    agagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGa        <  1:3088809/50‑1 (MQ=255)
                                     gagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGAt       >  1:3942189/1‑50 (MQ=255)
                                     gagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGAt       <  1:3459972/50‑1 (MQ=255)
                                      agTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTTATc      >  1:1389727/1‑50 (MQ=255)
                                      agTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc      <  1:668191/50‑1 (MQ=255)
                                       gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt     <  1:2046666/50‑1 (MQ=255)
                                       gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt     >  1:5006885/1‑50 (MQ=255)
                                        tttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt    <  1:1321050/50‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCttt   <  1:2332273/49‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTa  >  1:1460577/1‑50 (MQ=255)
                                              |                                             
CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTA  >  NC_000913/3477104‑3477195

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