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

GCGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGC  >  NC_000913/988452‑988545                                                |                                               gcgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACg                                              <  1:315008/50‑1 (MQ=255)  cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             >  1:1071213/1‑50 (MQ=255)  cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             <  1:8301460/50‑1 (MQ=255)  cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             >  1:5894633/1‑50 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  1:3936903/50‑1 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  1:7881153/50‑1 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  1:7131758/50‑1 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  1:3275865/50‑1 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            >  1:5986514/1‑50 (MQ=255)    gAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTg                                           <  1:7198222/50‑1 (MQ=255)      aCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaa                                         <  1:2537694/50‑1 (MQ=255)      aCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                          >  1:5446731/1‑49 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       >  1:5536782/1‑50 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       >  1:3840402/1‑50 (MQ=255)          ttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  1:7266552/1‑49 (MQ=255)           tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAccc                                     <  1:2355325/49‑1 (MQ=255)            cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  1:4649378/1‑47 (MQ=255)            cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCa                                    <  1:330313/49‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:5301324/1‑49 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   <  1:3554423/49‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:6265641/1‑49 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   <  1:4019637/49‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   <  1:3742422/49‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   <  1:2736371/49‑1 (MQ=255)              ccTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaa                                 <  1:1525730/50‑1 (MQ=255)               cTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                                <  1:4655046/50‑1 (MQ=255)               cTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:3315350/1‑47 (MQ=255)                tGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                                >  1:411039/1‑49 (MQ=255)                tGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                                >  1:6716661/1‑49 (MQ=255)                tGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                                >  1:3350729/1‑49 (MQ=255)                 gCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAgaagaa                              >  1:7144410/1‑50 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaaga                               <  1:4286097/48‑1 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAAt                             >  1:4370414/1‑50 (MQ=255)                   aCCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATc                            <  1:5969561/50‑1 (MQ=255)                    ccTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATcc                           <  1:3678397/50‑1 (MQ=255)                      ttCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         >  1:1050953/1‑50 (MQ=255)                      ttCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         >  1:6082359/1‑50 (MQ=255)                      ttCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         >  1:6351177/1‑50 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         <  1:6031955/48‑1 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         <  1:6276937/48‑1 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        <  1:7401571/49‑1 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        <  1:8298922/49‑1 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        <  1:1327411/49‑1 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGc                       >  1:6450138/1‑50 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGc                       <  1:6133151/49‑1 (MQ=255)                          gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     >  1:523653/1‑50 (MQ=255)                           tATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCg                    >  1:5265928/1‑50 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  <  1:804368/50‑1 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  <  1:3680909/50‑1 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   <  1:5465755/49‑1 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  1:1181402/1‑49 (MQ=255)                              cAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  >  1:2290838/1‑49 (MQ=255)                               aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTaa                <  1:5861712/50‑1 (MQ=255)                               aGACGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:1943432/1‑49 (MQ=255)                                gATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTaaa               <  1:8155786/50‑1 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              <  1:811529/50‑1 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           <  1:7736365/50‑1 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           >  1:3555027/1‑50 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaa            <  1:81701/49‑1 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          >  1:6413270/1‑50 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          <  1:5207699/50‑1 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          <  1:3640938/50‑1 (MQ=255)                                      gTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCg         <  1:2477429/50‑1 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        >  1:1494345/1‑50 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        >  1:8100571/1‑50 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        >  1:7682007/1‑50 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        >  1:4401243/1‑50 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        >  1:5275536/1‑50 (MQ=255)                                         aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGACGCGATg      <  1:8003670/50‑1 (MQ=255)                                         aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      <  1:7972345/50‑1 (MQ=255)                                         aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      <  1:192197/50‑1 (MQ=255)                                          tCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:2479920/1‑49 (MQ=255)                                          tCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGc     <  1:4335357/50‑1 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  1:4528540/49‑1 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  1:1343444/49‑1 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   >  1:5062407/1‑50 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   <  1:3464849/50‑1 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   >  1:3326164/1‑50 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   >  1:5715891/1‑50 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  >  1:4828863/1‑50 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  >  1:4224140/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 < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑