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

GCGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGC  >  NC_000913/988452‑988545                                                |                                               gcgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACg                                              <  1:2890223/50‑1 (MQ=255)  cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             >  1:343148/1‑50 (MQ=255)  cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             >  1:2331423/1‑50 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            >  1:3498921/1‑50 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  1:4400465/50‑1 (MQ=255)    gAACATTTCGCCTGCACCTTCGGTATCAGCTGCGGAAATCAGCGACGtt                                            >  1:137470/1‑49 (MQ=37)     aaCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTg                                           >  1:4690503/1‑49 (MQ=255)      aCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaa                                         <  1:2702891/50‑1 (MQ=255)      aCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                          >  1:4134327/1‑49 (MQ=255)      aCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                          >  1:598245/1‑49 (MQ=255)      aCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                          >  1:4182521/1‑49 (MQ=255)      aCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                          >  1:4353830/1‑49 (MQ=255)       cATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaaa                                        <  1:6367172/50‑1 (MQ=255)       cATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaaa                                        <  1:5520664/50‑1 (MQ=255)       cATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaaa                                        >  1:3469376/1‑50 (MQ=255)       cATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaaa                                        <  1:493229/50‑1 (MQ=255)       cATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaaa                                        <  1:437278/50‑1 (MQ=255)       cATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaaa                                        >  1:1468481/1‑50 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       >  1:4743483/1‑50 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       >  1:2939380/1‑50 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       <  1:1542056/50‑1 (MQ=255)         tttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      <  1:394795/50‑1 (MQ=255)          ttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  1:2760696/1‑49 (MQ=255)           tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAccc                                     <  1:6685346/49‑1 (MQ=255)           tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAccc                                     <  1:2603513/49‑1 (MQ=255)           tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAccc                                     <  1:5070885/49‑1 (MQ=255)           tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  1:2845518/1‑48 (MQ=255)           tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  1:5073618/1‑48 (MQ=255)            cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCa                                    <  1:6529688/49‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   <  1:888376/49‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:5001822/1‑49 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   <  1:5889324/49‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:2544992/1‑49 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:3077882/1‑49 (MQ=255)              ccTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaa                                 <  1:6186708/50‑1 (MQ=255)              ccTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:121733/1‑48 (MQ=255)               cTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:280076/1‑47 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaaga                               <  1:3564009/48‑1 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAgaagaa                              <  1:2538322/49‑1 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAAt                             >  1:4114943/1‑50 (MQ=255)                   aCCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATc                            <  1:3963939/50‑1 (MQ=255)                       tCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        >  1:3371527/1‑50 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         <  1:6320932/48‑1 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         <  1:4494527/48‑1 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         <  1:2712803/48‑1 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  1:2118883/1‑50 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      <  1:5250245/50‑1 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      <  1:442399/50‑1 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      <  1:3818422/50‑1 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  1:3699842/1‑50 (MQ=255)                          gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:4396535/50‑1 (MQ=255)                          gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     >  1:475483/1‑50 (MQ=255)                            aTCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  1:5850328/1‑50 (MQ=255)                            aTCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  1:4234714/1‑50 (MQ=255)                             tCAGATGTGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  <  1:2064195/50‑1 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   <  1:32950/49‑1 (MQ=255)                              cAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  >  1:4385215/1‑49 (MQ=255)                               aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 <  1:5432176/49‑1 (MQ=255)                                gATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTaaa               <  1:380032/50‑1 (MQ=255)                                gATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:1975044/1‑48 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:5965206/1‑47 (MQ=255)                                  tGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaaga             <  1:6120497/50‑1 (MQ=255)                                  tGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              >  1:1170359/1‑49 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           <  1:3478549/50‑1 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           >  1:6656405/1‑50 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaa            <  1:3245155/49‑1 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaa            <  1:6087547/49‑1 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          <  1:2959370/50‑1 (MQ=255)                                      gTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCg         <  1:5590111/50‑1 (MQ=255)                                      gTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCg         <  1:6095615/50‑1 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        <  1:2237269/50‑1 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        <  1:2234096/50‑1 (MQ=255)                                        aaTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGAt       >  1:1339066/1‑50 (MQ=255)                                         aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGAt       <  1:3953481/49‑1 (MQ=255)                                         aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGAt       <  1:4628735/49‑1 (MQ=255)                                         aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGAt       <  1:2341736/49‑1 (MQ=255)                                         aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:1423706/1‑50 (MQ=255)                                          tCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGc     <  1:3552518/50‑1 (MQ=255)                                           cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:4196811/1‑48 (MQ=255)                                           cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:824750/1‑48 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  1:6426988/49‑1 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   >  1:3048025/1‑50 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   <  1:3159805/50‑1 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  >  1:1508040/1‑50 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  >  1:3298362/1‑50 (MQ=255)                                                |                                               GCGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGC  >  NC_000913/988452‑988545

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