Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,799,169	T→G	10.3%	intergenic (‑227/+57)	pheS ← / ← pheM	phenylalanine tRNA synthetase, alpha subunit/phenylalanyl‑tRNA synthetase operon leader peptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,799,169	0	T	G	10.3%	117.8 / 7.3	39	intergenic (‑227/+57)	pheS/pheM	phenylalanine tRNA synthetase, alpha subunit/phenylalanyl‑tRNA synthetase operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (16/19);  new base G (2/2);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

TTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAA  >  NC_000913/1799022‑1799315                                                                                                                                                    |                                                                                                                                                   ttGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTg                                                                                                                                                   <  1:127560/149‑1 (MQ=255)        cGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGc                                                                                                                                             >  3:223233/1‑148 (MQ=255)           gCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCCAGTGGAGGCttt                                                                                                                                          <  4:181251/148‑1 (MQ=255)            cTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCCAGTGGAGGCttttt                                                                                                                                        <  1:379551/149‑3 (MQ=255)             tgtgCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAg                                                                                                                                       <  1:74046/149‑1 (MQ=255)                   cGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCt                                                                                                                                  >  3:103833/1‑148 (MQ=255)                    gAAATTCTGACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGt                                                                                                                                <  2:351385/149‑1 (MQ=255)                    gAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGt                                                                                                                                >  2:279962/1‑149 (MQ=255)                      aaTTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCACGCGCTGttt                                                                                                                              <  2:308215/149‑1 (MQ=255)                        ttCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTc                                                                                                                            <  2:379547/149‑1 (MQ=255)                                tACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTc                                                                                                                                         >  4:183112/1‑128 (MQ=255)                                tACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTc                                                                                                                                         <  3:183112/128‑1 (MQ=255)                                    tAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctcttc                                                                                                                 <  3:125642/148‑1 (MQ=255)                                             tCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCt                                                                                                        >  3:228462/1‑148 (MQ=255)                                              ccGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTa                                                                                                       >  4:194462/1‑148 (MQ=255)                                               cGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAg                                                                                                      <  3:194462/148‑1 (MQ=255)                                                    tctcCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAAGCTTTTTATTTGTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGcctcc                                                                                                 <  4:111707/148‑1 (MQ=255)                                                           gcTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGAGAGGCTTTTTATTTTTCAAAACGCGCATAC‑AAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCa                                                                                         <  4:266424/148‑1 (MQ=255)                                                                   tgtCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTa                                                                                  >  4:39276/1‑148 (MQ=255)                                                                      ccAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTACACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATCCAGGTGCTAAAg                                                                               >  3:137715/1‑148 (MQ=255)                                                                              gCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAaag                                                                      >  1:370806/1‑149 (MQ=255)                                                                                  cGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAaagaag                                                                   <  4:132880/148‑1 (MQ=255)                                                                                            acCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATagc                                                        >  1:304987/1‑149 (MQ=255)                                                                                             ccTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATagc                                                        <  4:158667/148‑1 (MQ=255)                                                                                                              tttATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCg                                       >  3:205492/1‑148 (MQ=255)                                                                                                              tttAATTTTCACAATGCGCATACAAAAAAATCATCCACTGGGAGGATTTCAGGCGGTGGTTTCCGTTTATCTTCTCAAGCGCTAGCGCCATGGATTCAGGCGCTAAAGTAAAAAAAGAAGCGGAAACTGACAGCATTCAttgattgcat                                      >  2:343363/1‑147 (MQ=255)                                                                                                                       cAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt                              <  4:228462/148‑1 (MQ=255)                                                                                                                           aCGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTa                          >  4:26940/1‑148 (MQ=255)                                                                                                                               gcATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCtt                     <  2:370806/149‑1 (MQ=255)                                                                                                                                cATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCtt                     >  3:371671/1‑148 (MQ=255)                                                                                                                                    cAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCaaaa                <  1:130906/149‑1 (MQ=255)                                                                                                                                    cAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCaaaa                <  1:279962/149‑1 (MQ=255)                                                                                                                                      ctgaaaGCCTCCCAGTGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGa              >  2:274697/4‑149 (MQ=255)                                                                                                                                        gaaaGCCTCCGAGTGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGAc             >  3:236967/2‑148 (MQ=255)                                                                                                                                             ccTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTAt       <  2:248517/149‑1 (MQ=255)                                                                                                                                                ccACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTg     >  3:323185/1‑148 (MQ=255)                                                                                                                                                 cACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGaa   >  2:53845/1‑149 (MQ=255)                                                                                                                                                 cACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGa    <  3:386985/148‑1 (MQ=255)                                                                                                                                                   cTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGaaa  <  4:107079/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAA  >  NC_000913/1799022‑1799315

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