Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 9,065 C→T 72.0% I276I (ATC→ATT talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092739,0650CT72.0% 54.1 / 21.0 35I276I (ATC→ATTtalBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base C (3/7);  new base T (11/14);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTC  >  CP009273/8918‑9197
                                                                                                                                                   |                                                                                                                                    
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)
                                                                            tccgaaagCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGtt                                                         <  2:112853/143‑1 (MQ=255)
                                                                             cGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTg                                                       <  1:147552/150‑1 (MQ=255)
                                                                                aaagCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGcga                                                    >  1:99271/3‑148 (MQ=255)
                                                                                  agCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGAtc                                                  <  1:134147/150‑2 (MQ=255)
                                                                                     gAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGAc                                               >  2:123044/1‑150 (MQ=255)
                                                                                          ggCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGa                                          <  1:72400/150‑1 (MQ=255)
                                                                                           gCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGaa                                         <  1:125336/150‑1 (MQ=255)
                                                                                               tCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGaaaaa                                      <  1:145340/149‑1 (MQ=255)
                                                                                                             gtatatacCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATc                       >  1:10501/7‑150 (MQ=255)
                                                                                                              tatatacCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCg                      <  1:138557/145‑1 (MQ=255)
                                                                                                                tatacCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCgac                    <  1:190215/147‑3 (MQ=255)
                                                                                                                  tacCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACg                   >  2:152619/2‑149 (MQ=255)
                                                                                                                     cGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATct               <  1:159742/150‑1 (MQ=255)
                                                                                                                      ggtgAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActg              >  1:71862/4‑150 (MQ=255)
                                                                                                                      ggtgAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActg              <  2:183261/147‑1 (MQ=255)
                                                                                                                        tgAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActgct            <  2:92404/149‑1 (MQ=255)
                                                                                                                        cGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATctgct            >  2:160600/1‑150 (MQ=255)
                                                                                                                         gAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActgctg           >  1:123210/1‑150 (MQ=255)
                                                                                                                          aaGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAACTGCTGt          <  2:82091/150‑1 (MQ=255)
                                                                                                                                  gcgcGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAACTGCTGTaagagaat  <  1:63609/150‑7 (MQ=255)
                                                                                                                                      gTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActgct            >  2:64725/1‑136 (MQ=255)
                                                                                                                                       tCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActgctg           <  1:64725/136‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                    
CCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTC  >  CP009273/8918‑9197

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