Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	5,806	G→A	F218F (TTC→TTT)	yaaA ←	DNA binding and peroxide stress response protein YaaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,806	0	G	A	94.1%	88.9 / ‑1.7	34	F218F (TTC→TTT)	yaaA	DNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (17/15);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

GATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATGATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATC  >  NC_000913/5569‑6036                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       gATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACtcatc                                                                                                                                                                                                                            <  1:191856/250‑1 (MQ=255)              cgcgTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCGGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCttt                                                                                                                                                                                                              >  1:162992/1‑251 (MQ=255)                        cAGGCCTACGGGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAg                                                                                                                                                                                                    <  1:111620/251‑1 (MQ=255)                         aGGCCTACGTGAATTCTGCAATATATTGAATCTGCCTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGt                                                                                                                                                                                                      >  1:248633/1‑248 (MQ=255)                          agcatACGTGAAATCTGCACTATATTGCATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGAATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACACGTTCGCCATTGCTGGAGGAAGCTTCATAAAAGAAGAAACCTTAGATATTAAAACAAGTAAGTAGCACTGGATTAGTAAGCCGATTTTCCATAATACACCGACTCATCAAACCGCGGGCTTTCTAAGCATAGa                                                                                                                                                                                                   >  2:233712/5‑250 (MQ=255)                                  tgaaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGC‑CGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCtgattat                                                                                                                                                                                         >  1:299986/3‑248 (MQ=255)                                           cAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTATCGTAGAAGATGATTATCTTAAAtt                                                                                                                                                                                 >  1:296708/1‑251 (MQ=255)                                                 atTGAATCTACATGCTTTGGTAGGACGGATACGGCGGTCGCGCCGCATCAGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACCCCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGATGGGTCAGCCGATTTTCAATAATAAAACGCCTCATCAGACCGCGTGCGTTATTAGCGTCGAAgatgattacatataatttgcagt                                                                                                                                                                           >  2:272168/1‑234 (MQ=255)                                                   tGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGttc                                                                                                                                                                         >  1:255503/1‑251 (MQ=255)                                                               cTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAgg                                                                                                                                                             >  1:204648/1‑251 (MQ=255)                                                                ttttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGCCTGCAAACTTAACGCTGCTCGTAGCGATTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCtt                                                                                                                                                                                                               >  1:328368/1‑200 (MQ=255)                                                                           ggATAAGGCTTTCAAGCCGCAACCGGCATTGACTGCAAACTTAACGCTGCGCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCGGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGGAGAAGCTGATGATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTg                                                                                                                                                 <  2:162992/251‑1 (MQ=255)                                                                                       cACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTAAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGACGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTGGGc                                                                                                                                      >  1:181504/1‑250 (MQ=255)                                                                                               aTCCGGCATTGACTGCAAACCTAACGCTGCACGGAGAGGTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTTATCAAAGAAGTAACCGTCGATATTAAAACCAGTCATTTGCGCTTGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCATACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAAttt                                                                                                                             <  2:213523/251‑1 (MQ=255)                                                                                                      aTTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGc                                                                                                                      >  1:14576/1‑251 (MQ=255)                                                                                                         gACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAAGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTc                                                                                                                   <  1:274553/251‑1 (MQ=255)                                                                                                              cAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGa                                                                                                              <  1:78423/251‑1 (MQ=255)                                                                                                                    tAACGCTGCTCGTAGCGATTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGAACGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGCTTTaa                                                                                                         >  1:169723/1‑250 (MQ=255)                                                                                                                           gctCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCAtt                                                                                                                                   <  1:100137/217‑1 (MQ=255)                                                                                                                                  gCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGGTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACg                                                                                          <  1:140239/251‑1 (MQ=255)                                                                                                                                         aaCACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGCTTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGAtt                                                                                   >  1:252386/1‑251 (MQ=255)                                                                                                                                                             ggaggtagcctcatctAAGAAATAACCTTCGCTATAAAAACCAGTCAGTTGCTCTGGGTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGAGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTGGCCGTTCTTCTCATCGATGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCt                                                                <  2:377747/240‑1 (MQ=255)                                                                                                                                                               aggaAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTtgt                                                             <  1:227727/251‑1 (MQ=255)                                                                                                                                                                              agaagTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCg                                              <  1:342192/251‑1 (MQ=255)                                                                                                                                                                              agaagTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCg                                              >  1:104307/1‑251 (MQ=255)                                                                                                                                                                                aagTAACCTTCGCTAGTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCGTCACCGATTTAAAATACTCGGCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGtt                                            <  1:218572/251‑1 (MQ=255)                                                                                                                                                                                           ccTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCATCCGATTTGCAATAATAAAACGACTCATCAGACCGCGTTCTTTCTTAGCGTAGAATCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGTATCACCGGCTTGATAATCTCGGCATTCAATTTCCTCGGCTTCTCCGATTTAAAATACTCATCTGACGCCATATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTgg                                 <  2:14576/250‑1 (MQ=255)                                                                                                                                                                                              aTTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGtg                               >  1:70530/1‑250 (MQ=255)                                                                                                                                                                                                   aaCCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTgatgat                          >  1:88161/1‑250 (MQ=255)                                                                                                                                                                                                       cgtcaggTGCTCTGGTTTGGTCTGCCGATTTTAAATAATATAACGACTCTGCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCTTTCTTCCCATCGAGGAACACCGGCTTGATAATCTCGGCATTCATTTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCTCCACCATATCGCCGTTTTTTGCGAGCGCCTCGTTCAGCTTGTTGGTTATGATAtctc                     <  2:255503/250‑1 (MQ=255)                                                                                                                                                                                                                   gggTTTGGGCAGCCGATTTGTAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTCAATGTGCCGTTTCTCGCATCGAGGAACACCGGCTTTATAATCTCGGCATGCAATTTCTTCGGCTTCTCCGATTCAAAAGACTCATCTTACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGTCTCGTTCATCTTGTTGGTGATGATATCTCCCCAGAATTGa          <  2:188050/249‑1 (MQ=255)                                                                                                                                                                                                                    ggTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATATTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATATCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACACCAGCTTAATCACCACATTAACGCCTTGTTCTGCGAACGCCTCGTTCAGCTTGTTTGTGATGATATCAACCCCGAATTGAtc        >  2:311906/1‑250 (MQ=255)                                                                                                                                                                                                                                     ttCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATc  <  2:191932/240‑1 (MQ=255)                                                                                                                                                                                                                                     ttCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATc  >  1:191932/1‑240 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       GATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATGATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATC  >  NC_000913/5569‑6036

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