Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,264,376 T→C I221I (ATT→ATC hemA → glutamyl tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,264,3760TC100.0% 148.9 / NA 43I221I (ATT→ATChemAglutamyl tRNA reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/23);  total (20/23)

CGCGGGCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCGCCAG  >  NC_000913/1264280‑1264468
                                                                                                |                                                                                            
cgcgGGCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTg                                                                                           >  2:378755/1‑100 (MQ=255)
     gCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCa                                                                                        >  2:250992/1‑98 (MQ=255)
     gCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGa                                                                                      <  1:339638/100‑1 (MQ=255)
      cGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGAt                                                                                     >  2:604239/1‑100 (MQ=255)
        aaaCTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGa                                                                                   <  1:695944/100‑1 (MQ=255)
           cTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGt                                                                                >  2:1120036/1‑100 (MQ=255)
            tATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTc                                                                               <  1:577275/100‑1 (MQ=255)
               cGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGgc                                                                            <  1:342443/100‑1 (MQ=255)
                   cTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCgg                                                                        <  2:974354/100‑1 (MQ=255)
                     ggtggCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGaa                                                                      >  1:341289/1‑100 (MQ=255)
                         gcgcGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGCAAGATGAAGTCGGCGCGGAAg                                                                     <  2:191974/97‑1 (MQ=255)
                            cgTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTg                                                               >  2:955267/1‑100 (MQ=255)
                             gTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGc                                                              <  2:635123/100‑1 (MQ=255)
                             gTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGc                                                              >  1:762411/1‑100 (MQ=255)
                                 tCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTg                                                          <  1:1097227/100‑1 (MQ=255)
                                    gcgcGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGt                                                       <  1:393847/100‑1 (MQ=255)
                                     cgcgAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTg                                                      >  1:182614/1‑100 (MQ=255)
                                        gAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGata                                                   >  2:576655/1‑100 (MQ=255)
                                           cacaAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATcg                                                <  2:1243310/100‑1 (MQ=255)
                                             naaaGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATcgac                                              >  1:338541/2‑100 (MQ=255)
                                             caAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAACCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATcgac                                              <  1:64999/100‑1 (MQ=255)
                                              aaaGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATcgac                                              >  2:243456/1‑99 (MQ=255)
                                                    cAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTc                                       <  1:1024208/100‑1 (MQ=255)
                                                           ccgatctCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGa                                <  1:593600/94‑1 (MQ=255)
                                                            gatTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGaa                               >  1:870479/1‑100 (MQ=255)
                                                                aTCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCg                           <  2:1017450/100‑1 (MQ=255)
                                                                  cGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAgatcgg                         >  2:593596/1‑95 (MQ=255)
                                                                   gCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGata                        <  1:495157/100‑1 (MQ=255)
                                                                    ccAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCAGCGCGAAGCCGatat                       >  1:1286270/1‑100 (MQ=255)
                                                                     cAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATatc                      <  2:240579/100‑1 (MQ=255)
                                                                        ccGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATatcatc                   <  2:174880/100‑1 (MQ=255)
                                                                           cACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATatcatcatn                <  2:227483/100‑2 (MQ=255)
                                                                             cTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAg              >  1:1287865/1‑100 (MQ=255)
                                                                             cTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAg              >  1:616120/1‑100 (MQ=255)
                                                                                gcgAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTc           <  1:244202/100‑1 (MQ=255)
                                                                                 cgAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTcc          >  1:908823/1‑100 (MQ=255)
                                                                                   aaCGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCa         <  2:979585/99‑1 (MQ=255)
                                                                                   aaCGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCAc        <  2:1267407/100‑1 (MQ=255)
                                                                                     cGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCg      <  2:613591/100‑1 (MQ=255)
                                                                                       tGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCGcc    >  2:166191/1‑100 (MQ=255)
                                                                                       tGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCGcc    <  1:811800/100‑1 (MQ=255)
                                                                                       tGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCGc     >  2:949490/1‑99 (MQ=255)
                                                                                         cccAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCGCCAg  >  1:174776/1‑100 (MQ=255)
                                                                                                |                                                                                            
CGCGGGCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCGCCAG  >  NC_000913/1264280‑1264468

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