Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,776,100	G→A	E172K (GAA→AAA)	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,100	0	G	A	97.1%	117.8 / NA	34	E172K (GAA→AAA)	ydiF	putative acetate‑CoA transferase YdiF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/18);  minor base C (1/0);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGCATTATGTGGTCACTGGGGACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTG  >  NC_000913/1775854‑1776324                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                 tGGGCATTATGTGGTCACTGGGGACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAat                                                                                                                                                                                                                               >  2:205850/1‑250 (MQ=255)             ggTCACTGGGGACACTCGCCGCGTATTTCTGCACTCGCAGAACAAAATAACATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGCCATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGACGACCTGATTAAACTGGTCGAGTTTGATCACAAACAATATCTCTATTACa                                                                                                                                                                                                                   >  2:88678/1‑250 (MQ=255)             ggTCACTGGGGACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACaa                                                                                                                                                                                                                  <  2:377444/251‑1 (MQ=255)              gTCACTGGGGACAATCGCCGCGTATTTCTGAACTCGCAGAACACAATAAAATTATTGCTTATAACTCCCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCCCCAGCCTGGCATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACaa                                                                                                                                                                                                                  >  1:314178/1‑250 (MQ=255)                          aaTCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCa                                                                                                                                                                                                     <  2:624829/251‑1 (MQ=255)                            tCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGcc                                                                                                                                                                                                      >  1:555965/1‑248 (MQ=255)                              gccgcGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGat                                                                                                                                                                                                  >  1:559964/1‑250 (MQ=255)                                                 gCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTac                                                                                                                                                                              >  2:156328/1‑251 (MQ=255)                                                     aaCAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGCATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTaccacc                                                                                                                                                                          >  1:315338/1‑251 (MQ=255)                                                                  tattGCTGATAACTACCCACAAGGTGTACTTACACAAACCTGACGCGCCGCCGCAGCCCACCAGCCTGGTATTAGTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATTAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGGTTGATAACAAAAAATAGCTCTATTACAAAGCGATTGCGCCAGATATGGCCTTCATTCGCGCTACCACCGGCGACAGTGaa                                                                                                                                                              <  2:711316/250‑1 (MQ=255)                                                                             aCTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATAGCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACt                                                                                                                                                  <  1:486999/251‑1 (MQ=255)                                                                                           gtACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGt                                                                                                                                    >  1:439013/1‑251 (MQ=255)                                                                                                 acacAAACCTTACGCGCCGCCGCAGCCCCCCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATACACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTGAAGAAGACCTGATTAAACTGGTCGAGTGTGATAACAAAAAAGATCTCTATTACAAAGCGATTGCGCCAGATAGTGCCTGCATTCGCGCTACCACCTGCGACAGGGAAGGCTACGCCACTTTTGAAGATGAGGGGATGt                                                                                                                               <  2:654379/250‑1 (MQ=255)                                                                                                   acaAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTAtc                                                                                                                            >  2:527992/1‑251 (MQ=255)                                                                                                               cgccgccgcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTTATGTATCTCGACGCATTgg                                                                                                                <  1:183901/251‑1 (MQ=255)                                                                                                                gccgccgcAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCTCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGGGGTGATGTATCTCGACGCATTGGt                                                                                                               >  1:329759/1‑251 (MQ=255)                                                                                                                         gCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGccc                                                                                                       >  1:208920/1‑250 (MQ=255)                                                                                                                          cccACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTAGTGCCCAg                                                                                                     <  2:171294/251‑1 (MQ=255)                                                                                                                           ccaccaGCCTGGCATTATTAGAGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAGGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAgg                                                                                                    >  2:282275/1‑251 (MQ=255)                                                                                                                                 gCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGGCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGc                                                                                              <  2:614419/251‑1 (MQ=255)                                                                                                                                     ggTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGcaca                                                                                           <  1:146987/250‑1 (MQ=255)                                                                                                                                               gTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGt                                                                                <  1:476161/251‑1 (MQ=255)                                                                                                                                                  atatTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATCAAACTGGTCGAGTTTGATAACAAAAACTATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATCACGGCGGtctt                                                                             >  1:370791/1‑248 (MQ=255)                                                                                                                                                      tGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATtgtg                                                                          <  1:348461/250‑1 (MQ=255)                                                                                                                                                             gggACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTCAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGGGCACAATAACGGCGGTATTTtgatggaggg                                                                  >  2:293818/1‑246 (MQ=255)                                                                                                                                                             gggACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCa                                                                  <  1:156328/251‑1 (MQ=255)                                                                                                                                                                                   aGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGaaa                                            <  1:694694/251‑1 (MQ=255)                                                                                                                                                                                          cggcAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGc                                     <  2:165585/251‑1 (MQ=255)                                                                                                                                                                                                    aatgaaGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGACACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTGTGATGATGGAGGGGCTGGAAATGGGTAAGAGAGCCACGcgtcagccta                            >  1:357699/1‑241 (MQ=255)                                                                                                                                                                                                       gaaGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAAGGGTTAAGAAAgaaaagagggagaagaaatc                        <  1:603607/251‑20 (MQ=255)                                                                                                                                                                                                                 aagaagACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTAtt              <  1:527992/251‑1 (MQ=255)                                                                                                                                                                                                                              tAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCtg  <  1:25630/250‑1 (MQ=255)                                                                                                                                                                                                                                        gAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAAt                                                                                         >  2:564587/1‑153 (MQ=255)                                                                                                                                                                                                                                        gAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTGTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAAt                                                                                         <  1:564587/153‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                 TGGGCATTATGTGGTCACTGGGGACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTG  >  NC_000913/1775854‑1776324

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