Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,093,223 A→G V72V (GTT→GTC rhaM ← L‑rhamnose mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,093,2230AG100.0% 201.8 / NA 52V72V (GTT→GTCrhaML‑rhamnose mutarotase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/29);  total (23/29)

GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATC  >  NC_000913/4093177‑4093268
                                              |                                             
ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGAc                                             <  1:91006/49‑1 (MQ=255)
  ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                           <  1:2387550/49‑1 (MQ=255)
    aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATc                                          <  1:1398342/48‑1 (MQ=255)
     gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                         <  1:163973/48‑1 (MQ=255)
     gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                         <  1:2610342/48‑1 (MQ=255)
     gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                         <  1:1471183/48‑1 (MQ=255)
      cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                       <  1:2558724/49‑1 (MQ=255)
      cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                       <  1:766984/49‑1 (MQ=255)
       gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                      <  1:1666507/49‑1 (MQ=255)
        ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                     >  1:633140/1‑49 (MQ=255)
        ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                     <  1:2623544/49‑1 (MQ=255)
         gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                    >  1:2845751/1‑49 (MQ=255)
         gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                    <  1:1876742/49‑1 (MQ=255)
          cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                   <  1:2846127/49‑1 (MQ=255)
          cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                   <  1:2518420/49‑1 (MQ=255)
          cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                   <  1:1183469/49‑1 (MQ=255)
            tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                  >  1:665463/1‑48 (MQ=255)
            tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                 <  1:1808533/49‑1 (MQ=255)
             aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                 >  1:2356238/1‑48 (MQ=255)
             aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                 >  1:2132370/1‑48 (MQ=255)
              aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCa                                <  1:5578/48‑1 (MQ=255)
                 tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                            >  1:1501984/1‑49 (MQ=255)
                 tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                            >  1:1333378/1‑49 (MQ=255)
                 tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                            >  1:939249/1‑49 (MQ=255)
                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                          >  1:2610114/1‑49 (MQ=255)
                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                          <  1:969635/49‑1 (MQ=255)
                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                          >  1:1746499/1‑49 (MQ=255)
                    gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                         <  1:1960256/49‑1 (MQ=255)
                    gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                         <  1:876815/49‑1 (MQ=255)
                    gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                         >  1:2764310/1‑49 (MQ=255)
                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTcc                      <  1:2698602/49‑1 (MQ=255)
                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                       >  1:2766857/1‑48 (MQ=255)
                        tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCa                     <  1:789653/49‑1 (MQ=255)
                         atatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                    >  1:1528408/1‑49 (MQ=255)
                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                   >  1:1593963/1‑49 (MQ=255)
                           atTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                  >  1:1558422/1‑49 (MQ=255)
                              tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc               >  1:1167187/1‑49 (MQ=255)
                              tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                <  1:10842/48‑1 (MQ=255)
                               ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc               <  1:1113456/48‑1 (MQ=255)
                               ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc               <  1:2537366/48‑1 (MQ=255)
                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt             <  1:437467/49‑1 (MQ=255)
                                 accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc            <  1:1480963/49‑1 (MQ=255)
                                 accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt             >  1:1195591/1‑48 (MQ=255)
                                   caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt             >  1:1827425/1‑46 (MQ=255)
                                   caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt             >  1:1441363/1‑46 (MQ=255)
                                   caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAg          <  1:2754636/49‑1 (MQ=255)
                                      cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt       >  1:1487088/1‑49 (MQ=255)
                                      cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt       <  1:1517631/49‑1 (MQ=255)
                                       gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt       >  1:1606169/1‑48 (MQ=255)
                                        ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc      >  1:1395407/1‑48 (MQ=255)
                                        ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc      >  1:399115/1‑48 (MQ=255)
                                           gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc  <  1:960614/49‑1 (MQ=255)
                                              |                                             
GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATC  >  NC_000913/4093177‑4093268

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