Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	8,428	A→G	56.0%	N64S (AAC→AGC)	talB →	transaldolase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	8,428	0	A	G	56.0%	6.5 / 27.7	25	N64S (AAC→AGC)	talB	transaldolase B
Reads supporting (aligned to +/- strand):  ref base A (5/6);  new base G (8/6);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-01

GCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAG  >  CP009273/8301‑8571                                                                                                                                |                                                                                                                                                gcgGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATcctc                                                                                                                           >  2:42749/1‑147 (MQ=255)     cAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATcctcgatg                                                                                                                       >  1:46430/1‑143 (MQ=255)      aaTGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGAAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGAcgc                                                                                                                      >  2:126440/1‑150 (MQ=255)                        gCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAg                                                                                                     >  1:145049/1‑149 (MQ=255)                         cAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGt                                                                                                    >  1:102498/1‑149 (MQ=255)                           ggATGCCACAACCAACCCTTCTCTCATGCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTaaa                                                                                                 <  2:97702/150‑1 (MQ=255)                                      ctaaccCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCtcg                                                                                      <  2:183315/148‑3 (MQ=255)                                                       ctcggcgcggctcagATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTgat                                                                      <  1:104692/144‑3 (MQ=255)                                                             gcggctcagATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGc                                                               >  2:58100/7‑150 (MQ=255)                                                              cggctcagATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGcc                                                              <  1:188405/145‑1 (MQ=255)                                                               agcgcagATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCg                                                             <  2:145049/150‑1 (MQ=255)                                                                             agtACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAg                                               >  2:75239/3‑150 (MQ=255)                                                                              gtACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGt                                              <  1:177441/149‑1 (MQ=255)                                                                               tACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGtt                                             <  2:149455/150‑1 (MQ=255)                                                                               tACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGtt                                             <  2:24705/150‑1 (MQ=255)                                                                                  cGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGac                                          <  1:44743/150‑2 (MQ=255)                                                                                      agctgattgacGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACgcac                                      <  1:48770/147‑3 (MQ=255)                                                                                              gatgatGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGCCCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTGGATGCGCGTCTTTc                               <  2:131046/149‑1 (MQ=255)                                                                                                      tGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGacac                      >  1:34287/1‑150 (MQ=255)                                                                                                            cTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGc                <  1:126440/150‑1 (MQ=255)                                                                                                            cTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGc                >  1:56898/1‑150 (MQ=255)                                                                                                            cTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGc                >  1:106107/1‑150 (MQ=255)                                                                                                                     tagccagAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGc       >  2:86026/5‑150 (MQ=255)                                                                                                                          agagCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAg  >  2:54146/1‑150 (MQ=255)                                                                                                                          agagCAGCGATCGCGCGCAGCAGAGCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAg  >  2:61782/1‑150 (MQ=255)                                                                                                                                |                                                                                                                                                GCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAG  >  CP009273/8301‑8571

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