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

TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048                                                |                                              tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             <  1:5854333/50‑1 (MQ=255) tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             >  1:5927393/1‑50 (MQ=255) tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             <  1:1177661/50‑1 (MQ=255)  cgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            >  1:2321024/1‑50 (MQ=255)  cgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            >  1:3540003/1‑50 (MQ=255)   ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCtt                                           >  1:114283/1‑50 (MQ=255)     aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:2729031/1‑49 (MQ=255)      gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCca                                        <  1:4167705/50‑1 (MQ=255)       tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       <  1:3143165/50‑1 (MQ=255)        cAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:6152275/1‑49 (MQ=255)        cAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGCCTTCcag                                       >  1:1781737/1‑49 (MQ=39)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagca                                     <  1:3819530/50‑1 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagc                                      <  1:3533611/49‑1 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:4030082/1‑48 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:6407260/1‑48 (MQ=255)          gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    <  1:1202086/50‑1 (MQ=255)            gACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  >  1:1686459/1‑50 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:3611791/50‑1 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:1207364/50‑1 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  <  1:4020216/49‑1 (MQ=255)              ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCg                                <  1:2979852/50‑1 (MQ=255)              ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCg                                <  1:5487001/50‑1 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:6086499/1‑49 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              <  1:5294533/50‑1 (MQ=255)                 cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              >  1:3526854/1‑49 (MQ=255)                 cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATa                             <  1:3958302/50‑1 (MQ=255)                 cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATa                             <  1:710340/50‑1 (MQ=255)                 cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATa                             <  1:1450389/50‑1 (MQ=255)                  aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            >  1:2030300/1‑50 (MQ=255)                   cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           >  1:4173149/1‑50 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          >  1:1267469/1‑50 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          <  1:470376/50‑1 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          <  1:3001809/50‑1 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          <  1:795755/50‑1 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          >  1:6141942/1‑50 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          >  1:1643025/1‑50 (MQ=255)                     tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         <  1:3269396/50‑1 (MQ=255)                       tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGtt                        <  1:6436153/49‑1 (MQ=255)                       tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       >  1:1129737/1‑50 (MQ=255)                        gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      >  1:4463806/1‑50 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     <  1:1801389/50‑1 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:6632881/1‑50 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:529862/1‑50 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      <  1:4671567/49‑1 (MQ=255)                          ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    <  1:1412973/50‑1 (MQ=255)                          ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    >  1:2036732/1‑50 (MQ=255)                            ttcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   <  1:6697741/49‑1 (MQ=255)                             tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGa                 >  1:1241358/1‑50 (MQ=255)                              cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGa                 <  1:3960591/49‑1 (MQ=255)                              cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:303576/1‑50 (MQ=255)                              cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:950869/1‑50 (MQ=255)                               ttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:5030997/1‑49 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAtt               <  1:3748416/48‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAtt               <  1:1789072/48‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAtt               <  1:5035059/48‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:2365055/47‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              <  1:6712239/49‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              <  1:5821155/49‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              <  1:5016685/49‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             >  1:6327079/1‑50 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             >  1:1114726/1‑50 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             >  1:4299947/1‑50 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           >  1:4792235/1‑50 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           >  1:3109281/1‑50 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          <  1:1246252/48‑1 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          <  1:2588113/48‑1 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:2809871/1‑48 (MQ=255)                                       cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:1025216/1‑47 (MQ=255)                                       cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAg       <  1:953637/50‑1 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      <  1:2852890/50‑1 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      <  1:3387073/50‑1 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      <  1:4683493/50‑1 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:5050530/1‑50 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      <  1:5370670/50‑1 (MQ=255)                                         aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     <  1:1121073/50‑1 (MQ=255)                                         aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  1:5249632/1‑50 (MQ=255)                                          tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    <  1:203259/50‑1 (MQ=255)                                          tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    <  1:6691335/50‑1 (MQ=255)                                            tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt   <  1:1784734/49‑1 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  1:6249820/49‑1 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:2454142/1‑49 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  1:1641044/49‑1 (MQ=255)                                                |                                              TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048

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