Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	167,798	C→A	Y105* (TAC→TAA)	fhuA →	ferrichrome outer membrane transporter/phage receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	167,798	0	C	A	96.7%	99.8 / ‑4.7	30	Y105* (TAC→TAA)	fhuA	ferrichrome outer membrane transporter/phage receptor
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (15/14);  total (16/14)
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.78e-01

GTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCA  >  NC_000913/167568‑168039                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  ttttATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTg                                                                                                                                                                                                                                <  1:391288/249‑1 (MQ=255)            ggCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCACAAAGAGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCGGTTGGTACGCGTGgcg                                                                                                                                                                                                                    >  2:16967/1‑251 (MQ=255)               aGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCAt                                                                                                                                                                                                                 <  1:402217/251‑1 (MQ=255)                  ggTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATcc                                                                                                                                                                                                               >  2:132413/1‑250 (MQ=255)                              agaagaCACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGGACGCGTGGCGCATCCAACACCTATGacc                                                                                                                                                                                                  >  2:87766/1‑251 (MQ=255)                                   acacTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCGGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTg                                                                                                                                                                                             <  1:318349/251‑1 (MQ=255)                                     acTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGAt                                                                                                                                                                                           >  1:186247/1‑251 (MQ=255)                                                      tGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCgg                                                                                                                                                                          >  2:416545/1‑251 (MQ=255)                                                                ccgcAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCa                                                                                                                                                                <  1:338540/251‑1 (MQ=255)                                                                   cAAGAAAGCGCATGGGGGCCTGCTCCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAACCCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAg                                                                                                                                                             <  2:72595/251‑1 (MQ=255)                                                                         ggcgcAGGGGGGCCTGGGGCACCTATTGCGGCGCGCCAGTCTGCTGCCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATGTCTGTTGTGACCGCCGAAGAGAGGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCGAAACGCCGGGTGTCTGTGTTGGTACGCGGGGCGCGTCCGACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGaa                                                                                                                                                       <  1:196267/250‑1 (MQ=255)                                                                              aTGGGGGCCTGCGGCAATAATGGCGGCGCGACAGTCGGCGACCGGCACTAAAACCGATACGCCGAATCAAAAATGGCCACAGTCTATTTCTGTTGGGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACt                                                                                                                                                  <  2:369811/251‑1 (MQ=255)                                                                                                                    cTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAAc                                                                                                            >  1:358656/1‑251 (MQ=255)                                                                                                                      aCCGGCCCTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACg                                                                                                           <  1:87766/250‑1 (MQ=255)                                                                                                                               aaacccGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTAACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCAt                                                                                                 <  2:432605/247‑1 (MQ=255)                                                                                                                               aaaaCCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCAt                                                                                                 >  2:46484/1‑251 (MQ=255)                                                                                                                                  aCCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTg                                                                                               >  1:158286/1‑250 (MQ=255)                                                                                                                                   ccGATACGCCGATTCAACAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGACGGACAAAGCCAGAATAACTATCTGAATGGCCTGAAGTGGCAGGGCAACTTCTATAACGATGCGGTCATTGAc                                                                                             >  2:441655/1‑251 (MQ=255)                                                                                                                                        aCGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGta                                                                                        <  1:416545/251‑1 (MQ=255)                                                                                                                                         cGCCGATTAAAAAAGTGCCACAGTCTATTCCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGtat                                                                                       <  2:368280/251‑1 (MQ=255)                                                                                                                                             gATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGc                                                                                   >  1:81060/1‑251 (MQ=255)                                                                                                                                                 cAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTgg                                                                                >  2:91791/1‑250 (MQ=255)                                                                                                                                                                         gttgtGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGccc                                                       >  1:65236/1‑251 (MQ=255)                                                                                                                                                                             tgACCGCCGAAGAGATGGCGCTGCATCAGCCGCAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGtt                                                   >  1:398450/1‑251 (MQ=255)                                                                                                                                                                                     gAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCt                                           >  2:310719/1‑251 (MQ=255)                                                                                                                                                                                          gaTGGCGCTGCCTCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGGTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTTCGGCAGAAGGTCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGTCCAGGTTTCCGTGCTTTACg                                      <  2:346287/251‑1 (MQ=255)                                                                                                                                                                                           aTGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACgg                                     >  2:151534/1‑251 (MQ=255)                                                                                                                                                                                                                      aaaGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGaa          <  2:424507/251‑1 (MQ=255)                                                                                                                                                                                                                              gcTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCTTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTc   >  2:74465/1‑250 (MQ=255)                                                                                                                                                                                                                              gcTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGtaa  <  2:53995/251‑3 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  GTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCA  >  NC_000913/167568‑168039

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