Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,068,254	G→A	*58* (TAG→TAA)	ymdF →	KGG family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,068,254	0	G	A	98.1%	163.1 / NA	52	*58* (TAG→TAA)	ymdF	KGG family protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (21/30);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAG  >  NC_000913/1068120‑1068385                                                                                                                                       |                                                                                                                                    gaCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCgg                                                                                                                                 >  2:95644/1‑139 (MQ=255)  aCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                                <  2:421224/139‑1 (MQ=255)  aCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                                <  1:621449/139‑1 (MQ=255)  aCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                                <  1:223498/139‑1 (MQ=255)  aCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGAATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCggg                                                                                                                                <  2:1030639/139‑1 (MQ=255)      cgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACCGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                                >  1:801613/1‑135 (MQ=255)      cgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                                <  2:801613/135‑1 (MQ=255)      cgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTaa                                                                                                                            >  1:344180/1‑139 (MQ=255)                tCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGa                                                                                                                  <  1:168544/139‑1 (MQ=255)                  aGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAt                                                                                                                <  2:745256/139‑1 (MQ=255)                   gAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAtc                                                                                                               <  2:900310/139‑2 (MQ=255)                       cAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CAttt                                                                                                          <  1:1314275/139‑1 (MQ=255)                           tAAAAAAGGTGGACAGCACAGCGGGGGTCATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTgc                                                                                                      <  1:788853/139‑2 (MQ=255)                                     ggACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CAttt                                                                                                          <  2:64055/125‑1 (MQ=255)                                     ggACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CAttt                                                                                                          >  1:64055/1‑125 (MQ=255)                                      gACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTc                                                                                          >  1:972556/1‑139 (MQ=255)                                       acagcacagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCa                                                                                         >  2:332155/1‑139 (MQ=255)                                        cagcacagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATAAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCaa                                                                                        <  2:779956/139‑1 (MQ=255)                                           cacagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCaaaaa                                                                                     <  2:344180/139‑1 (MQ=255)                                                   gggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCAc                                                                             <  1:565689/139‑1 (MQ=255)                                                     gTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCAccg                                                                           >  1:424105/1‑139 (MQ=255)                                                                  tGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAt                                                                                                                >  2:352280/1‑91 (MQ=255)                                                                  tGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAt                                                                                                                <  1:352280/91‑1 (MQ=255)                                                                        gCAGCGCGCATCTGAAGCATGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCaa                                                                                        <  2:877341/107‑1 (MQ=255)                                                                        gCAGCGCGCATCTGAAGCATGTAAAAAAGATGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCaa                                                                                        >  1:877341/1‑107 (MQ=255)                                                                        gCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGa                                                        <  1:1258986/139‑1 (MQ=255)                                                                          aGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                                >  2:717549/1‑67 (MQ=255)                                                                          aGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                                <  1:717549/67‑1 (MQ=255)                                                                            cgcgcATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACt                                                    <  1:1225840/139‑1 (MQ=255)                                                                              cgcATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATc                                                                                                                          <  2:1016227/69‑1 (MQ=255)                                                                              cgcATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATc                                                                                                                          >  1:1016227/1‑69 (MQ=255)                                                                              cgcATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtga                                                  >  1:622197/1‑139 (MQ=255)                                                                              cgcATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtga                                                  <  1:95644/139‑1 (MQ=255)                                                                                  tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                              <  2:972556/139‑1 (MQ=255)                                                                                  tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                              <  1:866919/139‑1 (MQ=255)                                                                                           ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAg                                     >  2:411928/1‑139 (MQ=255)                                                                                           ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAg                                     >  2:1160406/1‑139 (MQ=255)                                                                                            gTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGc                                    >  1:722234/1‑139 (MQ=255)                                                                                                      tggtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                <  2:495191/101‑1 (MQ=255)                                                                                                      tggtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                >  1:495191/1‑101 (MQ=255)                                                                                                        gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                        >  1:748747/1‑139 (MQ=255)                                                                                                          aaGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CAGTTTTG‑CTGTAGCTCGTCAAAAGGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGcgcg                      <  2:748747/139‑1 (MQ=255)                                                                                                           agagCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGa                     <  2:424105/139‑1 (MQ=255)                                                                                                            gagCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGAt                    <  2:906818/139‑1 (MQ=255)                                                                                                               cAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCt                 >  2:364255/1‑139 (MQ=255)                                                                                                               cAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCt                 >  1:678218/1‑139 (MQ=255)                                                                                                                  tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                              <  2:1204146/107‑1 (MQ=255)                                                                                                                  tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                              >  1:1204146/1‑107 (MQ=255)                                                                                                                    aCGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCa            >  2:261635/1‑139 (MQ=255)                                                                                                                      ggCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg          <  1:1160406/139‑1 (MQ=255)                                                                                                                           aaGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTAc     <  1:411928/139‑1 (MQ=255)                                                                                                                              cGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCag  <  1:261635/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    GACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAG  >  NC_000913/1068120‑1068385

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