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

TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048                                                |                                              tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             <  1:1837378/50‑1 (MQ=255)  cgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            >  1:4649255/1‑50 (MQ=255)   ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCtt                                           <  1:2930298/50‑1 (MQ=255)    gAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:1737315/1‑50 (MQ=255)    gAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          <  1:1549961/50‑1 (MQ=255)     aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:5126919/1‑49 (MQ=255)     aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:4920554/1‑49 (MQ=255)      gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCca                                        <  1:100414/50‑1 (MQ=255)       tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:983322/1‑50 (MQ=255)       tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       <  1:3283124/50‑1 (MQ=255)       tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       <  1:2489682/50‑1 (MQ=255)        cAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:4880244/1‑49 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagca                                     <  1:1232057/50‑1 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagca                                     <  1:2299770/50‑1 (MQ=255)         aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTCCcag                                       >  1:3580738/1‑48 (MQ=39)           tGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAgg                                   <  1:3861643/50‑1 (MQ=255)           tGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAgg                                   >  1:3938805/1‑50 (MQ=255)           tGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAgg                                   <  1:475340/50‑1 (MQ=255)            gACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  >  1:3708182/1‑50 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:1691638/50‑1 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:3776031/1‑50 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:425870/1‑50 (MQ=255)             aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  <  1:2866152/49‑1 (MQ=255)              ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:110333/1‑49 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:2557739/1‑49 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:4945564/1‑49 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:3737999/1‑49 (MQ=255)                gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              <  1:313643/50‑1 (MQ=255)                 cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATa                             <  1:3925241/50‑1 (MQ=255)                  aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            >  1:5082167/1‑50 (MQ=255)                   cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           >  1:1700211/1‑50 (MQ=255)                     tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  1:195342/1‑50 (MQ=255)                     tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  1:3805050/1‑50 (MQ=255)                      ttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  1:3583784/1‑49 (MQ=255)                      ttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  1:1148660/1‑49 (MQ=255)                       tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGtt                        <  1:1931032/49‑1 (MQ=255)                         aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     >  1:1396969/1‑50 (MQ=255)                          ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    <  1:262603/50‑1 (MQ=255)                          ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    >  1:1266352/1‑50 (MQ=255)                           tttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   >  1:2106155/1‑50 (MQ=255)                            ttcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   <  1:1459003/49‑1 (MQ=255)                               ttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:1619502/1‑49 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAtt               <  1:1027389/48‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:2270404/47‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:1217202/47‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:303851/47‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                <  1:3485836/47‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              <  1:3077221/49‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              <  1:12173/49‑1 (MQ=255)                                 cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             >  1:3790598/1‑50 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt             <  1:4819865/49‑1 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            <  1:1779921/50‑1 (MQ=255)                                  aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            <  1:158013/50‑1 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           <  1:575762/50‑1 (MQ=255)                                   cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           <  1:1467953/50‑1 (MQ=255)                                     aCCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:1502050/1‑49 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:67514/1‑48 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          <  1:904266/48‑1 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTc         <  1:313697/49‑1 (MQ=255)                                      ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCa        <  1:546899/50‑1 (MQ=255)                                       cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:2909255/1‑47 (MQ=255)                                       cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:605563/1‑47 (MQ=255)                                       cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAg       <  1:2130943/50‑1 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      <  1:739865/50‑1 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:1058807/1‑50 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:4569967/1‑50 (MQ=255)                                        gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:2675789/1‑50 (MQ=255)                                         aTTTCGGGTCTTTCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     <  1:4548711/50‑1 (MQ=39)                                         aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     <  1:361399/50‑1 (MQ=255)                                            tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt   <  1:2644207/49‑1 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:5136257/1‑49 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  1:289706/49‑1 (MQ=255)                                             cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  1:3169661/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 < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑