Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,610,440	C→T	E81K (GAA→AAA)	uxaB ←	tagaturonate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,610,440	0	C	T	96.4%	96.0 / NA	28	E81K (GAA→AAA)	uxaB	tagaturonate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (16/11);  minor base G (0/1);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTCCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTT  >  NC_000913/1610195‑1610669                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       aTCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTtcccc                                                                                                                                                                                                                                    <  1:306400/250‑1 (MQ=255) aTCCAGCGCACCGGTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCt                                                                                                                                                                                                                                   >  2:297638/1‑251 (MQ=255)       cgcACCGTTGAACTGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCAt                                                                                                                                                                                                                             >  1:88052/1‑251 (MQ=255)                           aGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAAtggtg                                                                                                                                                                                                         <  2:223358/251‑1 (MQ=255)                                      agcagACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGAcc                                                                                                                                                                                              <  1:162605/251‑1 (MQ=255)                                          gACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATcc                                                                                                                                                                                          >  1:120091/1‑251 (MQ=255)                                               gTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGt                                                                                                                                                                                     >  2:231125/1‑251 (MQ=255)                                                                      cGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGa                                                                                                                                                                >  1:74625/1‑249 (MQ=255)                                                                          ggctcgTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCATTGTAGACGCTGAGTTCACGATTAACCGATCGAATCAGACGCGCGTCGCGGACCGCTTTCCCCTGCTCATTCAGGCCGCGGATAATGGGGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTc                                                                                                                                                          <  1:36722/247‑1 (MQ=255)                                                                             gcgTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAAt                                                                                                                                                       >  1:217022/1‑251 (MQ=255)                                                                                               cccGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCAGCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCAc                                                                                                                                     <  1:206829/251‑1 (MQ=255)                                                                                                    aTGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAg                                                                                                                                                                               >  2:328396/1‑204 (MQ=255)                                                                                                    aTGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAg                                                                                                                                                                               <  1:328396/204‑1 (MQ=255)                                                                                                      ggTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTAACTGCAGACGCTGATTTCACGATTAACCGAGCGCATCAGACGCGCGTCGCTGACCGCTTTCACCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGGGCCGGCGGGAATGCAGTTTCAATCCGCCGAAAAACGACCACGACAGa                                                                                                                               >  2:334102/1‑250 (MQ=255)                                                                                                                         ggtggtGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGAACGGTGTGCTCAt                                                                                                           >  1:236807/1‑251 (MQ=255)                                                                                                                                  agagAAAACAAAGCGCATTTCCGGGTTGTGCTCCAGTTTCAGGAATGCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGGGTCGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGTTGTTCTCATTCAGgaga                                                                                                   <  2:240348/250‑1 (MQ=255)                                                                                                                                           aaaGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCG‑ATTTGcc                                                                                        >  2:160483/1‑251 (MQ=255)                                                                                                                                           aaaGCCCATGTCCGTGTTGTGAGCCAATGTCAGGAATTCATCGTATTCACTGTCGACGCTGATTGCACGAGTAACCGAGCGAATCGTACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCGGATCGGTGTGCTCATTCAGCAGATCG‑ATTTGcc                                                                                        <  2:313823/251‑1 (MQ=255)                                                                                                                                             aGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTGCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCG‑ATTTGCCAg                                                                                      <  2:66270/251‑1 (MQ=255)                                                                                                                                                    ttCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTAAGATCGGTGGGCTCATTCAGGAGATCG‑ATTTGCCAGTcaaaa                                                                                >  2:261787/1‑248 (MQ=255)                                                                                                                                                          gTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCG‑ATTTGCCAGTCAACAAAggggc                                                                         >  1:44039/1‑248 (MQ=255)                                                                                                                                                                        tCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCt                                                           >  1:184545/1‑251 (MQ=255)                                                                                                                                                                               ttCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGCGCGAAACAGACACGCGTCGCTGACCCCTTTCCCCTTCCCATTCAGGCCACGGATCATGGTGGTGTACAGACCATCCTCCCTGCTCAGTGACGGCGGGAATGACGTTTCAAAAGGACGAACAACGACCCCGCCAGCATTCAAGTCGGTGTCTTCCTTCAGCAGATCGCATTT‑CCAGTCACCAAACGCGCGCCGGAAGTTACCTTccccc                                                     >  1:20150/1‑246 (MQ=255)                                                                                                                                                                                    cGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTg                                                <  1:160483/250‑1 (MQ=255)                                                                                                                                                                                                aGACGCTGAGTTCACGATTATGCGAGCGAATCAGACGAGCGTCGCTGACCGCTTTCCCCTTCTCACTCAGGCCAAGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGGGGCAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAAGGAAGATCGGTGTGCGCATTTAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCt                                   <  1:334102/251‑1 (MQ=255)                                                                                                                                                                                                    gCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCTGTGTGCTCATTCAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGAt                               <  2:236807/251‑1 (MQ=255)                                                                                                                                                                                                                                   cgcgcgTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCCGGAGATCG‑ATTTGCCAGTCCACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGAttt  >  2:158583/1‑249 (MQ=255)                                                                                                                                                                                                                                        gTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGCATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATc                                                                                                 >  1:160624/1‑150 (MQ=255)                                                                                                                                                                                                                                        gTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATc                                                                                                 <  2:160624/150‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       ATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTCCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTT  >  NC_000913/1610195‑1610669

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