Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,324,911 T→G 8.9% 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.9% 1036.2 / 13.3 378E68A (GAG→GCG) obgEGTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base T (201/136);  new base G (37/2);  total (240/138)
Fisher's exact test for biased strand distribution p-value = 2.42e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 402 total aligned reads displayed.

ACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCA  >  CP009273/3324845‑3324978
                                                                  |                                                                   
 cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGctct                                                                   <  1:27355983/68‑1 (MQ=255)
 cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGctct                                                                   <  1:24784697/68‑1 (MQ=255)
    tAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCa                                                                <  1:14010989/68‑1 (MQ=255)
     aaTCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCGTGCGCCCTTCTGACCGCGCTCTGccc                                                               >  1:3458313/1‑66 (MQ=255)
        cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGGTTGGGCCATTATGACCGCGCGGTGCACGGa                                                            >  1:26085021/1‑68 (MQ=255)
        cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTTCGCCATTCTGACCGCGCGCTGGACGGa                                                            >  1:21717196/1‑68 (MQ=255)
        cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCAATATGACCCCGCTCTGGACggt                                                            >  1:5265330/1‑67 (MQ=255)
           cACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAg                                                         <  1:1098865/68‑1 (MQ=255)
           cACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAg                                                         <  1:11532683/68‑1 (MQ=255)
              gTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGAtt                                                      >  1:5862164/1‑68 (MQ=255)
                cTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGAtttt                                                    >  1:9188045/1‑68 (MQ=255)
                cTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGAtttt                                                    >  1:31496708/1‑68 (MQ=255)
                    aCCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTGTGACCGCGGTCTGCACGGAAAGATTTTTaaa                                                >  1:27149664/1‑68 (MQ=255)
                    aCCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACAGCGCTGTGCACGGAAAGATTTTTCaa                                                >  1:31040966/1‑68 (MQ=255)
                    aCCGCGCTTACCGGTACAGTCGCGGCGTGCGCCATTCTGACCGCGCTCGGCACGGAAAGATTTTTGaa                                                >  1:823481/1‑68 (MQ=255)
                     ccGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTATGACCGCGCTCTGCACGGAAAGATTTTTCaaa                                               >  1:32194270/1‑68 (MQ=255)
                     ccGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTATGACCGCGCGCTGCACGGAAAGATTTTTCaaa                                               >  1:30133212/1‑68 (MQ=255)
                      cgcgCTTACCGGTACAGTCGCGGCTTGGGCCATTCAGAACCCGCACTGCACGGAAAGAATTTTCaaa                                               >  1:21814574/1‑67 (MQ=255)
                      cgcgCTTACCGGTACAGTCGCGGCTTGCGGCATTCTGACCGCGGTCTGCCCGGTAAGAATTTTGaaa                                               >  1:15960919/1‑67 (MQ=255)
                      cgcgCTTACCGGTACAGTCGCGGCTTGCGCCCTTCTGACCGCGCGCTGGACGGAAAGATTTTTCaaa                                               >  1:25117055/1‑67 (MQ=255)
                      cgcgCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGGTCTGCACGGAAAGATTTTTCaaa                                               >  1:22680968/1‑67 (MQ=255)
                      cgcgCTTACCGGTACAGTCGCGGCTTGCGCCATTATGACCGCGCTCTGCACGGAAAGATTTTTCaaa                                               >  1:7347491/1‑67 (MQ=255)
                         gcTTACCGGTACAGTCGCGGGTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                                           >  1:19687340/1‑68 (MQ=255)
                         gcTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                                           <  1:15278999/68‑1 (MQ=255)
                             aCCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAAt                                       <  1:3246503/68‑1 (MQ=255)
                             aCCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAAt                                       <  1:20978757/68‑1 (MQ=255)
                             aCCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAAt                                       <  1:20532725/68‑1 (MQ=255)
                             aCCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAAt                                       <  1:19241353/68‑1 (MQ=255)
                             aCCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAAt                                       <  1:20355698/68‑1 (MQ=255)
                             aCCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAAt                                       <  1:7001889/68‑1 (MQ=255)
                               cGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCg                                     <  1:4226925/68‑1 (MQ=255)
                               cGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCg                                     >  1:28348832/1‑68 (MQ=255)
                               cGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTATGCACGGAAAGATTTTTCAAAACGATAATCg                                     <  1:13240984/68‑1 (MQ=255)
                               cGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCGCTGCACGGAAAGATTTTTCAAAACGATAATCg                                     >  1:3272131/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTTCGCCAATCTGCCCCCGCGATGCCCGGAAAGAATTTTCAAAACCCTAATCGATaa                                 >  1:5040644/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCCTTATGACCGCCCTCTGCACCGAAAGATTTTTGAAAACGATAATCGCTaa                                 >  1:2357373/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGGACGGAAAGATTTTTAAAAACGATAATCGATaa                                 >  1:12345054/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                 >  1:14406364/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                 >  1:17547765/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTAAAAACGATAATCGATaa                                 >  1:6605161/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTCTGACCGCGCGCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                 >  1:1886953/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTATGACCGCGCTCTGCACGGCAAGATTTTTCAAAACGATAATCGATaa                                 >  1:1501428/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTATGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                 >  1:27672019/1‑68 (MQ=255)
                                     aGTCGCGGGTTGCGCCATTCTGACCGCGCGCTGCACGGAAAGATTTTTAAAAACGATAATCGATAAGc                               >  1:3828406/1‑68 (MQ=255)
                                     aGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATAGATAAGc                               >  1:23996561/1‑68 (MQ=255)
                                     aGTCGCGGCTTGCGCCATTCTGACCCCGCGCTGCACGGAAAGACTTTTAAAAAAGATAATCGATAAGc                               >  1:22102296/1‑68 (MQ=255)
                                     aGTCGCGGCTTGCGCCAATCAGACCGCGCACTGGACGGAAAGATTTTTCAAAACGAAAATAGATAAGc                               >  1:28809828/1‑68 (MQ=255)
                                        cgcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAg                                <  1:29568226/64‑1 (MQ=255)
                                        cgcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAg                                <  1:33783309/64‑1 (MQ=255)
                                        cgcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAg                                <  1:6396349/64‑1 (MQ=255)
                                        cgcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAg                                <  1:22320186/64‑1 (MQ=255)
                                        cgcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAg                                <  1:6768002/64‑1 (MQ=255)
                                        cgcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgt                             <  1:20150215/67‑1 (MQ=255)
                                        cgcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCAATAAg                                <  1:1383815/64‑1 (MQ=255)
                                        cgcgGCTTGCGCCATTATGACCGCGCTCTGCACGGAAAGATTTTTAAAAAACATAATAGATAAGCgtg                            >  1:30391778/1‑68 (MQ=255)
                                         gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtgt                           <  1:30925457/68‑1 (MQ=255)
                                         gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtgt                           <  1:25942785/68‑1 (MQ=255)
                                         gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtgt                           <  1:23694054/68‑1 (MQ=255)
                                         gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtgt                           <  1:31619009/68‑1 (MQ=255)
                                         gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtgt                           <  1:11370485/68‑1 (MQ=255)
                                         gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtgt                           <  1:33113252/68‑1 (MQ=255)
                                           ggCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGtt                          <  1:11744076/67‑1 (MQ=255)
                                             cTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAg                       >  1:28119680/1‑68 (MQ=255)
                                             cTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAg                       >  1:16748452/1‑68 (MQ=255)
                                             cTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAg                       >  1:24189776/1‑68 (MQ=255)
                                              ttGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAgg                      <  1:10505580/68‑1 (MQ=255)
                                              ttGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAgg                      <  1:10209640/68‑1 (MQ=255)
                                                   ccATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTtctc                 >  1:14718709/1‑68 (MQ=255)
                                                   ccATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTtctc                 >  1:3916700/1‑68 (MQ=255)
                                                    cATTCTGACCGCGCGCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCtcg                <  1:14177695/68‑1 (MQ=255)
                                                       tCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCtcgtcg             >  1:16233867/1‑68 (MQ=255)
                                                       tCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCtcgtcg             >  1:16935467/1‑68 (MQ=255)
                                                       tCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCtcgtcg             >  1:34252444/1‑68 (MQ=255)
                                                        cTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCtcgtcg             >  1:4838929/1‑67 (MQ=255)
                                                         tGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGc           <  1:7169012/68‑1 (MQ=255)
                                                          gACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCt          <  1:33494583/68‑1 (MQ=255)
                                                          gACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCt          >  1:31926555/1‑68 (MQ=255)
                                                           aCCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCt          >  1:3627841/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:289715/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:33309021/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:8823725/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:7992457/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:7418126/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:25653526/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:24424270/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:566523/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:21375187/1‑67 (MQ=255)
                                                             cgcgCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtcc       <  1:14455665/68‑1 (MQ=255)
                                                             cgcgCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtcc       >  1:13541984/1‑68 (MQ=255)
                                                               cgcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtcca      >  1:19078092/1‑67 (MQ=255)
                                                                gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat     <  1:17292573/67‑1 (MQ=255)
                                                                gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat     >  1:17772614/1‑67 (MQ=255)
                                                                gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat     <  1:15733672/67‑1 (MQ=255)
                                                                gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat     <  1:18086724/67‑1 (MQ=255)
                                                                gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat     <  1:26721031/67‑1 (MQ=255)
                                                                gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat     <  1:18321196/67‑1 (MQ=255)
                                                                gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat     <  1:7751204/67‑1 (MQ=255)
                                                                gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat     <  1:19295213/67‑1 (MQ=255)
                                                                 ctctGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccatcc   >  1:30611510/1‑68 (MQ=255)
                                                                 ctctGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccatcc   >  1:23168278/1‑68 (MQ=255)
                                                                  |                                                                   
ACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCA  >  CP009273/3324845‑3324978

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: ‑