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

GCGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGC  >  NC_000913/988452‑988545                                                |                                               gcgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACg                                              >  2:2536910/1‑50 (MQ=255)  cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             <  1:4371729/50‑1 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  2:4441914/50‑1 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  2:4174039/50‑1 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            >  1:1483467/1‑50 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            >  2:1485096/1‑50 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  2:3977317/50‑1 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  2:2849691/50‑1 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  2:3086024/50‑1 (MQ=255)    gAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTg                                           <  2:3678980/50‑1 (MQ=255)       cATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaaa                                        >  1:2709709/1‑50 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       <  2:4043857/50‑1 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       >  1:1856071/1‑50 (MQ=255)          ttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  2:771482/1‑49 (MQ=255)           tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAccc                                     <  1:1801552/49‑1 (MQ=255)           tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAccc                                     <  1:2160869/49‑1 (MQ=255)           tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAccc                                     <  1:3956697/49‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   <  2:915620/49‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   <  2:1930671/49‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  2:2411913/1‑49 (MQ=255)              ccTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaa                                 <  1:3727837/50‑1 (MQ=255)              ccTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  2:4039739/1‑48 (MQ=255)               cTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                                <  2:1800265/50‑1 (MQ=255)               cTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                                <  1:4566556/50‑1 (MQ=255)               cTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                                <  2:4520612/50‑1 (MQ=255)               cTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:3638566/1‑47 (MQ=255)                 gCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAgaagaa                              >  1:1034130/1‑50 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaaga                               <  1:426087/48‑1 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAgaagaa                              <  2:2265688/49‑1 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAgaaaaa                              <  2:239878/49‑1 (MQ=39)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAAt                             >  2:3365049/1‑50 (MQ=255)                   aCCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATc                            >  1:2278342/1‑50 (MQ=255)                   aCCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATc                            <  1:984147/50‑1 (MQ=255)                   aCCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATc                            >  1:1399801/1‑50 (MQ=255)                   aCCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATc                            <  1:68047/50‑1 (MQ=255)                     cTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATccc                          <  2:405700/50‑1 (MQ=255)                     cTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATccc                          >  2:3236801/1‑50 (MQ=255)                     cTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATccc                          <  2:3454450/50‑1 (MQ=255)                     cTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATccc                          >  1:4445077/1‑50 (MQ=255)                     cTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATccc                          >  1:4061456/1‑50 (MQ=255)                      ttCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         >  1:1213176/1‑50 (MQ=255)                       tCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        >  1:643910/1‑50 (MQ=255)                       tCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        >  1:1347282/1‑50 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGc                       >  1:2684259/1‑50 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGc                       >  2:453132/1‑50 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGc                       <  1:1453376/49‑1 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGc                       <  1:39637/49‑1 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      <  1:4003385/50‑1 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  2:2035244/1‑50 (MQ=255)                          gTATCAGATGCGGTAATCAGTTACGTTGAAACCCAGAAGAATCCCTGCTc                     <  2:1877948/50‑1 (MQ=255)                          gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:411325/50‑1 (MQ=255)                          gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:2827389/50‑1 (MQ=255)                          gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  2:2108258/50‑1 (MQ=255)                           tATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCg                    >  2:3775748/1‑50 (MQ=255)                           tATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCg                    <  1:2533951/50‑1 (MQ=255)                           tATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCg                    >  2:4561736/1‑50 (MQ=255)                           tATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCg                    >  1:368376/1‑50 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  <  1:1554037/50‑1 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  <  1:3284592/50‑1 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  <  2:2497871/50‑1 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  2:3571720/1‑49 (MQ=255)                               aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTaa                <  1:2964403/50‑1 (MQ=255)                               aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 <  1:1017471/49‑1 (MQ=255)                               aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  2:3248066/1‑49 (MQ=255)                               aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 <  1:2320093/49‑1 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              <  1:1989312/50‑1 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              <  2:1161190/50‑1 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           <  1:443883/50‑1 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           <  1:3282045/50‑1 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           >  1:1648288/1‑50 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaa            <  2:2477480/49‑1 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          <  2:3535811/50‑1 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          <  2:377801/50‑1 (MQ=255)                                      gTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCg         <  1:2454339/50‑1 (MQ=255)                                      gTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCg         <  2:1916860/50‑1 (MQ=255)                                      gTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCg         <  1:770482/50‑1 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        <  2:689683/50‑1 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        <  2:1837793/50‑1 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        >  2:3308289/1‑50 (MQ=255)                                        aaTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGAt       >  2:259662/1‑50 (MQ=255)                                         aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGAt       <  2:1701653/49‑1 (MQ=255)                                          tCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:978815/1‑49 (MQ=255)                                          tCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGc     <  1:3675688/50‑1 (MQ=255)                                           cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:3933891/1‑48 (MQ=255)                                            aGTGACGTTGAAGCCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  1:924632/49‑1 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  2:4529049/49‑1 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   <  1:2287809/50‑1 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  >  2:1835383/1‑50 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  >  2:3401454/1‑50 (MQ=255)                                                |                                               GCGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGC  >  NC_000913/988452‑988545

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