Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	8,972	A→C	67.6%	A245A (GCA→GCC)	talB →	transaldolase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	8,972	0	A	C	67.6%	52.8 / 30.6	40	A245A (GCA→GCC)	talB	transaldolase B
Reads supporting (aligned to +/- strand):  ref base A (6/7);  new base C (17/10);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATA  >  CP009273/8828‑9117                                                                                                                                                 |                                                                                                                                                  ggCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGc                                                                                                                                              >  1:160555/1‑150 (MQ=255)  gCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGcg                                                                                                                                             <  2:88846/150‑2 (MQ=255)   cAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGcgc                                                                                                                                            >  2:127249/1‑148 (MQ=255)    agaagaTCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGct                                                                                                                                           >  2:107126/1‑150 (MQ=255)             gggCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGaaa                                                                                                                                    >  2:160757/1‑148 (MQ=255)                  tGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTgg                                                                                                                             >  2:145340/1‑150 (MQ=255)                          gtcagcgAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGc                                                                                                                     >  1:112853/7‑150 (MQ=255)                                gAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAggc                                                                                                               <  1:127249/150‑2 (MQ=255)                                gAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAggc                                                                                                               <  2:148319/150‑2 (MQ=255)                                  taTCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGgcgc                                                                                                             >  2:14118/2‑150 (MQ=255)                                  aaTCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGgcgc                                                                                                             >  1:163373/1‑150 (MQ=255)                                  aaTCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGgcgc                                                                                                             >  1:82091/1‑150 (MQ=255)                                        cgaGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGgcgccg                                                                                                           <  1:173326/144‑3 (MQ=255)                                        ccAGTACTACAAAGAGCACGGTTATGAAATCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCCGGCGGAGAGCGAAGGGGCTATCg                                                                                                        >  2:125336/1‑149 (MQ=255)                                            tactacAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACg                                                                                                    >  1:122902/1‑149 (MQ=255)                                               tacAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTaag                                                                                                >  2:55266/1‑149 (MQ=255)                                                 caagcagCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCt                                                                                              >  1:76814/6‑150 (MQ=255)                                                 cAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACt                                                                                              <  1:132840/150‑1 (MQ=255)                                                      agCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTggtat                                                                                         <  1:169550/150‑5 (MQ=255)                                                        cACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTac                                                                                       <  2:35894/150‑1 (MQ=255)                                                                  aaaCCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAg                                                                             >  2:72400/1‑150 (MQ=255)                                                                              tGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGt                                                                  <  2:167824/149‑1 (MQ=255)                                                                               gggCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTggtatat                                                                                       >  2:179947/1‑121 (MQ=255)                                                                                ggCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGACCGTAAGCTGTCTTCCTCCGCCGCAGTGAAAGCGCCTcc                                                                <  2:53826/149‑1 (MQ=255)                                                                                  cgcAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTTCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCg                                                               >  2:125500/1‑148 (MQ=255)                                                                                    cAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGgagc                                                           >  1:157288/1‑150 (MQ=255)                                                                                     aaGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGagcg                                                          <  1:114722/150‑1 (MQ=255)                                                                                      aGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGagcg                                                          >  1:54782/1‑149 (MQ=255)                                                                                        cTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTAt                                                       >  2:26170/1‑150 (MQ=255)                                                                                           ccGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATc                                                      <  1:160757/148‑1 (MQ=255)                                                                                                acgttggCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGaa                                                <  2:76814/144‑2 (MQ=255)                                                                                                 cgttggCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATc                                              <  2:176142/146‑1 (MQ=255)                                                                                                   tCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCg                                            <  1:5890/150‑1 (MQ=255)                                                                                                     ggCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCCTCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGCAGGGGCTATCGAACGTAAACTGTCCTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAg                                          >  2:159742/1‑150 (MQ=255)                                                                                                            ttCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCtgt                                   >  1:92404/1‑150 (MQ=255)                                                                                                              cTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCtgtg                                  <  1:139394/149‑1 (MQ=255)                                                                                                                                gACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATcc                >  2:173408/1‑149 (MQ=255)                                                                                                                                gACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATcc                >  1:10300/1‑149 (MQ=255)                                                                                                                                 aCCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCaa              <  2:122902/150‑1 (MQ=255)                                                                                                                                             tCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGaca  <  1:14118/150‑3 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATA  >  CP009273/8828‑9117

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