Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,241,035	G→A	6.3%	C312C (TGC→TGT)	uxaA ←	altronate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,241,035	0	G	A	6.3%	191.0 / 2.3	63	C312C (TGC→TGT)	uxaA	altronate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base A (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CATGGGCAATAAAGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGT  >  W3110S.gb/3240965‑3241105                                                                       |                                                                       cATGGGCAATAAAGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCgg                                                                        >  1:74632/1‑71 (MQ=255) cATGGGCAATAAAGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCgg                                                                        >  1:1702833/1‑71 (MQ=255)     ggCAATAAAGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCaa                                                                     <  1:2448467/70‑1 (MQ=255)     ggCAATAAAGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCaa                                                                     <  1:1619029/70‑1 (MQ=255)       cAATAAAGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATgg                                                                  >  1:1261674/1‑71 (MQ=255)       cAATAAAGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATg                                                                   >  1:2027960/1‑70 (MQ=255)       cAATAAAGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATg                                                                   >  1:547879/1‑70 (MQ=255)         aTAAAGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTc                                                                >  1:2293851/1‑71 (MQ=255)         aTAAAGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTc                                                                <  1:1872804/71‑1 (MQ=255)         aTAAAGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTc                                                                <  1:3829765/71‑1 (MQ=255)          tAAAGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCaa                                                                     <  1:3850243/65‑1 (MQ=255)           aaaGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTcc                                                               >  1:347533/1‑70 (MQ=255)           aaaGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCa                                                              >  1:3457713/1‑71 (MQ=255)           aaaGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCa                                                              >  1:2031564/1‑71 (MQ=255)           aaaGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCa                                                              >  1:2963929/1‑71 (MQ=255)            aaGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGt                                                                 >  1:3309019/1‑67 (MQ=255)            aaGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCAt                                                             >  1:1725623/1‑71 (MQ=255)            aaGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCAt                                                             >  1:3477212/1‑71 (MQ=255)               tACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGc                                                                       <  1:1334475/58‑1 (MQ=255)               tACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGc                                                                       <  1:2169385/58‑1 (MQ=255)               tACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGc                                                                       <  1:2173811/58‑1 (MQ=255)               tACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGc                                                                       <  1:2182152/58‑1 (MQ=255)               tACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGc                                                                       <  1:822934/58‑1 (MQ=255)                 cTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCaa                                                                     >  1:3711844/1‑58 (MQ=255)                 cTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGc                                                         <  1:2070101/70‑1 (MQ=255)                   gTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAAc                                                      >  1:2913578/1‑71 (MQ=255)                    tttGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATg                                                                   <  1:3131301/57‑1 (MQ=255)                      tGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTg                                                    <  1:1324784/70‑1 (MQ=255)                      tGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTg                                                    <  1:1127092/70‑1 (MQ=255)                      tGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTg                                                    <  1:2434720/70‑1 (MQ=255)                      tGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTg                                                    <  1:87802/70‑1 (MQ=255)                      tGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGc                                                   >  1:3620929/1‑71 (MQ=255)                          gTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCaa                                                                     <  1:570548/49‑1 (MQ=255)                           tCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGc                                                                       <  1:3568905/46‑1 (MQ=255)                             gTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTgcgc                                            >  1:3688496/1‑71 (MQ=255)                                gACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCa                                          >  1:880245/1‑70 (MQ=255)                                   cATGGTGACCAGTTTTTCAAATGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAAc                                       <  1:2631810/70‑1 (MQ=255)                                   cATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGctct                                                <  1:1126534/61‑1 (MQ=255)                                   cATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAAc                                       <  1:3222848/70‑1 (MQ=255)                                   cATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAAc                                       <  1:596008/70‑1 (MQ=255)                                     tGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGc                                                                       <  1:1020225/36‑1 (MQ=255)                                     tGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGc                                                                       <  1:1128018/36‑1 (MQ=255)                                     tGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACAt                                     >  1:1571380/1‑70 (MQ=255)                                       gTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATg                                                                   >  1:1015137/1‑38 (MQ=255)                                         gACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCaaa                                        >  1:1058054/1‑63 (MQ=255)                                               tttttCAAACGTTGCTTCGTCGCGGCAATGGTCCAt                                                             <  1:2124567/36‑1 (MQ=255)                                                 tttCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGctc                                                 >  1:1739174/1‑46 (MQ=255)                                                 tttCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCg                        >  1:2983487/1‑71 (MQ=255)                                                 tttCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCg                        >  1:107329/1‑71 (MQ=255)                                                 tttCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCg                        >  1:168892/1‑71 (MQ=255)                                                 tttCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCg                        >  1:2028270/1‑71 (MQ=255)                                                 tttCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCg                        >  1:2595723/1‑71 (MQ=255)                                                 tttCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCg                        >  1:3296873/1‑71 (MQ=255)                                                  ttCAAACGTTGCTTCGGCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCgg                               <  1:3627257/63‑1 (MQ=255)                                                       aCGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAg                   >  1:1639655/1‑70 (MQ=255)                                                       aCGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAg                   >  1:3109686/1‑70 (MQ=255)                                                       aCGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGt                  <  1:1755250/71‑1 (MQ=255)                                                        cGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTa                 <  1:706469/71‑1 (MQ=255)                                                            gCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACgg              >  1:601086/1‑70 (MQ=255)                                                             cTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACgg              >  1:3204835/1‑69 (MQ=255)                                                             cTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGgt             <  1:271535/70‑1 (MQ=255)                                                                cgtcgCGACAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGt                      >  1:3461815/1‑58 (MQ=255)                                                                cgtcgCGACAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGt                      >  1:2698641/1‑58 (MQ=255)                                                                 gtcgCGACAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTAc        <  1:2520031/71‑1 (MQ=255)                                                                 gtcgCGACAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTAc        <  1:3383069/71‑1 (MQ=255)                                                                       gCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGt  <  1:1666726/71‑1 (MQ=255)                                                                       |                                                                       CATGGGCAATAAAGTACTGTTTGAAGTCGTTGACCATGGTGACCAGTTTTTCAAACGTTGCTTCGTCGCGGCAATGGTCCATCAGCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGT  >  W3110S.gb/3240965‑3241105

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