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

CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTA  >  NC_000913/3477104‑3477195
                                              |                                             
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                             <  1:1311373/49‑1 (MQ=255)
cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                            >  1:1042793/1‑50 (MQ=255)
   aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:1916221/50‑1 (MQ=255)
    tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                        <  1:553507/50‑1 (MQ=255)
    tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                        <  1:650975/50‑1 (MQ=255)
      cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                       <  1:1791330/49‑1 (MQ=255)
      cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                      <  1:40462/50‑1 (MQ=255)
         gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaa                                     <  1:2192185/48‑1 (MQ=255)
         gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                   >  1:1812965/1‑50 (MQ=255)
           agagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                  >  1:543239/1‑49 (MQ=255)
               tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcc                             >  1:1401028/1‑50 (MQ=255)
               tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcc                             >  1:1898709/1‑50 (MQ=255)
               tCTGCTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcc                             <  1:1361467/50‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                            <  1:285607/47‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         <  1:1443955/50‑1 (MQ=255)
                   ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         <  1:1523998/50‑1 (MQ=255)
                     cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAac                        <  1:631835/49‑1 (MQ=255)
                       aTCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                      >  1:1207869/1‑49 (MQ=255)
                            ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                   <  1:1343005/47‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCg                 <  1:1134221/47‑1 (MQ=255)
                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt               >  1:1352192/1‑49 (MQ=255)
                               gTCGGAGAGTTTGCTTTTATCAGCAAATGCGTCCTGAACACCAGCGAt               >  1:917150/1‑48 (MQ=255)
                                 cGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc           >  1:1859668/1‑50 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc           <  1:451560/49‑1 (MQ=255)
                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          <  1:317465/50‑1 (MQ=255)
                                   gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTg         <  1:104179/50‑1 (MQ=255)
                                    agagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGa        >  1:2064867/1‑50 (MQ=255)
                                    agagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGa        >  1:267775/1‑50 (MQ=255)
                                    agagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGa        >  1:77885/1‑50 (MQ=255)
                                     gagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGAt       <  1:1500441/50‑1 (MQ=255)
                                      agTTTGCTTTTATCAGCAAATGCATCCGGAACACCAGCGATCAGCTGATc      <  1:1151961/50‑1 (MQ=38)
                                       gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt     <  1:1443013/50‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCttt   <  1:2014355/49‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCttt   <  1:1151314/49‑1 (MQ=255)
                                          tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTa  >  1:649530/1‑50 (MQ=255)
                                              |                                             
CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTA  >  NC_000913/3477104‑3477195

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