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	96.7%	85.6 / ‑4.3	30	D202A (GAT→GCT)	uxaB	tagaturonate reductase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (12/17);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

TTGCCACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTT  >  NC_000913/1609831‑1610323                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        ttGCCAATTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGAGCTAACCCAGTTCTTCTTCCAGGTTAGCGACTTCATCGCGCGGAGAACCGGTAACGATACGATCCACCCGCGTAGAACAGAAGCAGTTAGCCTGAGCCAGc                                                                                                                                                                                                                                                    >  2:176378/1‑251 (MQ=255)   gCCACTTTGCGTTCTTTGTACGGTTTAATATCGGCAACAATCAGCACGTTGAGCGGATATTGGTCCAGACGCAGGTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTGCTGCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAgccagcca                                                                                                                                                                                                                                                  <  2:218859/251‑1 (MQ=255)       cTTTGCGTTCTTTGTACGGTTTACTATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTCAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCCTTg                                                                                                                                                                                                                                               >  2:169511/1‑250 (MQ=255)          tGCGTTCTTTGTACGGTTTAATATCGACAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGAATTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAAt                                                                                                                                                                                                                                            >  1:208125/1‑250 (MQ=255)                         gTTTAATATCGTCAAAAAACAGCACGTTGAGCGGATATTTGTCCAGAAGCAGTTCAGTAGCTAATGATTTCGTTCCCGGAATCACAAACAGGTAAAAGTGTGCAGCGGTGTCGAGAAAACCATCTTGATAACCCAGTTCTTCTTCGAGTTGAGCCACGTCAGCGTGCGGATAGCCGGTAACGATACGGTTCAACAGCGTATAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGt                                                                                                                                                                                                                            <  2:208125/251‑1 (MQ=255)                                        caaTCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCa                                                                                                                                                                                                             <  1:292127/251‑1 (MQ=255)                                         aaTCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCTTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCAt                                                                                                                                                                                                            <  1:45065/251‑1 (MQ=255)                                            cAGCCCGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCACCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAgc                                                                                                                                                                                                         <  2:242744/251‑1 (MQ=255)                                                             atTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCACTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACg                                                                                                                                                                                         <  2:31227/250‑1 (MQ=255)                                                              tttGTCCAGAAGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCAAGTTTAGCCACTTCATCGCGCGGATAACCGGAAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGCGCCGGCATTTGAATAAATGCTTCTGGCAGTGCCCATTCTTCGGAAAAGCGCAGCACCAGTTCACGc                                                                                                                                                                                        >  2:384266/1‑250 (MQ=255)                                                               ttGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCa                                                                                                                                                                                       >  2:255562/1‑250 (MQ=255)                                                                 gTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAAc                                                                                                                                                                                     <  2:62272/250‑1 (MQ=255)                                                                                             ttCGGTCCCTGAAGTACAAACAGGTAAAAGTGTTCAGCGTTGTCAAGAAAACCGTCGTGATAACCCAGTTCGTCTTCCAGTTAAGCCCCGTCATCGCGCGGATAACCGGTAACGATACGATCCACCAGCGTAGAACAGAGGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCACCCTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTATTCAATCAACTcac                                                                                                                                                         <  2:305872/250‑1 (MQ=255)                                                                                                      tGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAatgat                                                                                                                                               >  1:95750/1‑251 (MQ=255)                                                                                                              aaaCAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAAc                                                                                                                                       <  2:275702/251‑1 (MQ=255)                                                                                                                    gTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTAt                                                                                                                                 >  1:3622/1‑251 (MQ=255)                                                                                                                      aaaaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATcc                                                                                                                               <  1:108800/251‑1 (MQ=255)                                                                                                                          gtgtTCAGCGGTGTCGATAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGTATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGTTGCACCAGTTCACGCAACGCGTCGCCATTTTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAgcg                                                                                                                           <  1:301128/251‑1 (MQ=255)                                                                                                                             ttCAGCGGAGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGCCAATCAACTCACACGGAATGCTGATCCAACCTTTATCCAGCGCAc                                                                                                                        >  2:347476/1‑251 (MQ=255)                                                                                                                              tCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCAc                                                                                                                        <  1:266096/250‑1 (MQ=255)                                                                                                                               cAGCGGTGCCTAGAAAACCATCGGGATAACCCAGTTCTTCTTCCATTTGAGCCACTTCATCGCGCTGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTCTGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCAct                                                                                                                       <  2:130278/250‑2 (MQ=255)                                                                                                                                       gcGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTGTATCCAGCGCACCGTTGAAATg                                                                                                              <  1:357625/250‑1 (MQ=255)                                                                                                                                                             ccAGTTCTTCTTCCCGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCCACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAAcagc                                                                                           >  1:61552/1‑248 (MQ=255)                                                                                                                                                                 ttcttcatcCAGTTTAGCCACTTCATCGCGCGGATAACCGGAAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCCATGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTTTAGTCAATCAACTCACACGGAATGATGCTCCAACCTTTATCCAGAGCACCGTTGAACTGGCTAAAGCGTTCGAACAGCAGACggg                                                                                    >  1:65928/1‑251 (MQ=255)                                                                                                                                                                      cttcCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGACAGt                                                                               >  1:296161/1‑251 (MQ=255)                                                                                                                                                                                         aTCGCGCGGATAACCGAGAACGAGACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTGATCCAGCGCACCGGTGAAATGTCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCg                                                            <  2:2153/251‑1 (MQ=255)                                                                                                                                                                                                ggATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCACCTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGgcgc                                                     >  2:341758/1‑251 (MQ=255)                                                                                                                                                                                                    aaCCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCATCTCACACGGAATGATGATCCAACCTATATCCATCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTc                                                  <  1:384266/250‑1 (MQ=255)                                                                                                                                                                                                                                     aGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCGCTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCAGTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGGTGAAATGGCTAAAGCGTGCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTATCTGGTAccc                <  2:408072/251‑1 (MQ=255)                                                                                                                                                                                                                                            gTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTTCCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACTCGCTTc           <  1:347476/249‑1 (MQ=255)                                                                                                                                                                                                                                                     aTCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGACATTGTAGTCAATCAACTCCCACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGATAAAGCGTTCGAACAGCCGACGGGTCAGTTTTGACGGATAGCTTAACGCTCGCGCGTCCTCGAATTTATCGCACGCATGGTAGCTGATACCCGCTTCGGTGGTGtt  >  2:314752/1‑249 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        TTGCCACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTT  >  NC_000913/1609831‑1610323

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