Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	988,499	G→A	P163S (CCA→TCA)	asnS ←	asparagine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	988,499	0	G	A	100.0%	281.6 / NA	73	P163S (CCA→TCA)	asnS	asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/39);  total (34/39)

GCGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGC  >  NC_000913/988452‑988545                                                |                                               gcgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACg                                              <  1:2396417/50‑1 (MQ=255) gcgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACg                                              <  1:2367163/50‑1 (MQ=255)  cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             <  1:6072312/50‑1 (MQ=255)  cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             <  1:6353368/50‑1 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  1:6448026/50‑1 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            >  1:4311897/1‑50 (MQ=255)    gAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTg                                           <  1:5397310/50‑1 (MQ=255)     aaCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                          <  1:2564568/50‑1 (MQ=255)     aaCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                          <  1:2409015/50‑1 (MQ=255)      aCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                          >  1:7471915/1‑49 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       <  1:6355230/50‑1 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       <  1:6879358/50‑1 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       <  1:492263/50‑1 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       >  1:6149692/1‑50 (MQ=255)         tttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  1:6702944/1‑50 (MQ=255)            cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCa                                    <  1:905411/49‑1 (MQ=255)            cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCa                                    <  1:2673993/49‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   <  1:6543648/49‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:5861379/1‑49 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:5679666/1‑49 (MQ=255)              ccTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:3087357/1‑48 (MQ=255)               cTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:4412298/1‑47 (MQ=255)                 gCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAgaagaa                              >  1:4382960/1‑50 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaaga                               <  1:4508386/48‑1 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAgaagaa                              <  1:7015627/49‑1 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAgaagaa                              <  1:4790278/49‑1 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAAt                             <  1:2732791/50‑1 (MQ=255)                    ccTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATcc                           >  1:4004112/1‑50 (MQ=255)                    ccTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATcc                           >  1:6812655/1‑50 (MQ=255)                     cTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATccc                          >  1:7257437/1‑50 (MQ=255)                      ttCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         >  1:6214318/1‑50 (MQ=255)                       tCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        >  1:295220/1‑50 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         <  1:1702489/48‑1 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         <  1:2608318/48‑1 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        <  1:3241123/49‑1 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  1:3822472/1‑50 (MQ=255)                          gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     >  1:4306646/1‑50 (MQ=255)                          gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:5163687/50‑1 (MQ=255)                           tATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCg                    <  1:1621954/50‑1 (MQ=255)                            aTCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  1:2608373/1‑50 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  <  1:7091090/50‑1 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  <  1:1707888/50‑1 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  1:6879212/1‑49 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   <  1:740238/49‑1 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  1:5642702/1‑49 (MQ=255)                              cAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  >  1:1830788/1‑49 (MQ=255)                               aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 <  1:6587864/49‑1 (MQ=255)                               aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 <  1:5358659/49‑1 (MQ=255)                               aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:1157231/1‑49 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              <  1:1216069/50‑1 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              <  1:3797646/50‑1 (MQ=255)                                  tGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaaga             <  1:5638671/50‑1 (MQ=255)                                  tGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              >  1:3346050/1‑49 (MQ=255)                                   gCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              >  1:64266/1‑48 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaa            <  1:5328856/49‑1 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          >  1:959553/1‑50 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          >  1:854144/1‑50 (MQ=255)                                      gTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCg         <  1:2957949/50‑1 (MQ=255)                                      gTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCg         <  1:1401022/50‑1 (MQ=255)                                        aaTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGAt       >  1:20835/1‑50 (MQ=255)                                          tCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:3941762/1‑49 (MQ=255)                                          tCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:1441627/1‑49 (MQ=255)                                           cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:378320/1‑48 (MQ=255)                                           cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:2439833/1‑48 (MQ=255)                                            aGTGACGTTGAAGCCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   <  1:5003367/50‑1 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  1:4010256/49‑1 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  1:3282366/49‑1 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   >  1:6721197/1‑50 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   >  1:2159498/1‑50 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   >  1:1584534/1‑50 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   <  1:4437/50‑1 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  <  1:3252373/50‑1 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  >  1:1839979/1‑50 (MQ=255)                                                |                                               GCGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGC  >  NC_000913/988452‑988545

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