Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,776,715	A→T	S377C (AGT→TGT)	ydiF →	putative acetate‑CoA transferase YdiF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,776,715	0	A	T	96.8%	91.3 / ‑6.5	31	S377C (AGT→TGT)	ydiF	putative acetate‑CoA transferase YdiF
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (14/16);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CATTATTCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTACTTTCAGCGGAAAAATCGC  >  NC_000913/1776470‑1776956                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  cATTATTCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGtttt                                                                                                                                                                                                                                              <  1:168912/251‑1 (MQ=255)  attattCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTg                                                                                                                                                                                                                                             >  1:187404/1‑251 (MQ=255)  attagtCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCTAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGAGCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCTCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTTTTATTTGTGTTTTg                                                                                                                                                                                                                                             <  1:63400/251‑1 (MQ=255)             tGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGAc                                                                                                                                                                                                                                  >  1:102182/1‑251 (MQ=255)                          ggTGGGGAAGGTCGGCGTCGGTATTTCTGACGGCATTGGCCCGGTCGCGCGAGAAGAATGTTGTGCTGATGACTTTATTTTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTGTTCAGTGGATCGCCTTTGGCGCGAACTTGAATACCCGTGCCTTTCTGTATATGACGTCTCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTAGTTGTGTTTTGCTGAAGTCGTCCAGCACGGTAACg                                                                                                                                                                                                                     <  2:8590/251‑1 (MQ=255)                             ggggAATGTCGGCGTCGGTATGGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCg                                                                                                                                                                                                                  <  2:305319/251‑1 (MQ=255)                               ggAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTGTGGCGCGAACGTGAATGCCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTACTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGc                                                                                                                                                                                                                <  1:116556/251‑1 (MQ=255)                                       ggcggcGGTATTGCTGACGGCATTGGCCTGGTCGCTCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATaa                                                                                                                                                                                                        <  1:261675/251‑1 (MQ=255)                                         cGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAAt                                                                                                                                                                                                      >  1:327025/1‑251 (MQ=255)                                                 ttGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGt                                                                                                                                                                                              >  1:309939/1‑251 (MQ=255)                                                   gCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCCCAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCg                                                                                                                                                                                                                                    >  1:293644/1‑211 (MQ=255)                                                   gCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCg                                                                                                                                                                                                                                    <  2:293644/211‑1 (MQ=255)                                                           gaTTGGCCTGGGCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATgg                                                                                                                                                                                    <  1:223521/250‑1 (MQ=255)                                                                 ccTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCg                                                                                                                                                                              >  1:225757/1‑251 (MQ=255)                                                                         cccGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGAttt                                                                                                                                                                      >  1:139902/1‑251 (MQ=255)                                                                                     gTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGt                                                                                                                                                          >  1:315466/1‑251 (MQ=255)                                                                                                                 aGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACTTGAATTCCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGc                                                                                                                               <  2:102182/250‑1 (MQ=255)                                                                                                                         gtctgATTGGCTGAATTATTTCTCAGGGGATCGCCTTTGACGCGAACTTTAATACCCGTGCCATTCTGGATATGTCGTCACAGTTTGATTTTTATCTCGGTGGCGGTCTGGATTGTTGTGTTTTGTGTTTTGCTCAAGTCGACCAGCACGGTAACTTCGGCGTGTATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTa                                                                                                                      <  1:1078/247‑1 (MQ=255)                                                                                                                                               acaGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAg                                                                                                <  1:223038/251‑1 (MQ=255)                                                                                                                                                caGGGGATCGCCTTTGGCGCGCACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGAACAATCATGGGAACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAAc                                                                                                  >  2:280121/1‑248 (MQ=255)                                                                                                                                                   gggATCTCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATAGTACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTTGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAAt                                                                                            <  2:225757/251‑1 (MQ=255)                                                                                                                                                       tCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTTGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTCTTTGTGTTTTGCTGAAGTCGCACAGCACGGAAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACAGGTGGATTTATTGATATCAGAGCCACTTAGCAGAAAAACATTTTCTTCGGCAGATTAACTACGGGCAGTTTAAAACAAGAAATTAcc                                                                                        >  1:5124/1‑251 (MQ=255)                                                                                                                                                         gCCTTTTGCGCGAACGTGAATTCCCGTGCCATTCTTGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGa                                                                                      <  2:309939/251‑1 (MQ=255)                                                                                                                                                         gCCTTTGGCGCGAGCGTGACTACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGCGGCGGTATGGCTGTTTGTTATTTGTGTTTTGCTGAAGTCGCCCAGCCCGGTAACGTCGGAGTGCATAAATTCAATGGGAAAATCATGGGCACCGGTGGATTTATTCATATCAGTTCCACTTCGAGGAAAATCTTTTTCTGCGGCACATTAACTGCGTGCAGTTTACAAACAGAAATTACCgc                                                                                      >  2:118169/1‑250 (MQ=255)                                                                                                                                                                     aacggtaaTACCCGTGCCATTCTGGATATGACGACGTAGTTTGATTTTGATCACGGGGGCGGTCTGGATGTTTGTTATTTGTGTTGTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAAGCATGGGCACCGGTGGATTGATTGATATCAGTGCCACTTCGAAGAAAAACATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGCCGGCAAATTaaa                                                                          <  2:144204/245‑1 (MQ=255)                                                                                                                                                                     aaCGTGAATCCCCGTGCCATTCTGGATATGTCGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTTGATGTTTGTTATTTGTGGTTTGCTGAAGTCGACCAGCACGGTATCGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTaaa                                                                          <  2:187404/251‑1 (MQ=255)                                                                                                                                                                     aaCGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCACATTaa                                                                           >  1:176579/1‑250 (MQ=255)                                                                                                                                                                               cccGTGCCAGTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCag                                                                <  2:298193/251‑1 (MQ=255)                                                                                                                                                                                   tGCCAGTCTGGATATGACGTCCCAGGTTGATTTTTATCACGGTGGCGGTCTGGCTGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGATAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACTGAAATTACCGACGGCAAATTAAATATCGTCCaggacg                                                            >  1:139821/1‑249 (MQ=255)                                                                                                                                                                                                                                        aTGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTACTTTCAGCGGaaa        >  1:122578/1‑250 (MQ=255)                                                                                                                                                                                                                                             tgttATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTACTTTCAGCGGAAAAATCGc  >  1:269664/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  CATTATTCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTACTTTCAGCGGAAAAATCGC  >  NC_000913/1776470‑1776956

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