Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	5,851	G→T	I203I (ATC→ATA)	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,851	0	G	T	96.9%	92.1 / ‑3.9	32	I203I (ATC→ATA)	yaaA	DNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (17/14);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

CAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATGATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATC  >  NC_000913/5611‑6096                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       cAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTATCGTAGAAGATGATTATCTTAAAtt                                                                                                                                                                                                                                              >  1:296708/1‑251 (MQ=255)         tGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGttc                                                                                                                                                                                                                                      >  1:255503/1‑251 (MQ=255)                     cTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAgg                                                                                                                                                                                                                          >  1:204648/1‑251 (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)                                                                                                                                    agaagTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCg                                                                                                          >  1:104307/1‑251 (MQ=255)                                                                                                                                    agaagTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCg                                                                                                          <  1:342192/251‑1 (MQ=255)                                                                                                                                      aagTAACCTTCGCTAGTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCGTCACCGATTTAAAATACTCGGCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCGtt                                                                                                        <  1:218572/251‑1 (MQ=255)                                                                                                                                                 ccTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCATCCGATTTGCAATAATAAAACGACTCATCAGACCGCGTTCTTTCTTAGCGTAGAATCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGTATCACCGGCTTGATAATCTCGGCATTCAATTTCCTCGGCTTCTCCGATTTAAAATACTCATCTGACGCCATATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCGTTCAGCTTGTTgg                                                                                             <  2:14576/250‑1 (MQ=255)                                                                                                                                                    aTTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCGTTCAGCTTGTTGGtg                                                                                           >  1:70530/1‑250 (MQ=255)                                                                                                                                                         aaCCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCGTTCAGCTTGTTGGTgatgat                                                                                      >  1:88161/1‑250 (MQ=255)                                                                                                                                                             cgtcaggTGCTCTGGTTTGGTCTGCCGATTTTAAATAATATAACGACTCTGCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCTTTCTTCCCATCGAGGAACACCGGCTTGATAATCTCGGCATTCATTTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCTCCACCATATCGCCGTTTTTTGCGAGCG‑CCTCGTTCAGCTTGTTGGTTATGATAtctc                                                                                 <  2:255503/250‑1 (MQ=255)                                                                                                                                                                         gggTTTGGGCAGCCGATTTGTAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTCAATGTGCCGTTTCTCGCATCGAGGAACACCGGCTTTATAATCTCGGCATGCAATTTCTTCGGCTTCTCCGATTCAAAAGACTCATCTTACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑TCTCGTTCATCTTGTTGGTGATGATATCTCCCCAGAATTGa                                                                      <  2:188050/249‑1 (MQ=255)                                                                                                                                                                          ggTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATATTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATATCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACACCAGCTTAATCACCACATTAACGCCTTGTTCTGCGAACG‑CCTCGTTCAGCTTGTTTGTGATGATATCAACCCCGAATTGAtc                                                                    >  2:311906/1‑250 (MQ=255)                                                                                                                                                                                           ttCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATc                                                              >  1:191932/1‑240 (MQ=255)                                                                                                                                                                                           ttCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATc                                                              <  2:191932/240‑1 (MQ=255)                                                                                                                                                                                                   taaaaggACTCCTCAGGCCGCGTGCTTTCTTAGCGTAGAAGCTCCTTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCATCACCGATTTAAAATACTCATCTGACCACAGATTAATCACCACATTATCTCCTTATGCTACGAGCGTCCAGATTC‑GCTAGTTGGTGACGATACCTCCCCAGTTTAAATACAGATCTTTCCCTTGAGCATtcacc                                           >  2:318170/7‑248 (MQ=255)                                                                                                                                                                                                         gACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGCTTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACgg                                     >  1:114350/1‑251 (MQ=255)                                                                                                                                                                                                                     cgcgTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTcc                         >  1:150442/1‑251 (MQ=255)                                                                                                                                                                                                                              tCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATa                 <  1:69868/250‑1 (MQ=255)                                                                                                                                                                                                                                ttAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAgg              >  1:54898/1‑251 (MQ=255)                                                                                                                                                                                                                                          acatgattaTCTTAAATGTGCCGTTCTTCTCATCGAGGAACCACGGCTTTATATTCTCGTCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCTCATTATCGCCTTGTGCTGCGAGCG‑CCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGTCGATAAGGCTGCATTaaa    <  1:92930/248‑1 (MQ=255)                                                                                                                                                                                                                                            cTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATc  >  1:281186/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       CAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATGATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCG‑CCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATC  >  NC_000913/5611‑6096

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: ‑