Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,969,337 G→T 49.7% E278* (GAG→TAG)  wecE → TDP‑4‑oxo‑6‑deoxy‑D‑glucose transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,969,3370GT49.7% ‑2.5 / 39.7 34E278* (GAG→TAG) wecETDP‑4‑oxo‑6‑deoxy‑D‑glucose transaminase
Reads supporting (aligned to +/- strand):  ref base G (11/6);  new base T (8/8);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

TGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCGGAAATCATGGCGGTGTTTCATTACATTCCGCTGCA  >  CP009273/3969192‑3969479
                                                                                                                                                 |                                                                                                                                              
tGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGc                                                                                                                                             <  1:322931/149‑1 (MQ=255)
   cTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGc                                                                                                                                          >  1:403079/1‑149 (MQ=255)
       cgcgagttGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCt                                                                                                                                            <  2:280899/137‑1 (MQ=255)
           aTTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTCTCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGAtt                                                                                                                                  <  1:385830/149‑1 (MQ=255)
              tGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCg                                                                                                                                     >  2:109227/1‑143 (MQ=255)
                 tGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGAt                                                                                                                            <  1:182210/149‑1 (MQ=255)
                  gTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGTCGTCGATTc                                                                                                                                 >  2:430022/1‑143 (MQ=255)
                   tCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGCCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGCTTcc                                                                                                                                >  2:399713/1‑143 (MQ=255)
                     cGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAAGGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGCAAGCCGGGCGTATCGAGCTGCCATCGATTCCCg                                                                                                                              >  2:247737/1‑143 (MQ=255)
                     cGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGCTGGCt                                                                                                                        >  1:318970/1‑149 (MQ=255)
                     cGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCt                                                                                                                        >  1:318934/1‑149 (MQ=255)
                            cAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCa                                                                                                                 >  1:501425/1‑149 (MQ=255)
                                           tgGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGtt                                                                                                  >  1:396470/1‑149 (MQ=255)
                                              gcgcAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTa                                                                                               <  1:477880/149‑1 (MQ=255)
                                              gcgcAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTAATACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTa                                                                                               >  1:521603/1‑149 (MQ=255)
                                                  aactcaaaGCAGCGCATCGAATCAACCAGCAACGCCTGGCGCTGTGGCAAAACCACTACGATGCGTTACCCCCTCCGCCGAAACCCGGGCCTACCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGtcc                                                                                                 <  2:433528/137‑3 (MQ=255)
                                                   aCTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTa                                                                                          <  1:78278/149‑1 (MQ=255)
                                                           ctgccgcatatgataaCTAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATTCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATAAGTTCTACATTcaa                                                                                        >  2:4372/15‑143 (MQ=255)
                                                           cagcGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTCCTCCGGCGCATGGTCGCCTCTGGCGATTGCTGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCgg                                                                                  <  1:232147/149‑1 (MQ=255)
                                                              cGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTg                                                                                     >  2:416242/1‑143 (MQ=255)
                                                                     tATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATtgatg                                                                        >  1:34383/1‑149 (MQ=255)
                                                                       tCAACCAGCAACGCCGGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGACAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGGACGCGCATATGTGCTACAGTAAAATGCGGGCTAtt                                                                            >  2:260108/1‑143 (MQ=255)
                                                                        cAACCAGCAACGTCTGGCGCGGTGGCAAAACTACTACGATGCGTGAGCGTCTCTGGCGAAAGGCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATtg                                                                           <  2:389937/143‑1 (MQ=255)
                                                                            cAGCAACGTCTGGCGCTGGGGCGAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATtgatga                                                                       <  2:252065/143‑1 (MQ=255)
                                                                                          gcTGGGGCAAAACCACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGa                                                         <  2:79520/143‑1 (MQ=255)
                                                                                                gCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGa                                             >  1:166041/1‑149 (MQ=255)
                                                                                                gCAAAACTACTACGATGCGTTAGCGCCTCAGGCGAAAGCCGGGCGTATCGAGCAGCCGTCGATTCCCCATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGATCGGAGCGCGTAGATTAACTTTCTGa                                             >  1:322036/1‑149 (MQ=255)
                                                                                                             gATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGCTTAACTTTGTGAAAGAAGCGGAAAt                                >  1:396079/1‑149 (MQ=255)
                                                                                                               tGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCGGAAATCa                              >  1:79392/1‑149 (MQ=255)
                                                                                                                    tAGCGCCCCTGGCGAAAGCCGGGTGTATCTAGCTGCCGTCGTTTCATGGATGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCGGAAATCATGGCg                         <  1:464780/149‑1 (MQ=255)
                                                                                                                          ctctGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGGGCGCGTTGATTAACTTTTTGAAAGATGCGGagttcgtggag                         >  2:230834/1‑133 (MQ=255)
                                                                                                                               gCGAAAGCCGGGCGTATCCAGCTGCCGCCGATTCCCGAGGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCGGAAATCATGGCGGTGTTTCATTa              <  1:2127/149‑1 (MQ=255)
                                                                                                                                    cgCCGGGCGTCTCGAGCAGCCGTCGAGTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAACGAAGCGGAAATCATGGCGGTGTTTCATTACATTc         <  1:278991/148‑1 (MQ=255)
                                                                                                                                           cGTAGCGAGCTGCCGTCGATTCCCGATGGCGGCCTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCGGAAATCATGGCGGTGTTTCATTACATTCCGCTGCa  <  1:70978/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCGGAAATCATGGCGGTGTTTCATTACATTCCGCTGCA  >  CP009273/3969192‑3969479

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: