Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,994 G→T 71.5% A253S (GCG→TCG) ‡ talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,9940GT71.5% 63.5 / 21.2 37A253S (GCG→TCG) ‡talBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base G (4/7);  new base T (16/10);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCC  >  CP009273/8845‑9135
                                                                                                                                                     |                                                                                                                                             
tGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTgg                                                                                                                                               >  2:145340/1‑150 (MQ=255)
        gtcagcgAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGc                                                                                                                                       >  1:112853/7‑150 (MQ=255)
              gAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAggc                                                                                                                                 <  2:148319/150‑2 (MQ=255)
              gAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAggc                                                                                                                                 <  1:127249/150‑2 (MQ=255)
                taTCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGgcgc                                                                                                                               >  2:14118/2‑150 (MQ=255)
                aaTCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGgcgc                                                                                                                               >  1:82091/1‑150 (MQ=255)
                aaTCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGgcgc                                                                                                                               >  1:163373/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)
                                                                                                                                ccGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATccga                                >  1:110610/1‑131 (MQ=255)
                                                                                                                                          tgAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTa     >  1:73313/1‑150 (MQ=255)
                                                                                                                                             aaGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATcc  <  1:26170/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
TGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCC  >  CP009273/8845‑9135

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: