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

GCGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGC  >  NC_000913/988452‑988545                                                |                                               gcgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACg                                              >  1:6353149/1‑50 (MQ=255)  cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             >  1:6101531/1‑50 (MQ=255)    gAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGAGGTTg                                           <  1:3658179/50‑1 (MQ=255)    gAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            >  1:735890/1‑49 (MQ=255)    gAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            >  1:5500467/1‑49 (MQ=255)    gAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTg                                           <  1:622010/50‑1 (MQ=255)     aaCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTg                                           >  1:2683101/1‑49 (MQ=255)     aaCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                          <  1:68942/50‑1 (MQ=255)       cATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaaa                                        <  1:2999714/50‑1 (MQ=255)       cATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaaa                                        <  1:2976755/50‑1 (MQ=255)        aTTTCGCCTGCGCCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       <  1:2616129/50‑1 (MQ=39)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       >  1:3282447/1‑50 (MQ=255)          ttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  1:4446662/1‑49 (MQ=255)           tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAccc                                     <  1:3946911/49‑1 (MQ=255)           tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  1:3327477/1‑48 (MQ=255)            cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  1:2532442/1‑47 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCaga                                  <  1:615502/50‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCaga                                  <  1:4082681/50‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:2523759/1‑49 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:3326212/1‑49 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:4386602/1‑49 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   <  1:601605/49‑1 (MQ=255)              ccTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaa                                 <  1:1582966/50‑1 (MQ=255)               cTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                                <  1:5293586/50‑1 (MQ=255)                tGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                                >  1:2513489/1‑49 (MQ=255)                   aCCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATc                            >  1:5812765/1‑50 (MQ=255)                    ccTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATcc                           <  1:4249163/50‑1 (MQ=255)                      ttCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         >  1:5371642/1‑50 (MQ=255)                      ttCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         >  1:5717380/1‑50 (MQ=255)                       tCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        >  1:478999/1‑50 (MQ=255)                       tCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        >  1:1766818/1‑50 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         <  1:6426162/48‑1 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGc                       >  1:413316/1‑50 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGc                       <  1:5769031/49‑1 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  1:1110038/1‑50 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  1:2245016/1‑50 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  1:4658000/1‑50 (MQ=255)                          gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:6281967/50‑1 (MQ=255)                          gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:5313118/50‑1 (MQ=255)                           tATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCg                    <  1:5744346/50‑1 (MQ=255)                            aTCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  1:2759475/1‑50 (MQ=255)                            aTCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  1:5124162/1‑50 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  <  1:4734683/50‑1 (MQ=255)                                gATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTaaa               <  1:2029888/50‑1 (MQ=255)                                gATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTaaa               <  1:1538874/50‑1 (MQ=255)                                gATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTaaa               <  1:4076822/50‑1 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:5750789/1‑47 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:617849/1‑47 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:5460876/1‑47 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              <  1:1051327/50‑1 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              <  1:4357685/50‑1 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              <  1:4176724/50‑1 (MQ=255)                                  tGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaaga             <  1:2344492/50‑1 (MQ=255)                                  tGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaaga             <  1:4895862/50‑1 (MQ=255)                                  tGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaaga             <  1:4592900/50‑1 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           >  1:1857669/1‑50 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           <  1:642164/50‑1 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           <  1:3589217/50‑1 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaa            <  1:1991123/49‑1 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          >  1:5972659/1‑50 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          >  1:5332713/1‑50 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          <  1:1476385/50‑1 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          <  1:2515414/50‑1 (MQ=255)                                      gTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          >  1:1850561/1‑49 (MQ=255)                                      gTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCg         <  1:3673206/50‑1 (MQ=255)                                         aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGAt       <  1:1790485/49‑1 (MQ=255)                                         aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:6136637/1‑50 (MQ=255)                                          tCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGc     <  1:4739757/50‑1 (MQ=255)                                          tCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGc     <  1:1285091/50‑1 (MQ=255)                                           cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:4409971/1‑48 (MQ=255)                                           cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:1375419/1‑48 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  1:3255783/49‑1 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   >  1:2086077/1‑50 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  <  1:4601272/50‑1 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  <  1:3467086/50‑1 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  >  1:3535203/1‑50 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  >  1:1351547/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: ‑