Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	707,102	A→G	N337S (AAC→AGC)	glnS →	glutamyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	707,102	0	A	G	97.8%	152.4 / ‑4.0	46	N337S (AAC→AGC)	glnS	glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (20/25);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

CGCGTATGCCGACCATTTCCGGTCTGCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAACGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAATGGTTACCATGCCGAACCATCCGAACAAACCGGAAATGGGCAGCCGTCAGGTGCCGTTTAGCGGTG  >  NC_000913/706958‑707245                                                                                                                                                 |                                                                                                                                                cgcgTATGCCGACCATTTCCGGTCTGCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAACGaaa                                                                                                                                            >  2:89126/1‑150 (MQ=255)  gcgTATGCCGACCATTTCCGGTCTGCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAAt                                                                                                                                          <  2:770858/151‑1 (MQ=255)    gTATGCCGACCATTTCCGGTCTGCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATgc                                                                                                                                        >  2:12201/1‑151 (MQ=255)      aTGCCGACCATTTCCGGTCTGCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATgcg                                                                                                                                       <  2:18169/150‑1 (MQ=255)       tGCCGACCATTTCCGGTCTGCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGcgcc                                                                                                                                     >  1:291848/1‑151 (MQ=255)        gCCGACCATTTCCGGTCTGCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGcgcc                                                                                                                                     <  2:543433/150‑1 (MQ=255)          cGACCATTTCCAGTATGCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCAGCGGCGTGACCAAGCAGGATAACTCCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGcgccta                                                                                                                                   <  1:336109/150‑3 (MQ=255)            aCCATTTCCGGTCTGCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCcgcgcg                                                                                                                                >  2:50996/1‑151 (MQ=255)             ccATTTCCGGTCTGCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCgcgcgc                                                                                                                               >  1:533081/1‑151 (MQ=255)               aTTTCCGGTCTGCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCaa                                                                                                                             <  1:333513/151‑1 (MQ=255)                         tGCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGtt                                                                                                                     >  1:563310/1‑149 (MQ=255)                          gCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTAt                                                                                                                   <  2:321362/150‑1 (MQ=255)                           cgtcgtcgtGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCg                                                                                                                 <  1:104899/151‑1 (MQ=255)                                        acacTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACt                                                                                                    <  2:8314/151‑1 (MQ=255)                                         cacTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGgaa                                                                                                      >  1:17733/1‑148 (MQ=255)                                          acTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTgg                                                                                                  <  2:351557/151‑1 (MQ=255)                                               ggCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATTCGGGGAAAATGGTTTTc                                                                                             >  2:126668/1‑151 (MQ=255)                                               ggCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAAATGGTTAt                                                                                              >  1:154750/1‑150 (MQ=255)                                                  ttCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGGGAAACTGGGTAt                                                                                              >  2:675497/1‑147 (MQ=255)                                                      aTTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGaaaa                                                                                        <  1:12173/149‑1 (MQ=255)                                                       ttCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACt                                                                                      <  1:305717/150‑1 (MQ=255)                                                              gTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGggc                                                                              <  2:77572/151‑1 (MQ=255)                                                                     aaaCGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTTTCGAAAACTAACAGggcgaagg                                                                         >  2:211340/1‑149 (MQ=255)                                                                      aaCGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGGTATCGAACCGGGGGAACTGGGTATTGAAAAATATCAGggggaagggg                                                                       >  1:513171/1‑148 (MQ=255)                                                                         gCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAAt                                                                   <  1:284974/151‑1 (MQ=255)                                                                          cATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAATg                                                                  <  2:324662/151‑1 (MQ=255)                                                                              ggCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACCGGGTATCGAAAACTATCAGGGCGGAGGCGAAATGGtt                                                               >  2:350338/1‑150 (MQ=255)                                                                                   gACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAATGGTTACCAt                                                          <  1:270412/150‑1 (MQ=255)                                                                                    aCCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGAACCGGTGGAACTGGTTATCGAAAAAAATTAGGGCGAAGGGGAAATGGTTACCa                                                           >  2:178695/1‑148 (MQ=255)                                                                                     ccAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGGGAAACTGGTTATCGAAAACTATCAGGGGGAAGGCGAAATGGTTACCaatcc                                                       >  2:404645/1‑147 (MQ=255)                                                                                             gACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAATGGTTACCATGCCGAACCATc                                               <  1:90715/151‑1 (MQ=255)                                                                                                   acCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAATGGTTACCATGCCGAACCATCCGAAc                                          <  2:610538/150‑1 (MQ=255)                                                                                                    ccATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAATGGTTACCATGCCGAACCATCCGAACa                                         <  2:635014/150‑1 (MQ=255)                                                                                                      aTTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTTTTGGTCCGGGGAAACTGGGTATTGAAAAATATCAGGGCGAAGGGGGAATGGGTAACATTGCGAACCATCCCAACaaa                                       >  1:334690/1‑150 (MQ=255)                                                                                                        tGAGATGGCGTGCCTGGAATCCTGCATCCGTGAAGACCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATACGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGCAATGGTTACCATGCCGAACCATCCGAACAAAcc                                     <  1:17515/150‑1 (MQ=255)                                                                                                         gagaTGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGGAGGCGAAATGGTTACCATGCCGAACCATCCGAACAAACCgg                                   <  2:391247/151‑1 (MQ=255)                                                                                                          agaTGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGGTATCGAAAACTATCAGGGCGAAGGCGAAATGGTTACCATGCCGAACCATCCGAACAAACCgg                                   >  2:627236/1‑150 (MQ=255)                                                                                                               gCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAATGGTTACCATGCCGAACCATCCGAACAAACCGGAAATgg                             <  2:433208/151‑1 (MQ=255)                                                                                                                 gTCGCCGGAATCCTGCATCCGTGGAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAATGGTTACCATGCCGAACCATCCGAACAAACCGGAAATGGGc                           <  1:405031/151‑1 (MQ=255)                                                                                                                               gCATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAATGGTTACCATGCCGAACCATCCGAACAAACCGGAAATGGGCAGCCGTCAGGTGcc             <  2:214310/151‑1 (MQ=255)                                                                                                                                cATCCGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTCTCGAACCCGTGAAACTGGTTTTTGAAAAAAATCAGGGCGAAGGCGCAATGGTTACCATACCGAAACCTCCGAAACAACCGGAAAGGGGCAGCCGTCAggaggc             >  2:411105/1‑146 (MQ=255)                                                                                                                                   ccGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAAATATCAGGGCGAAGGCGAAATGGGTACCATGCCGAACCATCCGAACAAACCGGAAAAGGGCAGCCGGCCGGTGCCGttt         >  2:333099/1‑151 (MQ=255)                                                                                                                                    cGTGAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGGGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAATGGTTACCATGCCGAAACATCCGAACAAACCGGAAAAGGGCAGCCGTCAGGTGCCGgtt         >  2:574110/1‑150 (MQ=255)                                                                                                                                       gAAGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAATGGTTACCATGCCGAACCATCCGAACAAACCGGAAATGGGCAGCCGTCAGGTGCCGTTTAGc      <  2:363625/150‑1 (MQ=255)                                                                                                                                        aaGATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAATGGTTACCATGCCGAACCATCCGAACAAACCGGAAATGGGCAGCCGTCAGGTGCCGTTTAGCgg    <  2:386343/151‑1 (MQ=255)                                                                                                                                          gATCTCAGCGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGGGAAAATGGTTATCGAAAACTATCAGGGCGAAGGCGAAATGGGTACCATGCCGAACCATCCGAACAAAACGGGAATGGGCAGCCGTCAGGTGCCGTTTAGCgggg  >  2:760186/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                CGCGTATGCCGACCATTTCCGGTCTGCGTCGTCGTGGTTACACTGCGGCTTCTATTCGTGAGTTCTGCAAACGCATCGGCGTGACCAAGCAGGACAACACCATTGAGATGGCGTCGCTGGAATCCTGCATCCGTGAAGATCTCAACGAAAATGCGCCGCGCGCAATGGCGGTTATCGATCCGGTGAAACTGGTTATCGAAAACTATCAGGGCGAAGGCGAAATGGTTACCATGCCGAACCATCCGAACAAACCGGAAATGGGCAGCCGTCAGGTGCCGTTTAGCGGTG  >  NC_000913/706958‑707245

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