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

TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048                                                |                                              tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             >  1:3583912/1‑50 (MQ=255)  cgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            >  1:2427199/1‑50 (MQ=255)    gAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          <  1:2489952/50‑1 (MQ=255)     aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTcc                                         <  1:2890403/50‑1 (MQ=255)     aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTcc                                         <  1:7405513/50‑1 (MQ=255)      gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:5182510/1‑48 (MQ=255)      gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCca                                        <  1:5073475/50‑1 (MQ=255)       tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:3810130/1‑50 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagca                                     <  1:4901435/50‑1 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagca                                     <  1:1909082/50‑1 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagca                                     <  1:6299830/50‑1 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagc                                      <  1:2122780/49‑1 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:6595648/1‑48 (MQ=255)          gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    <  1:2004703/50‑1 (MQ=255)          gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    <  1:1921716/50‑1 (MQ=255)          gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    <  1:3377018/50‑1 (MQ=255)          gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    >  1:2873844/1‑50 (MQ=255)           tGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAgg                                   <  1:7271633/50‑1 (MQ=255)           tGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAgg                                   <  1:234603/50‑1 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:3684171/1‑50 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:6793069/50‑1 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:7204496/50‑1 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:4379477/50‑1 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:6668638/1‑50 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:753202/50‑1 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:8147413/50‑1 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:2557185/50‑1 (MQ=255)              ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:2719337/1‑49 (MQ=255)              ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:7703265/1‑49 (MQ=255)              ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCg                                <  1:6067431/50‑1 (MQ=255)              ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCg                                <  1:7544431/50‑1 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              <  1:599894/50‑1 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGAAGGCCGa                               >  1:7524045/1‑49 (MQ=39)                  aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            >  1:5109627/1‑50 (MQ=255)                  aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            >  1:1361743/1‑50 (MQ=255)                   cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           >  1:6966412/1‑50 (MQ=255)                   cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           >  1:6750544/1‑50 (MQ=255)                   cGTTTGATCTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           <  1:5198247/50‑1 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          >  1:7551539/1‑50 (MQ=255)                     tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         <  1:7904805/50‑1 (MQ=255)                     tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         <  1:22742/50‑1 (MQ=255)                       tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGtt                        <  1:4073114/49‑1 (MQ=255)                       tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       >  1:3121584/1‑50 (MQ=255)                        gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       <  1:3615083/49‑1 (MQ=255)                        gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       <  1:2085795/49‑1 (MQ=255)                        gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       <  1:6999330/49‑1 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:2129005/1‑50 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:1185877/1‑50 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:4182733/1‑50 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:4080375/1‑50 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      <  1:3176292/49‑1 (MQ=255)                          ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    <  1:6027807/50‑1 (MQ=255)                          ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    <  1:5371622/50‑1 (MQ=255)                          ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    >  1:5551030/1‑50 (MQ=255)                          ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    <  1:1376434/50‑1 (MQ=255)                           tttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   >  1:8211468/1‑50 (MQ=255)                            ttcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   >  1:7596724/1‑49 (MQ=255)                            ttcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   <  1:5727102/49‑1 (MQ=255)                            ttcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   <  1:6453981/49‑1 (MQ=255)                               ttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:4043514/1‑49 (MQ=255)                               ttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:5930405/1‑49 (MQ=255)                               ttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:1031641/1‑49 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAtt               <  1:1618323/48‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:1165793/47‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:5836932/47‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:4891158/47‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              <  1:1655034/49‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              <  1:5964791/49‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              <  1:2364389/49‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             >  1:2131016/1‑50 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             <  1:5719469/49‑1 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             <  1:5351515/49‑1 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            <  1:4780318/50‑1 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            <  1:8130784/50‑1 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           <  1:7990533/50‑1 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           >  1:8188015/1‑50 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           >  1:1933507/1‑50 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           >  1:8123174/1‑50 (MQ=255)                                    tACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:7117857/1‑50 (MQ=255)                                    tACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:1227002/1‑50 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:3041209/1‑48 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          <  1:2245885/48‑1 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTc         <  1:4449562/49‑1 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCa        <  1:7292163/50‑1 (MQ=255)                                       cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAg       <  1:5975287/50‑1 (MQ=255)                                         aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  1:3777790/1‑50 (MQ=255)                                          tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    <  1:6183342/50‑1 (MQ=255)                                          tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    >  1:3718270/1‑50 (MQ=255)                                          tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    <  1:6417561/50‑1 (MQ=255)                                           ttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt   >  1:1591506/1‑50 (MQ=255)                                            tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt   <  1:1847017/49‑1 (MQ=255)                                            tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt   <  1:4668974/49‑1 (MQ=255)                                            tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt   <  1:904275/49‑1 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  1:7636689/49‑1 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:806662/1‑49 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:3153586/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 < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑