Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,610,076	T→G	D202A (GAT→GCT)	uxaB ←	tagaturonate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,610,076	0	T	G	90.3%	67.5 / ‑1.0	31	D202A (GAT→GCT)	uxaB	tagaturonate reductase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (14/14);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

GATTGCCACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTG  >  NC_000913/1609829‑1610321                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      gATTGCCACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCa                                                                                                                                                                                                                                                    <  2:209364/251‑1 (MQ=255)       cACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCCGCACGTTGAGCGGCTATTTGTCCAGACGCAGTTCAGACGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGACACTTCATAGCGCGGATAACCGGTAACGATCCGGTCCAACAGCGTAACACAGAAGCTGTTAGCCTGAGCCAGCCAtt                                                                                                                                                                                                                                              >  2:110160/1‑251 (MQ=255)            tGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGCTATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGGCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGTGAAAACCATCGTGATAACCCAGTTCTTCTTACAGTTTAGCCACTTCATCGCGCGGCTAACCGGTAACGATACGGTCCACCAGCGTAGAACAGACGCTGTTAGCCTGAGCCAGCCATTGAAt                                                                                                                                                                                                                                          >  2:89535/1‑250 (MQ=255)                 tCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATg                                                                                                                                                                                                                                     <  1:16077/250‑1 (MQ=255)                                     gctaacaaTCAGCACGTTGAGAGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAATGTGTTCAGCGGTGTCGAGAAAACCAGCTGGATAACCCCGTCCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATACCCGGTAACTATACGGTCCACCAGCTTAGAACCGAATCTGTGAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCtt                                                                                                                                                                                                                <  2:44600/248‑1 (MQ=255)                                     gTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCtt                                                                                                                                                                                                                <  1:141541/251‑1 (MQ=255)                                            aTCAGCACGTTGAGCGGATATTTGTTCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATa                                                                                                                                                                                                         <  1:17869/251‑1 (MQ=255)                                                             atatTTGTCCAGACGCAGTTCAGTCGCTAAGGCTTTCGGTCACTGAATCACAAACAGGTACAAGTGTTCAGCGGTGTCGAGACAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGACACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCAAGCCATTGCATAAATGCTTCTGGCAATGCCCACTTTTCCGCATCACGCAGCACCCGTTCAc                                                                                                                                                                                        >  2:392406/1‑251 (MQ=255)                                                                            cAGTTCAGTCGCTAAGGATTTCGGCCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCAGTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACTCAACGCGGCGCCAt                                                                                                                                                                         <  2:67021/251‑1 (MQ=255)                                                                                  aGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCCTTGTAGt                                                                                                                                                                   >  1:203192/1‑251 (MQ=255)                                                                                       cTAAGGATTTCGGTCCCGGAATCACAAGCTGGTACAAGTTTTCAGCGGTTTCGAGAAAGCCATCGAGATAACCCAGTTCTTCTACCCGTTTAGCCACTTCATCGCGCGGATAACAGGTCACGATACGGTCCAACATCGTCGAACAGAAGATGTTGGCCTCAGCCGGCTATTGAATAAATGCTTCTGGCAGTGCCCACTCTAGCGCATAGCGCAGCACCAGTTCACGccaccatgtgcccttgtagtccatc                                                                                                                                                              >  2:67864/1‑227 (MQ=255)                                                                                         aaGGATGTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGtcaatca                                                                                                                                                             <  1:174509/250‑1 (MQ=255)                                                                                                             cacaAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATcc                                                                                                                                        >  1:108620/1‑251 (MQ=255)                                                                                                              acaAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCa                                                                                                                                       <  1:303509/251‑1 (MQ=255)                                                                                                                    aGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCttt                                                                                                                                 >  2:116552/1‑251 (MQ=255)                                                                                                                         aaaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATcc                                                                                                                             >  1:361140/1‑250 (MQ=255)                                                                                                                                             gaAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTa                                                                                                        <  1:60341/251‑1 (MQ=255)                                                                                                                                                               ccAGTTCTTCTTCCAGATTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGTGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAAcagc                                                                                         >  1:377077/1‑248 (MQ=255)                                                                                                                                                                cAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGGGCCATCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTTGTCAATCAACTCACACGGAATGATTATCCAACCTTTATCCAGCGCACCTTCGAAATGGCTAAAGCGTTCGAACAGCAGAc                                                                                     <  2:405060/251‑1 (MQ=255)                                                                                                                                                                      ttctgcCAGTTTAGCCCCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTc                                                                                <  2:466364/250‑1 (MQ=255)                                                                                                                                                                           ccAGTTTAGCCACGTCATCGCGCGTATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGACGCTGTTAGCCTGCGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCACCGCGTCACCCTTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTACAGCGTACGCTCACCAGCCGGGTCAGtttt                                                                          >  2:74558/1‑251 (MQ=255)                                                                                                                                                                              gTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGcc                                                                       >  1:227330/1‑251 (MQ=255)                                                                                                                                                                                             cgcgcgGATAACCGGTATCGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCATGCGCATAGCGCAGCACCAGGTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCTCTCGGAATGATGATCCAACCTTCATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCTGATAGCTTACCGCt                                                        <  2:437371/251‑1 (MQ=255)                                                                                                                                                                                                                gATACGGTCCACCAGCGTAGAACAGAAGATGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTAt                                     >  1:157730/1‑251 (MQ=255)                                                                                                                                                                                                                      gTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCAAACCTTTATCCAGCGCACCGCTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCCACGACTTTATCGCCCg                               >  2:25647/1‑251 (MQ=255)                                                                                                                                                                                                                                    aaCAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGAt                  >  1:355431/1‑250 (MQ=255)                                                                                                                                                                                                                                     aCAGAAGCTGTTCGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTGCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCTCACCGTGGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATAc                <  1:300957/251‑1 (MQ=255)                                                                                                                                                                                                                                           gCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCtt          >  1:2632/1‑251 (MQ=255)                                                                                                                                                                                                                                               ttcgCCTTATCCCGCCATTGAATAAATGCTTCTGGCAGTGCCCCCTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTTGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGAt                  <  1:134366/236‑1 (MQ=255)                                                                                                                                                                                                                                               ttAGCCTGATCCCGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGAt                  >  2:134366/1‑239 (MQ=255)                                                                                                                                                                                                                                                   ccTGAGCCAGCCATTTAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGgtggtg  <  2:361140/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      GATTGCCACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTG  >  NC_000913/1609829‑1610321

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