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%	217.3 / NA	56	P163S (CCA→TCA)	asnS	asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/32);  total (24/32)

CGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGC  >  NC_000913/988453‑988545                                               |                                               cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             >  1:2576944/1‑50 (MQ=255) cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             >  1:2469951/1‑50 (MQ=255)  ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            >  1:136559/1‑50 (MQ=255)   gAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            >  1:3746383/1‑49 (MQ=255)   gAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTg                                           <  1:904409/50‑1 (MQ=255)     aCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaa                                         <  1:1126902/50‑1 (MQ=255)     aCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaa                                         <  1:1727736/50‑1 (MQ=255)      cATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaaa                                        >  1:3983417/1‑50 (MQ=255)         ttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  1:3486070/1‑49 (MQ=255)           cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCa                                    <  1:4287803/49‑1 (MQ=255)           cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCa                                    <  1:2082764/49‑1 (MQ=255)           cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCa                                    <  1:2175366/49‑1 (MQ=255)           cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCa                                    <  1:3272282/49‑1 (MQ=255)            gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCaga                                  <  1:3883889/50‑1 (MQ=255)            gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCaga                                  <  1:276781/50‑1 (MQ=255)            gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:3258922/1‑49 (MQ=255)               tGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                                >  1:122351/1‑49 (MQ=255)                 cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaaga                               <  1:4982317/48‑1 (MQ=255)                 cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAgaagaa                              <  1:3830102/49‑1 (MQ=255)                 cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAAt                             >  1:1701458/1‑50 (MQ=255)                 cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAAt                             <  1:4130928/50‑1 (MQ=255)                  aCCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATc                            >  1:1060393/1‑50 (MQ=255)                  aCCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATc                            <  1:1818015/50‑1 (MQ=255)                     ttCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         >  1:2662014/1‑50 (MQ=255)                     ttCGGTATCAGAAGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         >  1:4067472/1‑50 (MQ=255)                       cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        <  1:4088318/49‑1 (MQ=255)                       cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        <  1:1513069/49‑1 (MQ=255)                        ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      <  1:2254072/50‑1 (MQ=255)                         gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:5098859/50‑1 (MQ=255)                         gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     >  1:4019398/1‑50 (MQ=255)                         gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:3808424/50‑1 (MQ=255)                          tATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCg                    <  1:1679760/50‑1 (MQ=255)                          tATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCg                    <  1:1317933/50‑1 (MQ=255)                           aTCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  1:360013/1‑50 (MQ=255)                           aTCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  1:483234/1‑50 (MQ=255)                           aTCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  1:469264/1‑50 (MQ=255)                             cAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  >  1:3060500/1‑49 (MQ=255)                             cAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  >  1:2771629/1‑49 (MQ=255)                              aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:360018/1‑49 (MQ=255)                              aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:2222580/1‑49 (MQ=255)                              aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCACTGCTCGTTa                 <  1:4613360/49‑1 (MQ=255)                               gATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:756259/1‑48 (MQ=255)                                aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              <  1:1479166/50‑1 (MQ=255)                                   cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           <  1:1206004/50‑1 (MQ=255)                                   cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           >  1:3001933/1‑50 (MQ=255)                                   cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaa            <  1:5035223/49‑1 (MQ=255)                                    ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          <  1:3143515/50‑1 (MQ=255)                                      tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        >  1:2924738/1‑50 (MQ=255)                                      tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        >  1:784793/1‑50 (MQ=255)                                        aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGAt       <  1:1428599/49‑1 (MQ=255)                                        aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGAt       <  1:2480068/49‑1 (MQ=255)                                           aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  1:4343455/49‑1 (MQ=255)                                           aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   <  1:4512241/50‑1 (MQ=255)                                           aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   <  1:2208385/50‑1 (MQ=255)                                            gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  <  1:1714370/50‑1 (MQ=255)                                            gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  <  1:4322382/50‑1 (MQ=255)                                               |                                               CGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGC  >  NC_000913/988453‑988545

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