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_000913	3,513,003	0	A	G	100.0%	280.1 / NA	72	D212D (GAT→GAC)	trpS	tryptophan‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (40/32);  total (40/32)

GGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512958‑3513048                                              |                                              ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCtt                                           >  1:4367154/1‑50 (MQ=255)  gAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          <  1:4158821/50‑1 (MQ=255)  gAGTCAGTGACCGCACGTTTGATTTTCTCCACTACCGATTTCGGGTCTTc                                          <  1:3505587/50‑1 (MQ=255)   aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTcc                                         <  1:4003055/50‑1 (MQ=255)   aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTcc                                         <  1:2429250/50‑1 (MQ=255)   aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:1481423/1‑49 (MQ=255)   aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:1598582/1‑49 (MQ=255)   aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:1891734/1‑49 (MQ=255)    gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:2995510/1‑48 (MQ=255)    gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:424062/1‑48 (MQ=255)    gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCca                                        <  1:733526/50‑1 (MQ=255)     tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:2139785/1‑50 (MQ=255)       aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagca                                     <  1:440665/50‑1 (MQ=255)       aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagca                                     <  1:988969/50‑1 (MQ=255)       aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagc                                      <  1:4316341/49‑1 (MQ=255)        gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    >  1:2414303/1‑50 (MQ=255)        gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    <  1:5011356/50‑1 (MQ=255)        gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    <  1:5577912/50‑1 (MQ=255)        gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    <  1:3169858/50‑1 (MQ=255)           aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:5540787/50‑1 (MQ=255)            ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:2521914/1‑49 (MQ=255)             cGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:5845847/1‑50 (MQ=255)             cGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               <  1:5223415/50‑1 (MQ=255)             cGCACGTTTGATTTTCTTCACTACCGATTTCGGGGCTTCCAGCAGGCCGa                               >  1:4013495/1‑50 (MQ=39)               cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              >  1:1243781/1‑49 (MQ=255)               cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              >  1:6400020/1‑49 (MQ=255)                aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCTAGCAGGCCGATaa                            >  1:5033763/1‑50 (MQ=255)                aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            >  1:5396482/1‑50 (MQ=255)                aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            <  1:3034257/50‑1 (MQ=255)                aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            <  1:3001455/50‑1 (MQ=255)                aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            >  1:5108263/1‑50 (MQ=255)                aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            <  1:835289/50‑1 (MQ=255)                 cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           <  1:5677072/50‑1 (MQ=255)                  gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          >  1:653563/1‑50 (MQ=255)                   tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         <  1:4016294/50‑1 (MQ=255)                   tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  1:420558/1‑50 (MQ=255)                    ttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  1:6375920/1‑49 (MQ=255)                      gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      >  1:2406122/1‑50 (MQ=255)                       aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:3150748/1‑50 (MQ=255)                       aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      <  1:3329208/49‑1 (MQ=255)                        ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    >  1:993782/1‑50 (MQ=255)                        ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    >  1:6225284/1‑50 (MQ=255)                         tttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   >  1:5949241/1‑50 (MQ=255)                          ttcttcGCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCg                  <  1:5138315/50‑1 (MQ=255)                           tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCg                  <  1:1798152/49‑1 (MQ=255)                           tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCg                  <  1:4322892/49‑1 (MQ=255)                           tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCg                  <  1:3488319/49‑1 (MQ=255)                            cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:2604721/1‑50 (MQ=255)                            cctcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGa                 <  1:4398141/47‑1 (MQ=39)                             ttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:3751189/1‑49 (MQ=255)                               cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             >  1:6343787/1‑50 (MQ=255)                               cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             >  1:5886966/1‑50 (MQ=255)                                aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             <  1:319376/49‑1 (MQ=255)                                aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             <  1:605294/49‑1 (MQ=255)                                aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            >  1:45098/1‑50 (MQ=255)                                aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            <  1:1548797/50‑1 (MQ=255)                                aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            >  1:2122141/1‑50 (MQ=255)                                aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            >  1:2832999/1‑50 (MQ=255)                                 cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           >  1:2264463/1‑50 (MQ=255)                                  tACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:1308985/1‑50 (MQ=255)                                   aCCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:6038247/1‑49 (MQ=255)                                   aCCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:255117/1‑49 (MQ=255)                                   aCCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:930510/1‑49 (MQ=255)                                   aCCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:4571174/1‑49 (MQ=255)                                    ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTc         <  1:5026306/49‑1 (MQ=255)                                    ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTc         <  1:1787378/49‑1 (MQ=255)                                     cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:3981078/1‑47 (MQ=255)                                      gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      <  1:2611443/50‑1 (MQ=255)                                      gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:3159102/1‑50 (MQ=255)                                      gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:4624636/1‑50 (MQ=255)                                        tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    <  1:1815969/50‑1 (MQ=255)                                           cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  1:6321395/49‑1 (MQ=255)                                              |                                              GGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512958‑3513048

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