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,3470GA95.5% 280.4 / 1.9 88E337K (GAG→AAG) xylRxylose divergent operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (39/45);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TCCAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCGAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGATATCGCA  >  CP009273/3729284‑3729411
                                                               |                                                                
tCCAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGaa                                                              <  1:889102/68‑1 (MQ=255)
tCCAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGaa                                                              <  1:4346001/68‑1 (MQ=255)
tCCAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGaa                                                              >  1:3560970/1‑68 (MQ=255)
tCCAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGaa                                                              <  1:720093/68‑1 (MQ=255)
tCCAATCTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGaa                                                              <  1:5207182/68‑1 (MQ=255)
  ctttCTTGAGAAGCGTTTTAAAGAAGAGGTGGTTGAAACCATCCATGCCATGATTCATGCCATGAAGc                                                            <  1:815306/65‑1 (MQ=255)
      cTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGga                                                        >  1:339689/1‑68 (MQ=255)
      cTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGga                                                        <  1:3452328/68‑1 (MQ=255)
      cTTGAGAAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGga                                                        <  1:5677947/68‑1 (MQ=255)
          agaAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGaaa                                                    >  1:3486229/1‑68 (MQ=255)
           gaAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAg                                                   >  1:4509918/1‑68 (MQ=255)
           gaAGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAg                                                   <  1:2993298/68‑1 (MQ=255)
            aaGAGTTTTTAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAgc                                                  <  1:299180/68‑1 (MQ=255)
             aGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAgcg                                                 >  1:454175/1‑68 (MQ=255)
             aGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAgcg                                                 >  1:1666766/1‑68 (MQ=255)
             aGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAgcg                                                 >  1:3966653/1‑68 (MQ=255)
             aGCGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAgcg                                                 >  1:556043/1‑68 (MQ=255)
             aGCGTTTTAAAGAAGAGGTGGGGGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAgcg                                                 <  1:167842/68‑1 (MQ=255)
               cGTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCGAGAAGCTGGAGAAAgcgcg                                               >  1:3781724/1‑68 (MQ=255)
                gTTTTAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGTAAgcgcgc                                              >  1:830318/1‑68 (MQ=255)
                 ttttAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCa                                             <  1:2932592/68‑1 (MQ=255)
                  tttAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAg                                            <  1:2738360/68‑1 (MQ=255)
                   ttAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCGAGAAGCTGGAGAAAGCGCGCAGt                                           <  1:4842769/68‑1 (MQ=255)
                   ttAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGt                                           <  1:140232/68‑1 (MQ=255)
                   ttAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGt                                           <  1:1317788/68‑1 (MQ=255)
                   ttAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGt                                           <  1:1089425/68‑1 (MQ=255)
                   ttAAAGAAGAGGTGGGGGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGt                                           <  1:4100914/68‑1 (MQ=255)
                    tAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTc                                          <  1:249809/68‑1 (MQ=255)
                    tAAAGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTc                                          <  1:5382230/68‑1 (MQ=255)
                     aaaGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTct                                         <  1:4249385/68‑1 (MQ=255)
                     aaaGAAGAGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTct                                         <  1:3140538/68‑1 (MQ=255)
                       agaagaGGTGGGTGATACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTctgc                                       >  1:4906816/1‑68 (MQ=255)
                       agaagaGGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTctgc                                       <  1:289160/68‑1 (MQ=255)
                          agagGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGATGCTGGAGAAAGCGCGCAGTCTGCTGa                                    >  1:1954433/1‑68 (MQ=255)
                          agagGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGa                                    <  1:2429937/68‑1 (MQ=255)
                          agagGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGa                                    <  1:3846521/68‑1 (MQ=255)
                          agagGTGGGGGAAACCATCCATGCCATGATTCATGCCAAGAAGCTTGAGAAAGCGCGCAGTCTGCTGa                                    <  1:4700298/68‑1 (MQ=255)
                          agagGTGGGGGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGa                                    <  1:4700300/68‑1 (MQ=255)
                            agGTGGGTGAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGAtt                                  >  1:2321635/1‑68 (MQ=255)
                             ggtgggtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGCttt                                 <  1:1289279/68‑1 (MQ=255)
                             ggtgggtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGAttt                                 >  1:4565245/1‑68 (MQ=255)
                             ggtgggagAAACCATCCATGCCATGAATCATGCCGAGAAGCTGGAGAAACCGCGCAGTCTGCTGAttt                                 <  1:1463579/68‑1 (MQ=255)
                               tgggtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCa                               >  1:1553746/1‑68 (MQ=255)
                                  gtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAacc                            >  1:3522872/1‑68 (MQ=255)
                                  gtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAacc                            >  1:5527166/1‑68 (MQ=255)
                                  gtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAacc                            <  1:520060/68‑1 (MQ=255)
                                  gtgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAacc                            >  1:5601627/1‑68 (MQ=255)
                                  gggAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAacc                            <  1:5851824/66‑1 (MQ=255)
                                   tgAAACCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTCAaccac                          <  1:53149/68‑1 (MQ=255)
                                      aaCCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCt                        <  1:1006709/68‑1 (MQ=255)
                                       aCCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCTGTCTGCTGATTTCAACCACCtt                       >  1:2794263/1‑68 (MQ=255)
                                       aCCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCtt                       <  1:1555104/68‑1 (MQ=255)
                                       aCCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCtt                       >  1:1422285/1‑68 (MQ=255)
                                       aCCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCtt                       >  1:4924305/1‑68 (MQ=255)
                                       aCCATCCATGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCtt                       <  1:4881092/68‑1 (MQ=255)
                                           tccatGCCATGATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCg                   >  1:1330140/1‑68 (MQ=255)
                                            ccatgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGa                  >  1:4123474/1‑68 (MQ=255)
                                            ccaagccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGa                  <  1:1580310/68‑1 (MQ=255)
                                             catgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGAt                 <  1:2550402/68‑1 (MQ=255)
                                             catgcaatgATTCATGCCAAGAAGCTTGAGAAAGCGCGCAGTCTGCTGTTATCAACCACCTTGTCGAt                 <  1:3204923/68‑1 (MQ=255)
                                               tgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCa               <  1:45507/68‑1 (MQ=255)
                                               tgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCa               <  1:3994989/68‑1 (MQ=255)
                                               tgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCa               <  1:4901152/68‑1 (MQ=255)
                                               tgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCa               <  1:2761063/68‑1 (MQ=255)
                                               tgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCa               <  1:2129245/68‑1 (MQ=255)
                                               tgccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCa               <  1:166570/68‑1 (MQ=255)
                                               tgccatgAATCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATATCAACCACCTTGTCGATCa               <  1:585374/68‑1 (MQ=255)
                                                 ccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAAt             >  1:5180939/1‑68 (MQ=255)
                                                 ccatgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAAt             >  1:1091798/1‑68 (MQ=255)
                                                  catgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATg            >  1:4726737/1‑68 (MQ=255)
                                                  catgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATg            >  1:4735681/1‑68 (MQ=255)
                                                  catgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATg            >  1:4798336/1‑68 (MQ=255)
                                                  catgATTCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATg            >  1:5381058/1‑68 (MQ=255)
                                                     gaatCATGCCAAGAAGCTGGAGAAAGCGCGCAGGCTGCTGATTTCAACCACCTTGTCGATCAATgaga         <  1:3204010/65‑1 (MQ=255)
                                                       ttCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGata       >  1:980260/1‑68 (MQ=255)
                                                       ttCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGata       >  1:3749317/1‑68 (MQ=255)
                                                       ttCATGCCAAGAAGCTGGAGAAAGCGCGCAGGCTGCTGATTTCAACCACCTTGTCGATCAATGAGata       <  1:1105659/68‑1 (MQ=255)
                                                        tCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGatat      >  1:5121910/1‑68 (MQ=255)
                                                        tCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGatat      >  1:3994049/1‑68 (MQ=255)
                                                        tCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGatat      >  1:3070774/1‑68 (MQ=255)
                                                        tCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGatat      >  1:2485550/1‑68 (MQ=255)
                                                        tCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGatat      >  1:186291/1‑68 (MQ=255)
                                                        tCATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGatat      >  1:1160491/1‑68 (MQ=255)
                                                         cATGCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGATATc     <  1:6079451/68‑1 (MQ=255)
                                                            gCCGAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGATATCGCa  >  1:3147150/1‑68 (MQ=255)
                                                            gCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTGATTTCAACCACCTTGTCGATCAATGAGATATCGCa  >  1:4050074/1‑68 (MQ=255)
                                                            gCCAAGAAGCTGGAGAAAGCGCGCAGTCTGCTCATTTCAACCACCTTGTCGATCAATGAGATATCGCa  >  1:4322949/1‑68 (MQ=255)
                                                            gCCAAGAAGCGGGAGAACGCGCGCAGTCTGCTGATTTCAACCCCCTTGTCGATCAATGAGATTTCGCa  >  1:1056625/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 < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: