Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,324,911 T→G 8.0% E68A (GAG→GCG)  obgE ← GTPase involved in cell partioning and DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,324,9110TG 8.0% 848.5 / 11.0 322E68A (GAG→GCG) obgEGTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base T (167/124);  new base G (28/2);  total (196/126)
Fisher's exact test for biased strand distribution p-value = 4.78e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 342 total aligned reads displayed.

CTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCC  >  CP009273/3324846‑3324977
                                                                   |                                                                  
cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGctct                                                                  <  1:8082386/68‑1 (MQ=255)
cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGctct                                                                  <  1:16708194/68‑1 (MQ=255)
cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGctct                                                                  <  1:16002914/68‑1 (MQ=255)
cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGctct                                                                  <  1:25054611/68‑1 (MQ=255)
cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGctct                                                                  <  1:12358270/68‑1 (MQ=255)
  ttAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGc                                                                <  1:5055975/68‑1 (MQ=255)
    aaTCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGGTTGCGCCCTTAT‑‑GCCCCCGCGCTggcc                                                              >  1:3730316/1‑65 (MQ=255)
       cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGGTTGCGCCATTCT‑‑GACCGCGCTCTGGACGGa                                                           >  1:1374130/1‑68 (MQ=255)
       cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTGTGCACGGa                                                           >  1:6804853/1‑68 (MQ=255)
       cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACggt                                                           >  1:26700081/1‑67 (MQ=255)
       cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGa                                                           >  1:28588168/1‑68 (MQ=255)
       cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGa                                                           <  1:27968662/68‑1 (MQ=255)
       cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTAT‑‑GCCCGCGGGCTGCCCggc                                                           >  1:15291886/1‑67 (MQ=255)
       cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCGTGCGCCCGTAT‑‑GTCCGCGCGCTGGCCGGa                                                           >  1:17361382/1‑68 (MQ=255)
              tCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGGCATTCT‑‑GACCGCGCTCTGGACCGAAAGAttt                                                    >  1:11542643/1‑68 (MQ=255)
               cTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGGACGGAAAGAtttt                                                   >  1:24297363/1‑68 (MQ=255)
                   aCCGCGCTTACCGGTACAGTCGCGGGTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTGaa                                               >  1:19227063/1‑68 (MQ=255)
                   aCCGCGCTTACCGGTACAGTCGCGGCTTTGGCCCTTAT‑‑TACCCCGCGGTGCACGGGAAGATTTTTGaa                                               >  1:13192040/1‑68 (MQ=255)
                   aCCGCGCTTACCGGTACAGTCGCGGCTTGGGCCATTCT‑‑GGCCGCGCGCTGCACGGGAAGATTTTTaaa                                               >  1:3422204/1‑68 (MQ=255)
                   aCCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GCCCGCGCTCTGCACGGAAAGATTTTTCaa                                               >  1:19032088/1‑68 (MQ=255)
                   aCCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCaa                                               >  1:20331655/1‑68 (MQ=255)
                   aCCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCaa                                               >  1:1968998/1‑68 (MQ=255)
                   aCCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCaa                                               >  1:2226888/1‑68 (MQ=255)
                   aCCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCGCAGGCCGGAAAGATTTTTCaa                                               >  1:13492898/1‑68 (MQ=255)
                   aCCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTAT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCaa                                               >  1:5990977/1‑68 (MQ=255)
                   aCCGCGCTTACCGGTACAGTCGCGGCGTTCGCCAGTCTTAAACCGCGCTTACAGGAAAGATTTTTaaa                                               >  1:19942157/1‑68 (MQ=255)
                   aCCGCGCTTACCGGTACAGTCGCGGCGTGCGCCCTTCT‑‑GCCCGCGCGCTGCCCGGAAAGATTTTTcta                                               >  1:29171250/1‑66 (MQ=255)
                     cgcgCTTACCGGTACAGTCGCGGGTTGGGCCCGTAT‑‑GACCCCCCGCGGCCCCGTAAGAGTTTTGaaa                                              >  1:6632568/1‑67 (MQ=255)
                     cgcgCTTACCGGTACAGTCGCGGCTTGGGCCATTGT‑‑GACCGCGCTCTGGACGGAAAGATTTTTGaaa                                              >  1:5318779/1‑67 (MQ=255)
                     cgcgCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCaaaa                                             >  1:18088139/1‑68 (MQ=255)
                     cgcgCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCaaa                                              >  1:12007936/1‑67 (MQ=255)
                     cgcgCTTACCGGTACAGTCGCGGCTTGCGCCAATCT‑‑GACCGCGCTCGGCACGGAAAGATTTTTGaaa                                              >  1:26283426/1‑67 (MQ=255)
                        gcTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACg                                           <  1:27869847/67‑1 (MQ=255)
                        gcTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACg                                           <  1:7593529/67‑1 (MQ=255)
                          ttACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGAt                                         <  1:27286254/67‑1 (MQ=255)
                          ttACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGAt                                         <  1:27666059/67‑1 (MQ=255)
                            aCCGGTACTGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAAt                                      <  1:29722407/68‑1 (MQ=255)
                            aCCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAAt                                      <  1:2934816/68‑1 (MQ=255)
                            aCCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAAt                                      <  1:17692565/68‑1 (MQ=255)
                             ccGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATc                                     >  1:4978543/1‑68 (MQ=255)
                              cGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCg                                    <  1:18433961/68‑1 (MQ=255)
                              cGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCg                                    <  1:19478081/68‑1 (MQ=255)
                                gTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGCAAGATTTTTGAAAACGATAATAGAt                                  >  1:940604/1‑68 (MQ=255)
                                 tACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGAt                                  <  1:21540504/67‑1 (MQ=255)
                                  aCAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGGACGGAAAGATTTTTCAAAACGATAATCGATaa                                >  1:14152029/1‑68 (MQ=255)
                                  aCAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGTAAGATTTTTAAAAACGATAATCGATaa                                >  1:18710235/1‑68 (MQ=255)
                                  aCAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTGAAAACGATAATCGATaa                                >  1:26824562/1‑68 (MQ=255)
                                  aCAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                >  1:11386297/1‑68 (MQ=255)
                                  aCAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                >  1:18999706/1‑68 (MQ=255)
                                  aCAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                >  1:18928665/1‑68 (MQ=255)
                                  aCAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                >  1:11753588/1‑68 (MQ=255)
                                  aCAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                >  1:6249468/1‑68 (MQ=255)
                                  aCAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                >  1:5670702/1‑68 (MQ=255)
                                  aCAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTAAAAACGATTATCGATaa                                >  1:20783543/1‑68 (MQ=255)
                                  aCAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTAAAAACGATAATCGATaa                                >  1:28097991/1‑68 (MQ=255)
                                  aCAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTAAAAACGATAATCGATaa                                >  1:4788619/1‑68 (MQ=255)
                                  aCAGTCGCGGCTTGCGCCATTAT‑‑GACCGCGCGCTGCACGGAAAGATTTTTAAAAAAGATAATCGATaa                                >  1:18558191/1‑68 (MQ=255)
                                     gTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAg                               <  1:12609470/66‑1 (MQ=255)
                                       cgcgGCTTGCGCCATTCT‑‑GACGGCGCTCTGCAGGGAAAGATTTTTCAAAACCATTATCGATAAGCgt                            <  1:7422838/67‑1 (MQ=255)
                                       cgcgGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAg                               <  1:3917701/64‑1 (MQ=255)
                                       cgcgGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAg                               <  1:17170169/64‑1 (MQ=255)
                                       cgcgGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGc                              <  1:8870876/65‑1 (MQ=255)
                                       cgcgGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgt                            <  1:11464134/67‑1 (MQ=255)
                                          ggCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGtt                         <  1:23816123/67‑1 (MQ=255)
                                          ggCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTc                        >  1:11018186/1‑68 (MQ=255)
                                            cTTGCGCCATTCT‑‑GACCGCGCTCTGCCCGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAg                      >  1:1562605/1‑68 (MQ=255)
                                            cTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAg                      >  1:13045780/1‑68 (MQ=255)
                                            cTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAg                      >  1:4167829/1‑68 (MQ=255)
                                             ttGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAgg                     <  1:10745799/68‑1 (MQ=255)
                                              tGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGt                    <  1:18277830/68‑1 (MQ=255)
                                              tGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGt                    <  1:28452085/68‑1 (MQ=255)
                                              tGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGt                    <  1:4475288/68‑1 (MQ=255)
                                                 gcCATTCT‑‑GACGGCGCGCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTtct                 <  1:25859905/68‑1 (MQ=255)
                                                 gcCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTtct                 <  1:13925068/68‑1 (MQ=255)
                                                  ccATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTtctc                >  1:20400905/1‑68 (MQ=255)
                                                      tCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCtcgtcg            >  1:25539851/1‑68 (MQ=255)
                                                      tCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCtcgtcg            >  1:20877173/1‑68 (MQ=255)
                                                      tCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCtcgtcg            >  1:6446149/1‑68 (MQ=255)
                                                        t‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGc          <  1:24809335/68‑1 (MQ=255)
                                                        t‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGc          <  1:2002847/68‑1 (MQ=255)
                                                           gACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCt         >  1:24436179/1‑68 (MQ=255)
                                                           gACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCt         >  1:12202008/1‑68 (MQ=255)
                                                            aCCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCt         >  1:24208883/1‑67 (MQ=255)
                                                             ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt        >  1:2601226/1‑67 (MQ=255)
                                                             ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt        >  1:792870/1‑67 (MQ=255)
                                                             ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt        >  1:24021969/1‑67 (MQ=255)
                                                             ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt        >  1:24580732/1‑67 (MQ=255)
                                                             ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt        >  1:909480/1‑67 (MQ=255)
                                                             ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt        >  1:12840410/1‑67 (MQ=255)
                                                             ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt        >  1:8687682/1‑67 (MQ=255)
                                                             ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt        >  1:3182965/1‑67 (MQ=255)
                                                             ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTGCTCGTCGGCtt        >  1:21408283/1‑67 (MQ=255)
                                                               gcgcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtcca     >  1:14293463/1‑68 (MQ=255)
                                                               gcgcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtcca     <  1:20591884/68‑1 (MQ=255)
                                                                cgcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtcca     >  1:29528352/1‑67 (MQ=255)
                                                                 gcTCTGCTCGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGTTTCTCGTCGGCTtcaat    <  1:14596957/67‑1 (MQ=255)
                                                                 gcTCTGCACGGAAAGATTTTTCAAGACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat    <  1:10464830/67‑1 (MQ=255)
                                                                 gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat    <  1:2202931/67‑1 (MQ=255)
                                                                 gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat    <  1:10295086/67‑1 (MQ=255)
                                                                 gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat    <  1:9506609/67‑1 (MQ=255)
                                                                   |                                                                  
CTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCT‑‑GACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCC  >  CP009273/3324846‑3324977

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