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.5 / NA 69D212D (GAT→GACtrpStryptophan‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/34);  total (35/34)

CGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512957‑3513048
                                              |                                             
cgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            >  1:3215097/1‑50 (MQ=255)
  gAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:1886536/1‑50 (MQ=255)
   aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTcc                                         <  1:4357818/50‑1 (MQ=255)
   aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:3491783/1‑49 (MQ=255)
   aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:4115245/1‑49 (MQ=255)
   aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:5358675/1‑49 (MQ=255)
    gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:5892288/1‑48 (MQ=255)
    gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCca                                        <  1:4999423/50‑1 (MQ=255)
    gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCca                                        <  1:4493564/50‑1 (MQ=255)
    gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCca                                        <  1:1405966/50‑1 (MQ=255)
      cAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:4039496/1‑49 (MQ=255)
       aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagca                                     <  1:1613972/50‑1 (MQ=255)
       aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagca                                     <  1:5131353/50‑1 (MQ=255)
       aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagc                                      <  1:1722567/49‑1 (MQ=255)
       aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagc                                      <  1:1422897/49‑1 (MQ=255)
        gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    <  1:5499615/50‑1 (MQ=255)
        gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    >  1:5582273/1‑50 (MQ=255)
        gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    <  1:5220725/50‑1 (MQ=255)
         tGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAgg                                   <  1:2211804/50‑1 (MQ=255)
           aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:524521/50‑1 (MQ=255)
           aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:1505687/50‑1 (MQ=255)
           aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:3233617/50‑1 (MQ=255)
            ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:825034/1‑49 (MQ=255)
            ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:1811172/1‑49 (MQ=255)
              gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:306073/1‑49 (MQ=255)
              gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              <  1:4402759/50‑1 (MQ=255)
               cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              >  1:5096004/1‑49 (MQ=255)
               cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATa                             <  1:715846/50‑1 (MQ=255)
                aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            >  1:4064061/1‑50 (MQ=255)
                aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            >  1:5847155/1‑50 (MQ=255)
                 cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           >  1:5021121/1‑50 (MQ=255)
                 cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           <  1:2072275/50‑1 (MQ=255)
                  gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          <  1:21461/50‑1 (MQ=255)
                  gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          >  1:6011180/1‑50 (MQ=255)
                  gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          <  1:1666252/50‑1 (MQ=255)
                  gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          >  1:4741074/1‑50 (MQ=255)
                   tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         <  1:394944/50‑1 (MQ=255)
                   tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         <  1:5165654/50‑1 (MQ=255)
                   tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  1:3304249/1‑50 (MQ=255)
                     tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGtt                        <  1:3782505/49‑1 (MQ=255)
                     tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       >  1:3011381/1‑50 (MQ=255)
                     tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       >  1:4968885/1‑50 (MQ=255)
                      gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       <  1:4485502/49‑1 (MQ=255)
                      gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       <  1:545533/49‑1 (MQ=255)
                      gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      >  1:5031866/1‑50 (MQ=255)
                       aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:4767242/1‑50 (MQ=255)
                       aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:429740/1‑50 (MQ=255)
                       aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:3560267/1‑50 (MQ=255)
                        ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    <  1:5048831/50‑1 (MQ=255)
                        ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    >  1:3815228/1‑50 (MQ=255)
                        ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    <  1:113016/50‑1 (MQ=255)
                         tttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   >  1:4574839/1‑50 (MQ=255)
                             ttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:3833031/1‑49 (MQ=255)
                               cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAtt               <  1:665676/48‑1 (MQ=255)
                               cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              <  1:2663368/49‑1 (MQ=255)
                                aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            <  1:4992590/50‑1 (MQ=255)
                                 cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           >  1:4232564/1‑50 (MQ=255)
                                  tACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:1599725/1‑50 (MQ=255)
                                    ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTc         <  1:337143/49‑1 (MQ=255)
                                    ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCa        <  1:2433495/50‑1 (MQ=255)
                                     cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:5342873/1‑47 (MQ=255)
                                     cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:5035885/1‑47 (MQ=255)
                                     cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAg       <  1:1223722/50‑1 (MQ=255)
                                      gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:5459767/1‑50 (MQ=255)
                                      gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:2782636/1‑50 (MQ=255)
                                       aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  1:3353366/1‑50 (MQ=255)
                                          tCGGGTCTTCCAGCAGGCCGATTACGTTATTGCGATTATCGTCAGACtt   <  1:20702/49‑1 (MQ=38)
                                          tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:364516/1‑50 (MQ=255)
                                           cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  1:4838378/49‑1 (MQ=255)
                                              |                                             
CGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512957‑3513048

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: