Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,930,220 C→G 14.8% A118P (GCC→CCC)  zwf ← glucose‑6‑phosphate 1‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,930,2200CG14.8% 159.9 / 18.3 54A118P (GCC→CCC) zwfglucose‑6‑phosphate 1‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (27/19);  new base G (6/2);  total (33/21)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTATTCGCCAACCTGATCATTGATTTCCTGCGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCTTCATCAATGGTTTCTTTCATGAAAG  >  CP009273/1930072‑1930369
                                                                                                                                                    |                                                                                                                                                     
aaGTATTCGCCAACCTGATCATTGATTTCCTGCGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCa                                                                                                                                                      <  1:589211/150‑1 (MQ=255)
 aGTATTCGCCAACCTGATCATTGATTTCCTGCGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCa                                                                                                                                                      >  1:201285/1‑149 (MQ=255)
   tATTCGCCAACCTGATCATTGATTTCCTGCGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAag                                                                                                                                                   >  2:954624/1‑150 (MQ=255)
       cGCCAACCTGATCATTGATTTCCTGCGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAagtagt                                                                                                                                               >  1:59639/1‑150 (MQ=255)
           aaCCTGATCATTGATTTCCTGCGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAAt                                                                                                                                           >  2:542058/1‑150 (MQ=255)
                     ttGATTTCCTGCGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGGAAAGTAGTTAATGGTGATACGa                                                                                                                                 >  1:857543/1‑150 (MQ=255)
                      tGATTTCCTGCGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGAt                                                                                                                                >  1:1010527/1‑150 (MQ=255)
                      tGATTTCCTGCGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGAt                                                                                                                                >  1:816970/1‑150 (MQ=255)
                       gATTTCCTGCGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGGAAAGTAGTTAATGGTGATACGAt                                                                                                                                <  1:579805/149‑1 (MQ=255)
                         ttcccTGCGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGAtttt                                                                                                                             <  2:483290/147‑1 (MQ=255)
                           tCCTGCGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGAttttt                                                                                                                            <  1:830402/149‑1 (MQ=255)
                                cGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGa                                                                                                                        <  1:149602/148‑1 (MQ=255)
                                 gAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATc                                                                                                                      >  1:78279/1‑149 (MQ=255)
                                   ggTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAg                                                                                                                   <  1:799454/150‑1 (MQ=255)
                                    gTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAg                                                                                                                   >  1:528670/1‑149 (MQ=255)
                                      cGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCAt                                                                                                                >  2:158558/1‑150 (MQ=255)
                                               cGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGcc                                                                                                          >  1:1020596/1‑147 (MQ=255)
                                                   cccAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACg                                                                                                    <  1:443708/149‑1 (MQ=255)
                                                    ccAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGGAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGc                                                                                                  >  1:989007/1‑150 (MQ=255)
                                                        cGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGaa                                                                                              >  1:736652/1‑150 (MQ=255)
                                                         ggTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGTCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAAt                                                                                             >  1:802000/1‑150 (MQ=255)
                                                         ggTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAAt                                                                                             >  2:605439/1‑150 (MQ=255)
                                                                 cATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATgcagc                                                                                        <  1:746945/147‑1 (MQ=255)
                                                                 cATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAgt                                                                                     <  2:1010527/150‑1 (MQ=255)
                                                                  aTGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAgtg                                                                                    >  1:824883/1‑150 (MQ=255)
                                                                      cTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCCAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCAt                                                                                >  1:1036411/1‑150 (MQ=255)
                                                                      cTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCAt                                                                                >  2:248359/1‑150 (MQ=255)
                                                                               ccGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATc                                                                         >  1:1008050/1‑148 (MQ=255)
                                                                                       gCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACa                                                                <  2:201285/149‑1 (MQ=255)
                                                                                            cAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATc                                                          >  2:544726/1‑150 (MQ=255)
                                                                                             aGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATcc                                                         <  1:954624/150‑1 (MQ=255)
                                                                                             aGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATcc                                                         <  1:686387/150‑1 (MQ=255)
                                                                                                 ttGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGAc                                                     >  1:582351/1‑150 (MQ=255)
                                                                                                  tGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGAc                                                     <  2:1060677/149‑1 (MQ=255)
                                                                                                    ccTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGGAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTg                                                  >  1:318448/1‑150 (MQ=255)
                                                                                                        gCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACt                                              >  2:571090/1‑150 (MQ=255)
                                                                                                        gCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACt                                              >  1:27784/1‑150 (MQ=255)
                                                                                                          cAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGGAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACt                                              >  2:480320/1‑148 (MQ=255)
                                                                                                             gCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAggg                                         >  2:167588/1‑150 (MQ=255)
                                                                                                                      aaaTTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCAt                                 >  2:94927/1‑149 (MQ=255)
                                                                                                                       aaTTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATa                                >  2:647699/1‑149 (MQ=255)
                                                                                                                         ttGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAAc                             <  1:688744/150‑1 (MQ=255)
                                                                                                                            cgcCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCtt                          >  2:223691/1‑150 (MQ=255)
                                                                                                                                aaaaGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGTATCGCGCTGAGACGGCTGAATGCAGCAGCGCCATTGACATCGAGATTACAAAAATCCAGACGTGGATTCAGGCTGTCCCATAAACCtctat                       >  2:330893/1‑145 (MQ=255)
                                                                                                                                  aaGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCTTCATCaa                    <  2:383590/150‑1 (MQ=255)
                                                                                                                                   aGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCTTCATCAAt                   <  2:59639/150‑1 (MQ=255)
                                                                                                                                          ggcggcATGGGAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCTTCATCAATGGtttct            >  2:869127/1‑150 (MQ=255)
                                                                                                                                          ggcggcATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCTTCATCAATGGtttct            >  1:134953/1‑150 (MQ=255)
                                                                                                                                          ggcggcATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTAAAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCTTCATCAATGGtttc             <  2:623728/149‑1 (MQ=255)
                                                                                                                                           gcggcATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCTTCATCAATGGtttctt           <  2:27784/150‑1 (MQ=255)
                                                                                                                                            cggcATGGGAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCTTCATCAATGGtttcttt          <  2:989007/150‑1 (MQ=255)
                                                                                                                                            cggcATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCTTCATCAATGGtttcttt          <  2:328783/150‑1 (MQ=255)
                                                                                                                                             ggcATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCTTCATCAATGGtttctttc         <  1:248359/150‑1 (MQ=255)
                                                                                                                                              gcATGGGAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCTTCATCAATGGTTTCTTTCa        >  2:642915/1‑150 (MQ=255)
                                                                                                                                                    cAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCTTCATCAATGGTTTCTTTCATGAAAg  <  2:1036411/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGTATTCGCCAACCTGATCATTGATTTCCTGCGAGGTCGCCAGCGACGTCCCCAGCGGTTTCTCCATGACTACGCGTGCCGGTTTAGCATTCAGTTTTGCCTCGCCAAGCCCTTTGCAAATTGCGCCAAAAGTGCTGGGCGGCATGGCAAAGTAGTTAATGGTGATACGATTTTTTTGATCCAGCATCGCGCCGAGACGGCTGAATGCAGCAGTGTCATTGACATCGAGATTACAAAAATCCAGACGTGCACTCAGGGTGTCCCATAAACCTTCATCAATGGTTTCTTTCATGAAAG  >  CP009273/1930072‑1930369

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