Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 3,729,347 G→A E337K (GAG→AAG)  xylR → xylose divergent operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,729,3470GA97.4% 261.6 / ‑3.9 77E337K (GAG→AAG) xylRxylose divergent operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (39/36);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCGAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGATATCGCA  >  CP009273/3729284‑3729411
                                                               |                                                                
tCCAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCGAGaa                                                              >  1:2038947/1‑68 (MQ=255)
tCCAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGaa                                                              <  1:3184005/68‑1 (MQ=255)
tCCAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGaa                                                              <  1:3938859/68‑1 (MQ=255)
tCCAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGaa                                                              <  1:4097646/68‑1 (MQ=255)
tCCAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGaa                                                              <  1:2182083/68‑1 (MQ=255)
cccAATCCTGAGAAGCGTTTTAAAGAAGAGGTGGGAGAAACCATCCATGCCATGATTCATGCCAAGaa                                                              <  1:4316110/67‑1 (MQ=255)
  cAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGc                                                            <  1:4477948/68‑1 (MQ=255)
      cTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCGAGAAGCTGga                                                        >  1:2834079/1‑68 (MQ=255)
      cTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGga                                                        >  1:4194486/1‑68 (MQ=255)
      cTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGga                                                        >  1:2612504/1‑68 (MQ=255)
      cTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGga                                                        >  1:2635389/1‑68 (MQ=255)
      cTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGga                                                        >  1:833056/1‑68 (MQ=255)
        tGAGCAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGgaga                                                      >  1:949000/1‑68 (MQ=255)
        tGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGgaga                                                      >  1:725503/1‑68 (MQ=255)
           gaAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAg                                                   <  1:2772647/68‑1 (MQ=255)
           gaAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAg                                                   >  1:1937732/1‑68 (MQ=255)
           gaAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAg                                                   >  1:2576472/1‑68 (MQ=255)
           gaAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAACAAGCTGGAGAAAg                                                   >  1:2851884/1‑68 (MQ=255)
             aGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAgcg                                                 >  1:490791/1‑68 (MQ=255)
             aGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGAGTCATGCCAAGAAGCTGGAGAAAgcg                                                 >  1:3315954/1‑68 (MQ=255)
               cGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAgcgcg                                               >  1:998886/1‑68 (MQ=255)
                  tttAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAg                                            <  1:2144020/68‑1 (MQ=255)
                    tAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTc                                          <  1:3203431/68‑1 (MQ=255)
                     aaaGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTct                                         >  1:1885763/1‑68 (MQ=255)
                     aaaGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTct                                         <  1:1558092/68‑1 (MQ=255)
                     aaaGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTct                                         >  1:2910083/1‑68 (MQ=255)
                       agaagaGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTctgc                                       >  1:1772223/1‑68 (MQ=255)
                       agaagaGGTGGGTGAAACCATCCATGCCATGATACATGCCAAGAAGCTGGAGAAAGCGCGCAGTctgc                                       >  1:3413798/1‑68 (MQ=255)
                          agagGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGa                                    >  1:51905/1‑68 (MQ=255)
                          agagGTGGGTGAAACCATCCATGCCAAGATTAATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGa                                    >  1:4529650/1‑68 (MQ=255)
                           tagGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGAt                                   <  1:1280382/67‑1 (MQ=255)
                           gagGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGAt                                   <  1:3430158/68‑1 (MQ=255)
                             ggtgggtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAATGCGCGCAGTCTGCTGAttt                                 <  1:2676631/68‑1 (MQ=255)
                             ggtgggtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGactg                                 >  1:4279297/1‑65 (MQ=255)
                             ggtgggtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGAttt                                 <  1:1234471/68‑1 (MQ=255)
                               tgggtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCa                               >  1:3108129/1‑68 (MQ=255)
                               tgggtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCa                               <  1:3528147/68‑1 (MQ=255)
                               tgggtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCa                               >  1:3611141/1‑68 (MQ=255)
                                  gtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAacc                            >  1:2369630/1‑68 (MQ=255)
                                    gAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAaccac                          <  1:1874764/68‑1 (MQ=255)
                                    gAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAaccac                          >  1:926358/1‑68 (MQ=255)
                                     aaaCCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAaccacc                         <  1:1116069/68‑1 (MQ=255)
                                     aaaCCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAaccacc                         <  1:2604415/68‑1 (MQ=255)
                                     aaaCCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAaccacc                         <  1:4149425/68‑1 (MQ=255)
                                     aaaCCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAaccacc                         <  1:2885040/68‑1 (MQ=255)
                                      aaCCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCt                        <  1:4412220/68‑1 (MQ=255)
                                       accaccCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCtt                       <  1:543544/68‑1 (MQ=255)
                                       aCCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCtt                       >  1:3027918/1‑68 (MQ=255)
                                       aCCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCtt                       >  1:1771242/1‑68 (MQ=255)
                                        ccatccatGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTg                      >  1:39170/1‑68 (MQ=255)
                                        ccatccatGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTg                      >  1:4043788/1‑68 (MQ=255)
                                        ccatccatGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTg                      >  1:728785/1‑68 (MQ=255)
                                           tccatGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCg                   <  1:3613532/68‑1 (MQ=255)
                                            ccatgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCCGCTGATTTCAACCACCTTGTCGa                  >  1:2596372/1‑68 (MQ=255)
                                             catgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGAt                 <  1:1759349/68‑1 (MQ=255)
                                             catgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGAt                 <  1:96721/68‑1 (MQ=255)
                                             catgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGAt                 <  1:99848/68‑1 (MQ=255)
                                             catgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGAt                 <  1:43340/68‑1 (MQ=255)
                                             catgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGAt                 <  1:2009529/68‑1 (MQ=255)
                                             catgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGAt                 <  1:323799/68‑1 (MQ=255)
                                             catgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGAt                 <  1:3637795/68‑1 (MQ=255)
                                               tgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCa               >  1:538377/1‑68 (MQ=255)
                                               tgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCa               <  1:710520/68‑1 (MQ=255)
                                               tgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCa               <  1:767858/68‑1 (MQ=255)
                                               tgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCa               <  1:1705724/68‑1 (MQ=255)
                                                 ccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAAt             <  1:4704924/68‑1 (MQ=255)
                                                 ccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAAt             >  1:3075967/1‑68 (MQ=255)
                                                     gATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATgaga         >  1:844497/1‑68 (MQ=255)
                                                      aTTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATgagct        >  1:4210192/1‑66 (MQ=255)
                                                      aTTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGat        >  1:1718996/1‑68 (MQ=255)
                                                       ttCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGata       >  1:3362996/1‑68 (MQ=255)
                                                       ttCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGata       >  1:4250860/1‑68 (MQ=255)
                                                        tCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGatat      >  1:3386725/1‑68 (MQ=255)
                                                         cATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGATATc     <  1:1695924/68‑1 (MQ=255)
                                                         cATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGATATc     <  1:4519616/68‑1 (MQ=255)
                                                            gCCAAGAAGCTGGAGAAGGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGATATCGCa  >  1:3397496/1‑68 (MQ=255)
                                                            gCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGATATCGCa  >  1:3513525/1‑68 (MQ=255)
                                                               |                                                                
TCCAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCGAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGATATCGCA  >  CP009273/3729284‑3729411

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