Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,339,708	G→A	A722T (GCC→ACC)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,339,708	0	G	A	90.4%	149.1 / 7.2	53	A722T (GCC→ACC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base G (0/5);  new base A (9/38);  total (9/44)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

CGAAAAGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTAC  >  W3110S.gb/1339642‑1339771                                                                   |                                                                cGAAAAGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCaa                                                             >  1:2978704/1‑71 (MQ=255) cGAAAAGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCaa                                                             >  1:2151813/1‑71 (MQ=255)         cTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCa                                                     <  1:2765069/71‑1 (MQ=255)         cTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCa                                                     <  1:918253/71‑1 (MQ=255)         cTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCa                                                     <  1:542574/71‑1 (MQ=255)         cTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCa                                                     <  1:2198097/71‑1 (MQ=255)               cTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTa                                               <  1:1878535/71‑1 (MQ=255)               cTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTa                                               <  1:3648111/71‑1 (MQ=255)                     cGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCTCCTTCACCAATATTCGCATCCGTAATGa                                           <  1:1261590/69‑1 (MQ=255)                     cGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGa                                           <  1:1393025/69‑1 (MQ=255)                     cGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGa                                           <  1:1097830/69‑1 (MQ=255)                     cGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGa                                           <  1:3224714/69‑1 (MQ=255)                           gcggcgTGGTAACCATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCa                                                     <  1:2993525/53‑1 (MQ=255)                                        cATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGt                                     >  1:510213/1‑56 (MQ=255)                                        cATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAgg                      >  1:4046056/1‑71 (MQ=255)                                        cATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAgg                      >  1:2623658/1‑71 (MQ=255)                                        cATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAgg                      >  1:3094825/1‑71 (MQ=255)                                        cATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAgg                      >  1:1728092/1‑71 (MQ=255)                                        cATGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAgg                      >  1:887082/1‑71 (MQ=255)                                         aTGAAGTGATGATGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAggg                     <  1:3075341/71‑1 (MQ=255)                                         aTGAAGTGATGATGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAggg                     <  1:696774/71‑1 (MQ=255)                                         aTGAAGTGATGATGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAggg                     <  1:99408/71‑1 (MQ=255)                                         aTGAAGTGATGATGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAggg                     <  1:212832/71‑1 (MQ=255)                                         aTGAAGTGATGATGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAggg                     <  1:1241028/71‑1 (MQ=255)                                         aTGAAGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAggg                     <  1:989803/71‑1 (MQ=255)                                          tGAAGTGATGATGCGCGGCACATTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAgggg                    >  1:812808/1‑71 (MQ=255)                                            aaGTGATGATTCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTg                                    <  1:3716619/53‑1 (MQ=255)                                            aaGTGATGATGCGCGGCACCTTCCCCAATATTCGCATCCGTAATGAAATGGTg                                    <  1:3792438/53‑1 (MQ=255)                                            aaGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTg                                    <  1:4121433/53‑1 (MQ=255)                                            aaGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTg                                    <  1:453172/53‑1 (MQ=255)                                            aaGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTg                                    <  1:2299006/53‑1 (MQ=255)                                            aaGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTg                                    <  1:2395690/53‑1 (MQ=255)                                            aaGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTg                                    <  1:2696366/53‑1 (MQ=255)                                            aaGTGATGATGCGCGGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTg                                    <  1:1166090/53‑1 (MQ=255)                                                  tgatgCGCGGCACCTTCACCAATATTCGCATCCGTAATGa                                           <  1:1635299/40‑1 (MQ=255)                                                       cgcgGCACCTTCACCTATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCa       <  1:3212364/71‑1 (MQ=255)                                                       cgcgGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACTCGGCa       <  1:3852064/71‑1 (MQ=255)                                                       cgcgGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCa       <  1:1673338/71‑1 (MQ=255)                                                       cgcgGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCa       <  1:636171/71‑1 (MQ=255)                                                       cgcgGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCa       <  1:2012837/71‑1 (MQ=255)                                                       cgcgGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCa       <  1:290573/71‑1 (MQ=255)                                                       cgcgGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCa       <  1:2326981/71‑1 (MQ=255)                                                        gcgGCACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCa       <  1:3280248/70‑1 (MQ=255)                                                            cACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTAc  <  1:3837723/71‑1 (MQ=255)                                                            cACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTAc  <  1:2018206/71‑1 (MQ=255)                                                            cACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTAc  <  1:2255163/71‑1 (MQ=255)                                                            cACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTAc  <  1:3527024/71‑1 (MQ=255)                                                            cACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTAc  <  1:3479495/71‑1 (MQ=255)                                                            cACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTAc  <  1:78623/71‑1 (MQ=255)                                                            cACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTAc  <  1:3418222/71‑1 (MQ=255)                                                            cACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTAc  <  1:905988/71‑1 (MQ=255)                                                            cACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTAc  <  1:2312121/71‑1 (MQ=255)                                                            cACCTTCACCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTAc  <  1:2739467/71‑1 (MQ=255)                                                                   |                                                                CGAAAAGACTTTAACTCCTACGGTTCGCGGCGTGGTAACCATGAAGTGATGATGCGCGGCACCTTCGCCAATATTCGCATCCGTAATGAAATGGTGCCTGGCGTTGAAGGGGGGATGACGCGGCATTTAC  >  W3110S.gb/1339642‑1339771

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