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	98.9%	357.3 / NA	87	P163S (CCA→TCA)	asnS	asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (45/41);  minor base T (0/1);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGC  >  NC_000913/988452‑988545                                                |                                               gcgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACg                                              >  1:1114428/1‑50 (MQ=255) gcgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACg                                              >  1:825700/1‑50 (MQ=255)  cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             >  1:5752684/1‑50 (MQ=255)  cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             >  1:119784/1‑50 (MQ=255)  cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             <  1:2628258/50‑1 (MQ=255)  cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                             >  1:4742725/1‑50 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  1:3413100/50‑1 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  1:843773/50‑1 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            >  1:5243264/1‑50 (MQ=255)   ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            >  1:2128758/1‑50 (MQ=255)    gAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            >  1:1101934/1‑49 (MQ=255)     aaCATTTCGCCTGCTCCTTCGGTATCAGATGCGGTAATCAGTGTCGTTGa                                          <  1:8070085/50‑1 (MQ=255)     aaCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTg                                           >  1:1218955/1‑49 (MQ=255)     aaCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTg                                           >  1:742205/1‑49 (MQ=255)     aaCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                          <  1:1858339/50‑1 (MQ=255)     aaCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                          <  1:4120860/50‑1 (MQ=255)     aaCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                          <  1:3007445/50‑1 (MQ=255)       cATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaaa                                        >  1:4685662/1‑50 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       <  1:8093241/50‑1 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       <  1:2727422/50‑1 (MQ=255)        aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       <  1:1211405/50‑1 (MQ=255)         tttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      <  1:3578792/50‑1 (MQ=255)          ttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  1:2623187/1‑49 (MQ=255)          ttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  1:4218994/1‑49 (MQ=255)           tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAccc                                     <  1:5724810/49‑1 (MQ=255)           tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAccc                                     <  1:2814825/49‑1 (MQ=255)            cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCa                                    <  1:1966006/49‑1 (MQ=255)            cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCa                                    <  1:2675208/49‑1 (MQ=255)            cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCa                                    <  1:5173408/49‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCaga                                  <  1:3823745/50‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCaga                                  <  1:6437060/50‑1 (MQ=255)             gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:6940143/1‑49 (MQ=255)               cTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:4965572/1‑47 (MQ=255)                tGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                                >  1:993209/1‑49 (MQ=255)                 gCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAgaagaa                              >  1:6014359/1‑50 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaaga                               <  1:3312652/48‑1 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaaga                               <  1:3036678/48‑1 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaaga                               <  1:2734266/48‑1 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAAt                             <  1:384460/50‑1 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAAt                             >  1:6375198/1‑50 (MQ=255)                  cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAAt                             <  1:5648166/50‑1 (MQ=255)                   aCCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATc                            >  1:4177669/1‑50 (MQ=255)                    ccTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATcc                           >  1:7570249/1‑50 (MQ=255)                     cTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATccc                          >  1:1714032/1‑50 (MQ=255)                     cTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATccc                          >  1:5528471/1‑50 (MQ=255)                     cTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATccc                          >  1:6128418/1‑50 (MQ=255)                      ttCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         >  1:3349086/1‑50 (MQ=255)                      ttCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         >  1:3382557/1‑50 (MQ=255)                       tCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        >  1:770346/1‑50 (MQ=255)                        cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         <  1:2877924/48‑1 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGc                       <  1:8319966/49‑1 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      <  1:6361998/50‑1 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      <  1:5922351/50‑1 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  1:1669787/1‑50 (MQ=255)                         ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  1:247902/1‑50 (MQ=255)                          gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:2511911/50‑1 (MQ=255)                          gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     >  1:85084/1‑50 (MQ=255)                          gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:4587118/50‑1 (MQ=255)                           tATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCg                    <  1:4686732/50‑1 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  1:1714358/1‑49 (MQ=255)                             tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   <  1:722405/49‑1 (MQ=255)                               aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:4562067/1‑49 (MQ=255)                               aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:8021986/1‑49 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:6133595/1‑47 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:4937249/1‑47 (MQ=255)                                 aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              <  1:4449195/50‑1 (MQ=255)                                   gCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              >  1:2641500/1‑48 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           <  1:3341981/50‑1 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaa            <  1:3708670/49‑1 (MQ=255)                                    cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaa            <  1:7252414/49‑1 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          <  1:4918280/50‑1 (MQ=255)                                     ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          >  1:4010390/1‑50 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        <  1:7154238/50‑1 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        >  1:7959098/1‑50 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        >  1:4725611/1‑50 (MQ=255)                                       tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        <  1:1507651/50‑1 (MQ=255)                                         aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGAt       <  1:7597195/49‑1 (MQ=255)                                         aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:5837349/1‑50 (MQ=255)                                         aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:1357687/1‑50 (MQ=255)                                           cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:5610352/1‑48 (MQ=255)                                           cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:3926340/1‑48 (MQ=255)                                           cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:471041/1‑48 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  1:5386972/49‑1 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  1:1195990/49‑1 (MQ=255)                                            aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   >  1:7101307/1‑50 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  <  1:1113674/50‑1 (MQ=255)                                             gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  >  1:4314673/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 < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑