Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,603,050 A→C 100% E160A (GAG→GCG)  RCAP_rcc02427 → M23 family peptidase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,603,0500A→C100.0% 75.2 20E160A (GAG→GCG) RCAP_rcc02427M23 family peptidase
Reads supporting (aligned to +/- strand):  new base (5/15):  ref base (0/0):  total (5/15)

GCGGGCGGGGCTGCCGCCGCGGGCAGCGAGCGGATCGACATCACCTCGCTGGCGGGGGCCGCGATCGACCGCGCCGCCGTCACCCAGCCGCCCGCCGAGACGCCCGCGGCCCCCGCCACCGCGCCCCGTCCCGCCGCGCCAAGTGCCAGCCCGGTGGCGAAACCGGCCGCCCC.CGCCCAGACCGGCGCCGA  >  NC_014034/2602953‑2603143
                                                                                                 |                                                                                              
gggggggggGCTGCCGCCGCGGGCAGCGAGCGGATCGACATCCCCGCGCTGGGGGGGGCCGCGATCGACCGCGCCGCCGTCCCCCAGCCGCCCGCCGCGAc                                                                                             <  1:3854148/99‑1 (MQ=255)
   gggggggCGGCCGCCGCGGGCAGCGAGCGGATCGACATCCCCTCGCTGGCGGGGGCCGCGATCGACCGCGCCGCCGTCACCCAGCCGCCCGCCGCGACGcc                                                                                          <  1:4171490/98‑1 (MQ=255)
   ggcgggGCTGCCGCCGCGGGCAGCGAGCGGATCGACATCACCTCGCTGGCGGGGGCCGCGATCGACCGCGCCGCCGTCACCCAGCCGCCCGCCGCGACGcc                                                                                          <  2:3338504/101‑1 (MQ=255)
    gcgggGCTGCCGCCGCGGGCAGCGAGCGGATCGACATCACCTCGCTGGCGGGGGCCGCGATCGACCGCGCCGCCGTCACCCAGCCGCCCGCCGCGACGccc                                                                                         <  2:1207084/101‑1 (MQ=255)
    gcgggGCTGCCGCCGCGGGCAGCGAGCGGATCGACATCACCTCGCTGGCGGGGGCCGCGATCGACCGCGCCGCCGTCACCCAGCCGCCCGCCGCGACGccc                                                                                         <  1:2625638/101‑1 (MQ=255)
      gggcCTGCCGCCGCGGGCAGCGAGCGGATCGACATCACCTCGCTGGCGGGGGCCGCGATCGACCGCGCCGCCGTCACCCAGCCGCCCGCCGCGACGCCcgc                                                                                       <  2:3181054/101‑1 (MQ=255)
         gCTGCCCCCGCGGGCAGCGAGCGGATCGACATCACCTCGCTGTCGGGGGCCGCGTTCGACCGCGCCGCCGTCACCCAGCCCCCCGCCGCGACGCCCGCGGc                                                                                    <  2:2867819/101‑1 (MQ=255)
             ccgccgCGGGCAGCGAGCGGATCGACATCCCCTCGCTGGCGGGGGCCGCGATCGACCGCGCCGCCGTCACCCAGCCGCCCGCCGCGACGCCCGCGGccccc                                                                                <  1:1502124/101‑1 (MQ=255)
               gccgcGGGCAGCGAGCGGATCGACATCCCCTCGCTGGCGGGGGCCGCGATCGACCGCGCCGCCGTCACCCAGCCGCCCGCCGCGACGCCCGCGGCCCCCGc                                                                              <  1:5587443/101‑1 (MQ=255)
                 cgcgGGCAGCGAGCGGATCGACATCCCCTCGCTGGCGGGGGCCGCGATCGACCGCCCCCCCGTCCCCCAGCCGCCCGCCGCGACGCCCGCGGCCCCCGCCa                                                                            <  1:5779433/101‑1 (MQ=255)
                     ggCAGCGAGCGGATCGACATCCCCCCGCTGGCGGGGGCCGCGATCGACCGCGCCGCCGTCCCCCAGCCGCCCGCCGCGACGCCCGCGGCCCCCGCCACcgc                                                                        <  2:1535251/101‑1 (MQ=255)
                        agcgagcgGATCGACATCACCTCGCTGGCGGGGGGCGCGGGCGACCGCGCCGCCCCCACCCACCCGCCCGCCGCCACGCCCGCGGCCCCCCCCACCCCGcc                                                                     >  2:2532352/1‑101 (MQ=255)
                         gcgaccgGATCGACATCCCCTCGCTGGCGGGGGCCGCGATCGACCGCGCCGCCGTCACCCAGCCGCCCGCCGCGACGCCCGCGGCCCCCGCCACCGCGccc                                                                    <  1:2868384/101‑1 (MQ=255)
                              cgGATCGACATCACCTCGCTGGCGGGGGGCGCGGTCGACCGCGCCGCCGCCACCCCCCCCCCCCCCGCGGCCCCCCCGGCCCCCCCCCCCCCGcccccgcc                                                               >  2:3558863/1‑97 (MQ=255)
                               ggATCGACATCACCTCGCTGGCGGGGGCCGCGCCCCGCCGGGCCCCCCGGCCCCCCACGCCCGCCGCGCCGCCCCCCGGCCCCGCCCCCCCCccccccccc                                                              >  2:1989559/1‑96 (MQ=255)
                               ggATCGACATCACCTCGCTGGCGGGGGCCGCGATCGACCGCGCCCCCGGCACCCACCCCACCCCCCCCGCGCCCGCCGCCCCCCCCCCCCCCcccccgccc                                                              >  2:5737616/1‑96 (MQ=255)
                                                cgggCGGGGCCCGGTACCGGCCGCGCCGCCTGCACCCAGCCGCCCGCCGCGACGCCCGCGGCCCCCGCCACCGCGCCCCGTCCCGCCGCGCCAAGTGCCAg                                             <  1:5737911/99‑1 (MQ=255)
                                                           cgcgATCGACCGCGCCGCCGTCACCCAGCCGCCCGCCGCGACGCCCGCGGCCCCCGCCACCGCGCCCCGTCCCGCCGCGCCAAGTGCCAGCCCGGTGGCGa                                  >  1:1707814/1‑101 (MQ=255)
                                                             cccctccccGCCCCCCCGTCCCCCCCCCGCCCGCCGCGACGCCCGCGGCCCCCGCCACCGCGACCCGTCCCGCCGCGCCAAGTGCCAGCCCGGTGGCGaaa                                <  1:1653439/94‑1 (MQ=255)
                                                                  gACCGCGCCGCCGTCACCCAGCCGCCCGCCGCGACGCCCGCGGCCCCCGCCACCGCGACCCGTCCCGCCGCGCCAAGTGCCAGCCCGGTGGCGaaaagatc                           >  2:1653345/1‑96 (MQ=255)
                                                                   cccGCCCCGCCGCCCCCCACCCCCCCGCCGCGACGCCCGCGGCCCCCGCCACCGCGCCCCGTCCCGCCGCGCCAAGTGCCAGCCCGGTGGCGAAACCGgcc                          <  1:3559080/100‑1 (MQ=255)
                                                                             cgTCACCCAGCCGCCCGCCGCGACGCCCGCGGCCCCCGCCACCGCGCCCCGTCCCGCCGCGCCAAGTGCCAGCCCGGTGGCGAAACCGGccgcccc.cgccc               >  2:3645115/1‑101 (MQ=255)
                                                                                 aCCCAGCCGCCCGCCGCGACGCCCGCGGCCCCCGCCACCGCGCCCCGTCCCGCCGCGCCAAGTGCCAGCCCGGTGGCGAAACCGGCCGCCCC.CGCCCAGAc           >  2:2442023/1‑101 (MQ=255)
                                                                                   cccgcccccCGCCGCGACGCCCGCGGCCCCCGCCACCGCGCCCCGTCCCGCCGCGCCAAGTGCCAGCCCGGTGGCGAAACCGGCcgcccctcgccctcctc          <  2:1707715/98‑6 (MQ=255)
                                                                                          cccgccGCGACGCCCGCGGCCCCCGCCACCGCGCCCCGTCCCGCCGCGCCAAGTGCCAGCCCGGTGGCGAAACCGGCCGCCCC.CGCCCAGACCGGCgccga  >  1:1766838/1‑101 (MQ=255)
                                                                                                 |                                                                                              
GCGGGCGGGGCTGCCGCCGCGGGCAGCGAGCGGATCGACATCACCTCGCTGGCGGGGGCCGCGATCGACCGCGCCGCCGTCACCCAGCCGCCCGCCGAGACGCCCGCGGCCCCCGCCACCGCGCCCCGTCCCGCCGCGCCAAGTGCCAGCCCGGTGGCGAAACCGGCCGCCCC.CGCCCAGACCGGCGCCGA  >  NC_014034/2602953‑2603143

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 39 ≤  ATCG