Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,513,003 A→G D212D (GAT→GAC trpS ← tryptophan‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,513,0030AG100.0% 277.3 / NA 68D212D (GAT→GACtrpStryptophan‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/33);  total (35/33)

TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048
                                               |                                             
tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             >  1:3035379/1‑50 (MQ=255)
tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             <  1:5381797/50‑1 (MQ=255)
tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             >  1:5870359/1‑50 (MQ=255)
  ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            <  1:3811153/49‑1 (MQ=255)
  ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTGCGGGTCtt                                           >  1:2918975/1‑50 (MQ=255)
    aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTcc                                         <  1:661332/50‑1 (MQ=255)
    aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:4223133/1‑49 (MQ=255)
     gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCca                                        <  1:5576613/50‑1 (MQ=255)
      tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:73329/1‑50 (MQ=255)
        aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagc                                      <  1:2477608/49‑1 (MQ=255)
        aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:2741954/1‑48 (MQ=255)
         gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    <  1:1906663/50‑1 (MQ=255)
         gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    >  1:1533022/1‑50 (MQ=255)
          tGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAgg                                   <  1:381340/50‑1 (MQ=255)
          tGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAgg                                   >  1:2565016/1‑50 (MQ=255)
           gACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  >  1:6170545/1‑50 (MQ=255)
           gACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  >  1:6646478/1‑50 (MQ=255)
            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:3293334/1‑50 (MQ=255)
            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:3579618/50‑1 (MQ=255)
            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:2332162/50‑1 (MQ=255)
            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:1896541/50‑1 (MQ=255)
            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:1593237/1‑50 (MQ=255)
            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:4533817/50‑1 (MQ=255)
            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  <  1:1857198/49‑1 (MQ=255)
            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  <  1:359467/49‑1 (MQ=255)
             ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:3777120/1‑49 (MQ=255)
               gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:4685507/1‑49 (MQ=255)
                cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATa                             <  1:815131/50‑1 (MQ=255)
                 aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            >  1:6497570/1‑50 (MQ=255)
                 aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            >  1:4952474/1‑50 (MQ=255)
                 aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            >  1:1971403/1‑50 (MQ=255)
                  cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           <  1:5505987/50‑1 (MQ=255)
                  cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           <  1:2797038/50‑1 (MQ=255)
                  cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           <  1:4328679/50‑1 (MQ=255)
                   gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          >  1:2526068/1‑50 (MQ=255)
                   gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          <  1:6596419/50‑1 (MQ=255)
                      tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       >  1:3294066/1‑50 (MQ=255)
                      tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       >  1:3004487/1‑50 (MQ=255)
                       gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       <  1:906191/49‑1 (MQ=255)
                       gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       <  1:459402/49‑1 (MQ=255)
                       gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      >  1:3415912/1‑50 (MQ=255)
                        aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:868078/1‑50 (MQ=255)
                        aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     <  1:5681288/50‑1 (MQ=255)
                           ttcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   >  1:1899245/1‑49 (MQ=255)
                              ttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:589924/1‑49 (MQ=255)
                               tcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              >  1:6461814/1‑50 (MQ=255)
                                cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAtt               <  1:5600177/48‑1 (MQ=255)
                                cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAtt               <  1:5734497/48‑1 (MQ=255)
                                cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:3104264/47‑1 (MQ=255)
                                 aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            >  1:5678334/1‑50 (MQ=255)
                                 aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            >  1:3665339/1‑50 (MQ=255)
                                  cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           <  1:3204809/50‑1 (MQ=255)
                                   tACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:5173953/1‑50 (MQ=255)
                                     ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          <  1:1741305/48‑1 (MQ=255)
                                     ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCa        <  1:3353417/50‑1 (MQ=255)
                                     ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCa        <  1:1005437/50‑1 (MQ=255)
                                      cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:1783541/1‑47 (MQ=255)
                                       gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:4699303/1‑50 (MQ=255)
                                        aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  1:5814211/1‑50 (MQ=255)
                                        aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  1:105453/1‑50 (MQ=255)
                                         tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    <  1:4564961/50‑1 (MQ=255)
                                         tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    >  1:3636432/1‑50 (MQ=255)
                                           tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt   <  1:5131910/49‑1 (MQ=255)
                                           tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt   <  1:4368984/49‑1 (MQ=255)
                                            cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:23929/1‑49 (MQ=255)
                                            cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  1:3075712/49‑1 (MQ=255)
                                            cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  1:714048/49‑1 (MQ=255)
                                            cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:5393923/1‑49 (MQ=255)
                                               |                                             
TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048

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