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%	285.9 / NA	71	D212D (GAT→GAC)	trpS	tryptophan‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (40/31);  total (40/31)

CGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512957‑3513048                                               |                                              cgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            >  1:2609022/1‑50 (MQ=255) cgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            >  1:5163553/1‑50 (MQ=255)  ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            <  1:1280898/49‑1 (MQ=255)    aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:958523/1‑49 (MQ=255)    aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:2847501/1‑49 (MQ=255)     gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:4508572/1‑48 (MQ=255)      tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       <  1:1166802/50‑1 (MQ=255)       cAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:4234176/1‑49 (MQ=255)       cAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:3981973/1‑49 (MQ=255)       cAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:2931732/1‑49 (MQ=255)        aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagc                                      <  1:4426901/49‑1 (MQ=255)        aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagc                                      <  1:3758733/49‑1 (MQ=255)        aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:695078/1‑48 (MQ=255)        aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:3127930/1‑48 (MQ=255)            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:2218948/1‑50 (MQ=255)            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:4347168/50‑1 (MQ=255)            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:4448730/50‑1 (MQ=255)            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:4462172/50‑1 (MQ=255)            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:5574666/50‑1 (MQ=255)             ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:5179988/1‑49 (MQ=255)              cGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:4583019/1‑50 (MQ=255)              cGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:3847613/1‑50 (MQ=255)               gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:4697924/1‑49 (MQ=255)               gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:1400832/1‑49 (MQ=255)               gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              <  1:500713/50‑1 (MQ=255)                cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATa                             <  1:4671046/50‑1 (MQ=255)                  cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           >  1:276149/1‑50 (MQ=255)                  cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           >  1:4400311/1‑50 (MQ=255)                      tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGtt                        <  1:3592768/49‑1 (MQ=255)                      tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGtt                        <  1:3667476/49‑1 (MQ=255)                       gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       <  1:3107215/49‑1 (MQ=255)                       gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       <  1:5493772/49‑1 (MQ=255)                       gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      >  1:4925188/1‑50 (MQ=255)                        aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     <  1:3708011/50‑1 (MQ=255)                        aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      <  1:1941904/49‑1 (MQ=255)                         ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    >  1:4870766/1‑50 (MQ=255)                         ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    <  1:969086/50‑1 (MQ=255)                         ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    >  1:2369994/1‑50 (MQ=255)                            tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCg                  <  1:1633691/49‑1 (MQ=255)                             cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:2198088/1‑50 (MQ=255)                             cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:57285/1‑50 (MQ=255)                             cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:3564221/1‑50 (MQ=255)                                cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAtt               <  1:303141/48‑1 (MQ=255)                                cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:4228225/47‑1 (MQ=255)                                cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:3825744/47‑1 (MQ=255)                                cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:1631172/47‑1 (MQ=255)                                cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              <  1:1133849/49‑1 (MQ=255)                                cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             >  1:866888/1‑50 (MQ=255)                                 aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             <  1:761084/49‑1 (MQ=255)                                 aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             <  1:5516220/49‑1 (MQ=255)                                 aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            <  1:3806783/50‑1 (MQ=255)                                  cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           >  1:3155815/1‑50 (MQ=255)                                   tACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:2103171/1‑50 (MQ=255)                                    aCCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:4443693/1‑49 (MQ=255)                                    aCCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:3222565/1‑49 (MQ=255)                                     ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          <  1:698254/48‑1 (MQ=255)                                     ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:5585145/1‑48 (MQ=255)                                     ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTc         <  1:1131519/49‑1 (MQ=255)                                     ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCa        <  1:5255901/50‑1 (MQ=255)                                     ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCa        <  1:5032686/50‑1 (MQ=255)                                      cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:5231408/1‑47 (MQ=255)                                      cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:4029173/1‑47 (MQ=255)                                       gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:2109635/1‑50 (MQ=255)                                       gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:3915345/1‑50 (MQ=255)                                        aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  1:5385911/1‑50 (MQ=255)                                        aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  1:3906526/1‑50 (MQ=255)                                         tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    <  1:248748/50‑1 (MQ=255)                                         tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    >  1:4438974/1‑50 (MQ=255)                                         tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    >  1:3915670/1‑50 (MQ=255)                                          ttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt   >  1:4118065/1‑50 (MQ=255)                                            cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:2769169/1‑49 (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: ‑