Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,324,911 T→G 7.5% 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 7.5% 1305.6 / 11.8 478E68A (GAG→GCG) obgEGTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base T (247/182);  new base G (41/2);  total (294/184)
Fisher's exact test for biased strand distribution p-value = 1.07e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 530 total aligned reads displayed.

ACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCA  >  CP009273/3324845‑3324978
                                                                  |                                                                   
 cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGctct                                                                   <  1:1003954/68‑1 (MQ=255)
 cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGctct                                                                   <  1:40024803/68‑1 (MQ=255)
 cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGctct                                                                   <  1:38546408/68‑1 (MQ=255)
 cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGctct                                                                   <  1:31011204/68‑1 (MQ=255)
 cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGctct                                                                   <  1:29059952/68‑1 (MQ=255)
 cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGctct                                                                   <  1:2864781/68‑1 (MQ=255)
    tAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCa                                                                <  1:42754836/68‑1 (MQ=255)
     aaTCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCAc                                                               >  1:18460971/1‑68 (MQ=255)
      aTCGTCACGTCTTTACCGCGCTTTCCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCTCTCTGCACg                                                              <  1:28248043/68‑1 (MQ=255)
      aTCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACg                                                              <  1:16553664/68‑1 (MQ=255)
        cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGGTTTCGCCAATATGACCCCCCGCTGCACCGa                                                            >  1:25353927/1‑68 (MQ=255)
        cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTGGGCACGGa                                                            >  1:4742993/1‑68 (MQ=255)
        cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGa                                                            >  1:3497830/1‑68 (MQ=255)
        cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCGCTGCACGGa                                                            >  1:11633769/1‑68 (MQ=255)
        cGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCCCGCGCTGGACGGa                                                            >  1:9543015/1‑68 (MQ=255)
          tCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGaaa                                                          >  1:40806731/1‑68 (MQ=255)
           cACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAg                                                         <  1:36993259/68‑1 (MQ=255)
            aCGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGa                                                        <  1:34712765/68‑1 (MQ=255)
               tCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTATGACCGCGCGATGGACGGAAAGAttt                                                     >  1:2142452/1‑68 (MQ=255)
                cTTTACCGCTCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGa                                                            <  1:42409287/60‑1 (MQ=255)
                cTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGAtttt                                                    >  1:24214060/1‑68 (MQ=255)
                    aCCGCGCTTACCGGTACAGTCGCGGCTTGGGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCaa                                                >  1:23071453/1‑68 (MQ=255)
                    aCCGCGCTTACCGGTACAGTCGCGGCTTGGGCCATTCTGACCCCGCTCTGCACGGAAAGATTTTTCaa                                                >  1:10977973/1‑68 (MQ=255)
                    aCCGCGCTTACCGGTACAGTCGCGGCTTGCGCCCTTCTGACCGCGCTCTGCCCGGAAAGATTTTTCaa                                                >  1:31857071/1‑68 (MQ=255)
                      cgcgCTTACCGGTACAGTCGCGGGTTGCGGCATTCTGACCGCGCGCGGCACGGGAAGATTTTTaaaa                                               >  1:6010538/1‑67 (MQ=255)
                      cgcgCTTACCGGTACAGTCGCGGCTTGCGCCCATATGACCCCGCTCTGCACGGAAAGATTTTTaaaa                                               >  1:14063066/1‑67 (MQ=255)
                      cgcgCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCaaaa                                              >  1:33647910/1‑68 (MQ=255)
                      cgcgCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCGTCACGGAAAGATTTTTAaaaa                                              >  1:3061455/1‑68 (MQ=255)
                      cgcgCTTACCGGTACAGTCGCGGCTTGCGCCAATCTGACCGCGGGCTGCACCGAAAGATTTTTCaaa                                               >  1:1506592/1‑67 (MQ=255)
                       gcgcTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAAc                                             <  1:7331582/68‑1 (MQ=255)
                         gcTTACCGGTACAGTCGCGGGGTGCGCCATTCTGACCGCGCGGTGCCCCGAAACATTTTTCAAAAc                                             >  1:35873551/1‑66 (MQ=255)
                         gcTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                                           >  1:9200704/1‑68 (MQ=255)
                           ttACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGAt                                          <  1:19827595/67‑1 (MQ=255)
                           ttACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATa                                         <  1:12248305/68‑1 (MQ=255)
                             aCCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAAt                                       <  1:35659048/68‑1 (MQ=255)
                             aCCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAAt                                       <  1:27575112/68‑1 (MQ=255)
                              ccGGTACAGTCGCGGGTTTGGCCCTTCTGACCCCGCGGGGCCCCGAAAAATTTTTAAAAACGATAAta                                      >  1:14508461/1‑67 (MQ=255)
                              ccGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATc                                      >  1:8236863/1‑68 (MQ=255)
                              ccGGTACAGTCGCGGCTTGCGCCAATCTGACCGCGCGCTGCACCGAAAGATTTTTCAAAAAGATAATc                                      >  1:31702733/1‑68 (MQ=255)
                               cGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATc                                      >  1:29529912/1‑67 (MQ=255)
                               cGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCg                                     <  1:25012950/68‑1 (MQ=255)
                               cGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCg                                     <  1:4373059/68‑1 (MQ=255)
                               cGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCg                                     <  1:24556766/68‑1 (MQ=255)
                                ggTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGa                                    >  1:17969062/1‑68 (MQ=255)
                                   aCAGTCGCGGGTTGCGCCCTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACCATAATAGATaa                                 >  1:22666386/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTTCGCCATTAATACCGCCCACTGCAACTAAACATTTTTAAAAAAGAAAATACATaa                                 >  1:37481195/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                 >  1:8014774/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                 >  1:41718913/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                 >  1:18679506/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACCGAAAGATTTTTCAAAACGATAATCGATaa                                 >  1:30114658/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTCTGACCGCGCTCAGCACGGAAAGATTTTTAAAAACCATAATAGATaa                                 >  1:5321186/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTCTGACCCCCCTCTGCACGGAAAGAATTTTCAAAACCATAAACGCTaa                                 >  1:26002164/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTCAGACCCCGCGCTGCCCCGAAACATTTTTAAAAACGATAATCGATaa                                 >  1:3908873/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTATGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                                 >  1:27989141/1‑68 (MQ=255)
                                   aCAGTCGCGGCTTGCGCCATTATGACCGCGCTCTGCACGGAAAGATTTTTAAAAACGATAATAGATaa                                 >  1:24004661/1‑68 (MQ=255)
                                     aGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGAATTTTCAAAACGATAATCGATAAGc                               >  1:39619730/1‑68 (MQ=255)
                                        cgcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAg                                <  1:40518189/64‑1 (MQ=255)
                                        cgcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAg                                <  1:35272577/64‑1 (MQ=255)
                                        cgcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgt                             <  1:8918429/67‑1 (MQ=255)
                                        cgcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgt                             <  1:40337821/67‑1 (MQ=255)
                                        cgcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgt                             <  1:23446513/67‑1 (MQ=255)
                                        cgcgGCTTGCGCCATTCTGACCCCGCTCTGCACGGAAAGATTTTTCAAAACGATAATAGATAACCgtg                            >  1:43435186/1‑68 (MQ=255)
                                        cgcgGCTTGCGCCATTATGACCCCGCTCTGCACGGAAAGATTTTTAAAAACCATAATAGATAATCgtg                            >  1:1051403/1‑68 (MQ=255)
                                         gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtgt                           <  1:6251765/68‑1 (MQ=255)
                                         gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtgt                           <  1:21130417/68‑1 (MQ=255)
                                         gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtgt                           <  1:20164170/68‑1 (MQ=255)
                                         gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtgt                           <  1:32328164/68‑1 (MQ=255)
                                           ggCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGtt                          <  1:21845532/67‑1 (MQ=255)
                                             cTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAg                       >  1:43102366/1‑68 (MQ=255)
                                             cTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAAACGATAAGCGTGTTCAg                       >  1:37855068/1‑68 (MQ=255)
                                             cTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTACAAACGATAATCGATAATCGTGTTCAg                       >  1:23570684/1‑68 (MQ=255)
                                              ttGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAgg                      <  1:37296629/68‑1 (MQ=255)
                                                       tCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCtcgtcg             >  1:3444068/1‑68 (MQ=255)
                                                       tCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCtcgtcg             >  1:22119406/1‑68 (MQ=255)
                                                       tCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCtcgtcg             >  1:39373599/1‑68 (MQ=255)
                                                       tCTGACCGCCCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCtcgtcg             >  1:28738468/1‑68 (MQ=255)
                                                         tGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGc           <  1:1938951/68‑1 (MQ=255)
                                                         tGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGc           <  1:6622480/68‑1 (MQ=255)
                                                         tGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGc           <  1:13774136/68‑1 (MQ=255)
                                                         tGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGc           <  1:17276696/68‑1 (MQ=255)
                                                         tGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGc           <  1:25670140/68‑1 (MQ=255)
                                                         tGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGc           <  1:32912185/68‑1 (MQ=255)
                                                         tGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGc           <  1:9436318/68‑1 (MQ=255)
                                                          gACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCt          >  1:7668874/1‑68 (MQ=255)
                                                          gACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCt          >  1:5106961/1‑68 (MQ=255)
                                                          gACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCt          >  1:36475238/1‑68 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGcat         >  1:11370558/1‑65 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:15522378/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:32648662/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:41208595/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:7027293/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCtt         >  1:12868898/1‑67 (MQ=255)
                                                            ccGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtc        <  1:16951149/68‑1 (MQ=255)
                                                            ccGCGCGCTGCACGGATAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTGGGCtt         <  1:17669798/67‑1 (MQ=255)
                                                             cgcgCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtcc       <  1:20527861/68‑1 (MQ=255)
                                                             cgcgCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtcc       >  1:26354573/1‑68 (MQ=255)
                                                                gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat     <  1:16076946/67‑1 (MQ=255)
                                                                gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat     <  1:7228591/67‑1 (MQ=255)
                                                                gcTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccat     <  1:27112832/67‑1 (MQ=255)
                                                                  tctGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTtccatcc   <  1:36264762/67‑1 (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: ‑