Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,989	C→G	S213S (TCC→TCG)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,989	0	C	G	96.4%	74.7 / ‑6.4	28	S213S (TCC→TCG)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (2/25);  minor base T (1/0);  total (3/25)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTG  >  NC_000913/71743‑72224                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             aaTGCTGGATCTGCCGGTGGCAAAAGTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCGGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTcgcgcg                                                                                                                                                                                                                                          <  2:240090/251‑1 (MQ=255)             gccggtgatacaaTTCATTACGCCGAATATTATCGGCTGCCGGCTGTGGCCACCGGTTTACTTCCTGCCGGGGATACTGGCGGGCGCGGCGATCGATATTCCAGGCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTGTGGCGCAGCGGTAAAGCACCCGATCGTTTAAGTCATTGTTTGTCGCGCGGCCG‑Attgtt                                                                                                                                                                                                                               <  2:268400/238‑1 (MQ=255)                       aaaTTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTTTGGCGCAGCGGGAAAGCAACCCATCGTTTTAGTCATTATTTTTTGCGCGGGCG‑TTTTTTGTTGCttgcgcc                                                                                                                                                                                                                   >  1:217954/1‑245 (MQ=255)                                  gCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGGGCAGCGGGAAAGGAACCGATCGTTTTAGGCATTATTTTTCTCGCGGCCGGTTTGT‑GTGGGTGACGCCggtgttttttg                                                                                                                                                                                                        >  1:196775/1‑241 (MQ=255)                                                   gcctgcTGTGGCCACCGTTTTTCTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGgc                                                                                                                                                                                        <  2:211933/250‑1 (MQ=255)                                                          gtgGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtggt                                                                                                                                                                                  <  2:144728/249‑1 (MQ=255)                                                              ccaccgTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGtt                                                                                                                                                                            <  2:85772/248‑1 (MQ=255)                                                               caccgTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTa                                                                                                                                                                           <  2:197781/249‑1 (MQ=255)                                                                                   gggATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATg                                                                                                                                                         <  2:196775/249‑1 (MQ=255)                                                                                    ggATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCg                                                                                                                                                      <  2:165959/251‑1 (MQ=255)                                                                                    ggATTCTGGCGGGCGCGGCGATCGATATGCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTGGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCg                                                                                                                                                      <  2:130313/251‑1 (MQ=255)                                                                                                gcgcGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGat                                                                                                                                          <  2:31144/251‑1 (MQ=255)                                                                                                    ggcgATCGGTATGCCGGGCGGAATGCAGAGCGGCGGGTTTAAATGGTTGCTGCTGGCGACGGCGGTAGTTGTGTGGGTTGGGGGGTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTGTAAGTCATTATTTGGCGCGCGGCCG‑TTTGTTGTGGCTGACTCCGTTGATTTCTGCCATCGGCGTGGTGTCGCTGGTGGTTTTAATTCGCCGCTCGTTGATGCCGGTGTATATTGATAttt                                                                                                                                      <  2:44097/251‑1 (MQ=255)                                                                                                      cgatcgatATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTg                                                                                                                                    <  1:257885/251‑1 (MQ=255)                                                                                                              atTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGt                                                                                                                            <  1:139366/251‑1 (MQ=255)                                                                                                                 ccagCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGAGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTTTCGCGCGGCCG‑TTTGTTGTGGCTGACGGCGTTGATGTCTGCCAGCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTTGTGCCGGTGTATATCGATATTTTGCGTAAAGTGGt                                                                                                                         <  1:241092/248‑1 (MQ=255)                                                                                                                                 gcggcgAGTTTAAATGGTGGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGAtt                                                                                                         <  1:71777/251‑2 (MQ=255)                                                                                                                                                        gtggcgacggCGGTATTTGTGTGGGTTGGGGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGGTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATAGTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGc                                                                                  <  2:23784/242‑1 (MQ=255)                                                                                                                                                                   ggtatttgtgtgGGTTGGTGGCTGGCTTTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTTCCTGATGCGACGCtgaccg                                                                       <  2:235612/243‑6 (MQ=255)                                                                                                                                                                                      ggctggctGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACaa                                                    <  2:213352/251‑1 (MQ=255)                                                                                                                                                                                       gctggctgTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACaaa                                                   <  1:83658/251‑1 (MQ=255)                                                                                                                                                                                              gtgCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACt                                            <  2:79230/251‑1 (MQ=255)                                                                                                                                                                                                   ggCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACtgtc                                         <  2:156982/249‑4 (MQ=255)                                                                                                                                                                                                      ggTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAg                                    <  1:116095/251‑1 (MQ=255)                                                                                                                                                                                                                       aaagcaaccgatcGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCa                   <  1:111874/239‑1 (MQ=255)                                                                                                                                                                                                                                  tcGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCTCCCGTTGATGCCGGGGTATATCGATATTTTGCGTAAAGTGGTTGGGTTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATTCGACGCTGACCGCGTCTTATCAGGCCTACAAATCGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATct         <  1:44153/249‑2 (MQ=255)                                                                                                                                                                                                                                        aaGTCATTATTTGTCGGGCGGCTG‑TTTGTTGTGGCTGCCGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCCCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAACCGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTg  <  2:4450/250‑1 (MQ=255)                                                                                                                                                                                                                                        aaGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTTTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCCTAGGTCGGATAAGGTGTTCACCCCCCCTCTGACAGTg  >  2:39902/2‑251 (MQ=255)                                                                                                                                                                                                                                        aaGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTg  <  2:66891/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             AATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTG  >  NC_000913/71743‑72224

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