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	97.9%	148.4 / ‑4.5	48	*58* (TAG→TAA)	ymdF	KGG family protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (21/26);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGG  >  NC_000913/1068120‑1068377                                                                                                                                       |                                                                                                                            gaCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCgg                                                                                                                         >  2:761427/1‑139 (MQ=255) gaCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCgg                                                                                                                         >  2:364952/1‑139 (MQ=255)  aCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                        <  2:826746/139‑1 (MQ=255)  aCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                        <  2:369580/139‑1 (MQ=255)  aCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                        <  1:364952/139‑1 (MQ=255)           gagCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAAt                                                                                                               <  2:1418982/139‑1 (MQ=255)                           tAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCa                                                                                                                 >  2:731338/1‑121 (MQ=255)                           tAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCa                                                                                                                 <  1:731338/121‑1 (MQ=255)                                     ggACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                        >  1:302195/1‑104 (MQ=255)                                     ggACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                        <  2:302195/104‑1 (MQ=255)                                      gACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGa                                                                                                             >  1:1250547/1‑114 (MQ=255)                                      gACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGa                                                                                                             <  2:1250547/114‑1 (MQ=255)                                      gACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTc                                                                                  >  1:1237248/1‑139 (MQ=255)                                      gACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTc                                                                                  >  1:955361/1‑139 (MQ=255)                                                cGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCAt                                                                        >  1:872924/1‑139 (MQ=255)                                                   gggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCAc                                                                     <  2:1237248/139‑1 (MQ=255)                                                                  tGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTCTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                        >  2:1471591/1‑139 (MQ=255)                                                                  tGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggt                                                      >  1:571047/1‑139 (MQ=255)                                                                    aCCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCaaa                                                                               >  1:248071/1‑112 (MQ=255)                                                                    aCCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCaaa                                                                               <  2:248071/112‑1 (MQ=255)                                                                          aGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATa                                              <  1:99705/139‑1 (MQ=255)                                                                                 aTCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtgatga                                       <  2:571047/139‑1 (MQ=255)                                                                                  tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                        >  1:1450974/1‑59 (MQ=255)                                                                                  tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                        <  2:1450974/59‑1 (MQ=255)                                                                                  tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                      <  2:457691/139‑1 (MQ=255)                                                                                  tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                      <  1:313008/139‑1 (MQ=255)                                                                                  tCTGAAGCAGGCAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                      <  2:955361/139‑1 (MQ=255)                                                                                           ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                        <  2:400296/50‑1 (MQ=255)                                                                                           ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCggg                                                                                                                        >  1:400296/1‑50 (MQ=255)                                                                                           ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                          <  1:799102/110‑1 (MQ=255)                                                                                           ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                          >  2:799102/1‑110 (MQ=255)                                                                                                        gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                        <  1:653820/99‑1 (MQ=255)                                                                                                        gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                        >  2:653820/1‑99 (MQ=255)                                                                                                        gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                <  2:1008221/139‑1 (MQ=255)                                                                                                          aaGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGcgcg              >  2:1234503/1‑139 (MQ=255)                                                                                                           agagCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                          <  2:1257073/94‑1 (MQ=255)                                                                                                           agagCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                          >  1:1257073/1‑94 (MQ=255)                                                                                                               cAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtggt                                                   >  2:1433283/1‑97 (MQ=255)                                                                                                               cAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtggt                                                   <  1:1433283/97‑1 (MQ=255)                                                                                                               cAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACa                                     >  2:800874/1‑111 (MQ=255)                                                                                                               cAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACa                                     <  1:800874/111‑1 (MQ=255)                                                                                                                  tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCgg                                                               >  1:361492/1‑82 (MQ=255)                                                                                                                  tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCgg                                                               <  2:361492/82‑1 (MQ=255)                                                                                                                  tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg           >  2:509787/1‑134 (MQ=255)                                                                                                                  tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg           <  1:509787/134‑1 (MQ=255)                                                                                                                      ggCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg  <  1:1234503/139‑1 (MQ=255)                                                                                                                               gACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                      >  1:557639/1‑94 (MQ=255)                                                                                                                               gACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                      <  2:557639/94‑1 (MQ=255)                                                                                                                                       |                                                                                                                            GACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGG  >  NC_000913/1068120‑1068377

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