Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,657 C→T 79.3% I140I (ATC→ATT talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,6570CT79.3% 114.0 / 24.9 53I140I (ATC→ATTtalBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base C (6/5);  new base T (21/21);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

CGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAA  >  CP009273/8511‑8803
                                                                                                                                                  |                                                                                                                                                  
cGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGt                                                                                                                                                 <  2:56898/150‑1 (MQ=255)
       ccaccgAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGcgccgct                                                                                                                                           >  1:26209/6‑148 (MQ=255)
               gTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTg                                                                                                                                  >  1:159872/1‑150 (MQ=255)
                   aTGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGaa                                                                                                                               >  2:158248/1‑149 (MQ=255)
                      cacgTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGaaa                                                                                                                              >  2:169713/3‑147 (MQ=255)
                        cgTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAaagaa                                                                                                                         <  1:163992/150‑1 (MQ=255)
                         gTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGATGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCCCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAaagaag                                                                                                                        <  1:86026/150‑1 (MQ=255)
                             ctCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCAt                                                                                                                    >  2:42370/2‑150 (MQ=255)
                                cTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTaa                                                                                                                 <  2:47724/150‑1 (MQ=255)
                                   tGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAAc                                                                                                                <  1:144326/148‑1 (MQ=255)
                                     acacCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGt                                                                                                             <  1:58100/149‑1 (MQ=255)
                                     acacCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTa                                                                                                            <  2:85418/150‑1 (MQ=255)
                                       cccGAAGCGTCTATCGCGAAAGCCACGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAAc                                                                                                          <  1:61782/149‑1 (MQ=255)
                                          gAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCt                                                                                                        >  2:97631/1‑149 (MQ=255)
                                          gAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCt                                                                                                        >  1:15725/1‑149 (MQ=255)
                                            aGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCCGAc                                                                                                     <  1:167412/150‑1 (MQ=255)
                                             gCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGAcg                                                                                                    >  1:59341/1‑149 (MQ=255)
                                               gTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGct                                                                                                  >  1:140583/1‑150 (MQ=255)
                                                tctatcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGctg                                                                                                 >  1:19050/7‑150 (MQ=255)
                                                      gCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGtt                                                                                            <  2:71850/149‑1 (MQ=255)
                                                      gCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTtc                                                                                           >  2:38708/1‑150 (MQ=255)
                                                       cGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTtct                                                                                          >  2:134773/1‑150 (MQ=255)
                                                       cGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTtc                                                                                           <  2:159872/149‑1 (MQ=255)
                                                                  cgcatcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAg                                                                               >  1:38725/7‑150 (MQ=255)
                                                                              ctgtACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTtgc                                                                   <  1:56999/147‑2 (MQ=255)
                                                                                 tACAACGACCCGGGCATCAGCAACGATCGCATTCTGATCAAACTTGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTgcgcc                                                                >  2:137168/1‑149 (MQ=255)
                                                                                       gATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGGGCGGAAGCg                                                          <  2:9711/150‑1 (MQ=255)
                                                                                            gggcatcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCgt                                                     >  1:126244/8‑150 (MQ=255)
                                                                                             ggcatcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCgtc                                                    >  1:145644/7‑149 (MQ=255)
                                                                                               catcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCtt                                                  >  1:156430/5‑150 (MQ=255)
                                                                                                atcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTc                                                 >  1:194250/4‑150 (MQ=255)
                                                                                                   aGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTg                                              <  1:1544/150‑1 (MQ=255)
                                                                                                        cGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGAtttc                                         >  2:132185/1‑150 (MQ=255)
                                                                                                         gATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCg                                        <  2:59341/150‑1 (MQ=255)
                                                                                                           tCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGcc                                      <  1:29056/150‑1 (MQ=255)
                                                                                                            cgcaTTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCg                                     >  2:115272/4‑150 (MQ=255)
                                                                                                               aTTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGttc                                  <  1:10952/150‑2 (MQ=255)
                                                                                                                     aTCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTgg                             <  1:88225/149‑1 (MQ=255)
                                                                                                                      tCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTGCGCTCAGGCTCGCGCTTGCACCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGcc                           <  2:166432/150‑1 (MQ=255)
                                                                                                                       cAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCg                          >  1:145374/1‑150 (MQ=255)
                                                                                                                       cAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGATGGCCg                          >  2:179167/1‑150 (MQ=255)
                                                                                                                           cTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTAtc                      <  1:132185/150‑2 (MQ=255)
                                                                                                                              gcctccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCtc                   <  1:37766/144‑2 (MQ=255)
                                                                                                                                  ccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACt               >  1:140175/3‑150 (MQ=255)
                                                                                                                                   caCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTg              <  1:42370/149‑1 (MQ=255)
                                                                                                                                     ccTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACt               <  2:87062/147‑1 (MQ=255)
                                                                                                                                     ccTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGt            <  1:169713/150‑1 (MQ=255)
                                                                                                                                      cTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTa           >  2:101798/1‑150 (MQ=255)
                                                                                                                                       tGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTa           >  1:88846/1‑149 (MQ=255)
                                                                                                                                       tGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTa           >  2:114722/1‑149 (MQ=255)
                                                                                                                                        ggCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTAc          >  2:190167/1‑149 (MQ=255)
                                                                                                                                         gCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACa         <  1:4235/149‑1 (MQ=255)
                                                                                                                                               tATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGaa  <  2:15725/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAA  >  CP009273/8511‑8803

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