Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,629,247	C→T	I11I (ATC→ATT)	ydfZ →	putative selenoprotein YdfZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,629,247	0	C	T	96.9%	86.8 / ‑4.7	32	I11I (ATC→ATT)	ydfZ	putative selenoprotein YdfZ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (20/11);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCT  >  NC_000913/1629014‑1629483                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             cAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCgtgt                                                                                                                                                                                                                             >  1:138879/1‑251 (MQ=255)    cGACTGGTTTATCGAAGATTAATTCATCTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGTATCAAAAAAGATGTGAGCTTTATCAAAATCAAAAAATATTTCTCTCGACATGAGTATTTATACTTCGTCCGTTACGTGGGCTTCGACTTTAAATTAGAAAGGAGAAAACTCCTATGACGACCTGCGATCGGAACCGTAATGCAATTAGCACTGGCCGCCGTGTTAt                                                                                                                                                                                                                          <  2:204513/251‑1 (MQ=255)         ggTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATTGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGtt                                                                                                                                                                                                                      <  1:429564/250‑1 (MQ=255)          gTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCAGAGATATGGtt                                                                                                                                                                                                                      >  1:261557/1‑249 (MQ=255)                               ttCCTCTCCTCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCCTTCATCAAAAAAGATGTTAGCTTGATCAAAAACAAAAAATATTTCCCTCGCCAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGCCTGTAAATCCGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCCCCCTGGCa                                                                                                                                                                                                >  2:166218/1‑250 (MQ=255)                                      catcccatccGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCCTGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGAAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACAATGGCAAGATCCTg                                                                                                                                                                                        >  2:315126/1‑251 (MQ=255)                                      catcccatccGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTg                                                                                                                                                                                        >  2:386998/1‑251 (MQ=255)                                                        agTCTTTTCCCCCGCCTTATGGCGCATGCATGCAGCAAAAAAGAGGTGAGCTTGATCAAAAACAAAAAATATTTTACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTTGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCGACGATCGTAACCGTAACGCAATTACCACGGGCAGGCGTGTTATGGTTAGCGGCACATGTCACTCTGTCAAGATCCTGTCGATTGATACTGAAGGt                                                                                                                                                                      <  2:41070/251‑1 (MQ=255)                                                                            ggCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATcc                                                                                                                                                   <  1:303956/250‑1 (MQ=255)                                                                             gCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATccgc                                                                                                                                                 >  1:47311/1‑251 (MQ=255)                                                                                catgcatgCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGGTGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATccgcc                                                                                                                                                <  1:230418/249‑1 (MQ=255)                                                                                           cAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGAGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGt                                                                                                                                   >  1:344128/1‑251 (MQ=255)                                                                                                                caaaaacaaaaaATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTGTGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTAAAGGTCTGACCGATGAGCAAATCAGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGa                                                                                                                 >  1:435843/1‑248 (MQ=255)                                                                                                                    aacaaaaaATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACATACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGACGTTCTGCCCGCTGAGCAACTCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTtgtgtagcaaa                                                                                                           >  1:233130/1‑243 (MQ=255)                                                                                                                          acatatTTCACTAGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGTGAAAACACCTATGACGACCTACGATCGTAACAGTAACGCGATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCACGATCCTGTCGATTGAGACTGACGGTCTGACCGCTGCGCAAATCCGCCGCGGACAAACCGTAGTTGTTTAAGGTTGTGACGAGACACTGGcc                                                                                                    >  2:114631/3‑250 (MQ=255)                                                                                                                          aaaTATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGTCAAGATCCTGTCGTTTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTTGGAAGAGAAACTGGcc                                                                                                    >  2:192303/1‑250 (MQ=255)                                                                                                                            atatTTCACTCGACAGGATGATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGATAAAACCCATATGACGACCTACGATCGTAACCGTAACGCAATTACCCCTGGCAGCCGGGTTATTGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGGCCGCTGAGCAAAACCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGcacc                                                                                                  <  2:98443/251‑1 (MQ=255)                                                                                                                                ttCACTCGACAGGATGATTTATATTGCGCCCGTTACGTGGGCTTCGACGGTAAAGCAGAAAGGAGAAAACACCTAGGACGACCTACGATCGTAACCGTAACGCAAACACCCCTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACt                                                                                               <  2:259318/250‑1 (MQ=255)                                                                                                                                                          cgcCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACAACTGGCAGCCGTGATATGGATAGAGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGATCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGATGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTATCGGCATGAACt                                                                     >  1:466516/1‑250 (MQ=255)                                                                                                                                                             ccGTTACGTGGGCTTCGACTGTAAATCAGACAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGCTTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCATCTCGGCCTGAACTatgc                                                                 >  2:20953/1‑248 (MQ=255)                                                                                                                                                              cGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCACGATCCTGTCGCTTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAACCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCg                                                                >  2:245437/1‑251 (MQ=255)                                                                                                                                                                ttACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGCCGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGATGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCgtg                                                              >  1:478785/1‑251 (MQ=255)                                                                                                                                                                 tACGTGGGATTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCgtg                                                              <  1:81143/250‑1 (MQ=255)                                                                                                                                                                                tGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATgg                                              >  1:190241/1‑251 (MQ=255)                                                                                                                                                                                 ttAAATAAGAAAGGAAAACACCCCTATGACGCCCTCCGATCGTAACCCTAACACAATTACCACGGGCAGCCGTGGTGGGGTTGGCGGCACCGGTCACACTGGAAAGATCCTGTCGATTGATACTGAAGTGCTGACCGCTGAGCAAACCCGCCGCGTAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGGGTGAATGCCGCCGATg                                               <  2:298283/248‑1 (MQ=255)                                                                                                                                                                                       cAGAAAGGAGAAAACACCTATGACGACCGACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCAt                                        <  1:245437/250‑1 (MQ=255)                                                                                                                                                                                            aGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCttt                                  <  1:208218/251‑1 (MQ=255)                                                                                                                                                                                                  aaaCACCTATGACGACCTACGATCGTAAACGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGTTCACACTGGCAAGATCCTGTCGATTGCTCCTAAAGGCCTGACCGCTGAGCAAATCCGCCGCGGAAACACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAGGCGTGTGACTGCCGCCGGTGGCGGCTTAGCTTTTTTACt                            >  2:377043/1‑251 (MQ=255)                                                                                                                                                                                                    acacCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTc                          >  1:306753/1‑251 (MQ=255)                                                                                                                                                                                                       ccTATGACGACCTACGATCGTAACCGTAACGCCATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAGGCGTGTGAATGCCGCCGATGGCGGCCTTGCTTTTTTACTTCACg                       >  2:474252/1‑251 (MQ=255)                                                                                                                                                                                                             actacCTACGATCGGACCCTAAACGCAATTACCACTGGCAGCCGTGTTAGGGTAAGCGGCCCCGTTCACACGGGCAAGATCCGGTCGATTGATACTGAAGGTCGGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGTTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGTGTGTGAATTCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAata                  <  2:325902/247‑1 (MQ=255)                                                                                                                                                                                                                            aaccgtaacGCACTTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGTGAAACTGGCCCCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGACTGCCGCCGATGGCGCCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCt  >  2:52642/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             CAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCT  >  NC_000913/1629014‑1629483

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