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

TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048                                                |                                              tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             <  1:3970167/50‑1 (MQ=255) tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             >  1:3536702/1‑50 (MQ=255)  cgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            >  1:2810317/1‑50 (MQ=255)   ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCtt                                           >  1:1205992/1‑50 (MQ=255)   ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCtt                                           <  1:3427169/50‑1 (MQ=255)     aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:5952009/1‑49 (MQ=255)      gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:4313286/1‑48 (MQ=255)      gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:366444/1‑48 (MQ=255)      gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCca                                        <  1:3279449/50‑1 (MQ=255)       tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       <  1:956014/50‑1 (MQ=255)       tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       <  1:5797030/50‑1 (MQ=255)       tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       <  1:4451138/50‑1 (MQ=255)       tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:4023712/1‑50 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagca                                     <  1:7171022/50‑1 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagc                                      <  1:6673153/49‑1 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagc                                      <  1:5254590/49‑1 (MQ=255)          gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    <  1:1684407/50‑1 (MQ=255)            gACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  >  1:2581912/1‑50 (MQ=255)            gACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  >  1:2779595/1‑50 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:4868280/50‑1 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:907451/1‑50 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:4101098/1‑50 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:5120312/1‑50 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  <  1:4986214/49‑1 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  <  1:637192/49‑1 (MQ=255)              ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:6648227/1‑49 (MQ=255)               cGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:4667535/1‑50 (MQ=255)               cGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:2249748/1‑50 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:3755387/1‑49 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:4997273/1‑49 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:4865545/1‑49 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              <  1:5821728/50‑1 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              <  1:6671959/50‑1 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              <  1:1040678/50‑1 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              <  1:4792079/50‑1 (MQ=255)                 cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              >  1:4175536/1‑49 (MQ=255)                 cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATa                             <  1:5859680/50‑1 (MQ=255)                 cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATa                             <  1:7453840/50‑1 (MQ=255)                  aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            <  1:6174997/50‑1 (MQ=255)                  aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            <  1:2608032/50‑1 (MQ=255)                  aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            >  1:3100780/1‑50 (MQ=255)                   cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           <  1:2567214/50‑1 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          <  1:2311927/50‑1 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          <  1:6163228/50‑1 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          <  1:5712614/50‑1 (MQ=255)                    gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          <  1:4355633/50‑1 (MQ=255)                       tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       >  1:71919/1‑50 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:51912/1‑50 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     <  1:5864434/50‑1 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     <  1:2012390/50‑1 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:3335979/1‑50 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     <  1:1445502/50‑1 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     <  1:5585054/50‑1 (MQ=255)                         aTTTACTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     <  1:6871268/50‑1 (MQ=255)                            ttcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   <  1:1797000/49‑1 (MQ=255)                             tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCg                  <  1:5646417/49‑1 (MQ=255)                             tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGa                 >  1:5995538/1‑50 (MQ=255)                              cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:4442336/1‑50 (MQ=255)                              cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:5838941/1‑50 (MQ=255)                               ttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:5370144/1‑49 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAtt               <  1:3733665/48‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:5393883/47‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:3827391/47‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:87384/47‑1 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             <  1:4557365/49‑1 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            <  1:637750/50‑1 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            <  1:5926057/50‑1 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            <  1:7331541/50‑1 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            >  1:3836153/1‑50 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            >  1:1541610/1‑50 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           <  1:238297/50‑1 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           >  1:7170100/1‑50 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           <  1:7045094/50‑1 (MQ=255)                                    tACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:6360902/1‑50 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:1210774/1‑48 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:7176934/1‑48 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTc         <  1:3313785/49‑1 (MQ=255)                                       cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAg       <  1:5907655/50‑1 (MQ=255)                                       cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAg       <  1:4768009/50‑1 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      <  1:7098324/50‑1 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:1913978/1‑50 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:4317383/1‑50 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:2876957/1‑50 (MQ=255)                                         aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  1:7036153/1‑50 (MQ=255)                                         aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  1:5568781/1‑50 (MQ=255)                                         aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  1:5718830/1‑50 (MQ=255)                                            tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt   <  1:4036129/49‑1 (MQ=255)                                            tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:4051641/1‑50 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:7160897/1‑49 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  1:5565010/49‑1 (MQ=255)                                                |                                              TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048

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