Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,729,405 C→G 8.8% intergenic (‑34/+38) xylB ← / ← xylA xylulokinase/xylose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,729,4050CG 8.8% 153.7 / 9.1 57intergenic (‑34/+38)xylB/xylAxylulokinase/xylose isomerase
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base G (2/3);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-02

GCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGG  >  NC_000913/3729263‑3729540
                                                                                                                                               |                                                                                                                                       
gccggCGCGCAACCGTCAGCTGTTCCGTTCGCGCAGCAACCACCTCACCCTGCTCGTTGCGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGc                                                                                                                                         <  2:437259/140‑1 (MQ=255)
   ggCGCGAAGCGGTCAGCTTTTCCGTTCGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCGCGGTATCGCTAc                                                                                                                                      <  1:506321/143‑1 (MQ=255)
         aaaCGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGAGCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATaa                                                                                                                                >  1:93074/1‑143 (MQ=255)
             ggTCAGCTTGTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGATaccgg                                                                                                                             <  1:577824/143‑5 (MQ=255)
             ggTCAGCTTCTCCGTCTGCGCAGCAACCACCTCACGCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGATaccgg                                                                                                                             <  1:397808/143‑5 (MQ=255)
             gctCATCTTTGCCGTTGGCACACCATCCCCCACACTCTTCGCGTTTAGCAATATAACTTTTACGCCCGAGGCCGCAGGATCTAGTTCGTGATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGATaccgg                                                                                                                             <  2:26327/141‑5 (MQ=255)
               tCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAGGATCTATCCCGATATACATATCGGGAGTTACTTACCAAAATGCTCGCTATCGCTACCGATAACCggtc                                                                                                                          >  1:77138/1‑141 (MQ=255)
               tCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGGTATACATATCGATCGTTCCCTAAAAAAATGCTCGGTGTCGCTACCGATAActcggc                                                                                                                          >  2:449417/1‑138 (MQ=255)
                cAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCAAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGACACCCGGGcc                                                                                                                         >  2:598917/1‑143 (MQ=255)
                   cTTTTTCTTTTGCGCAGCGACCCCCTCCCCCTGCTCGTTGTGCAAATTCACTTTTCCGCCCGTGGTGCCGCGATCTATCCCGTTATACATATCGATCGTTCCTTAATCGAATGCCCGGTATCGCTACCGATCACCGGGCCacc                                                                                                                      <  2:388166/143‑3 (MQ=255)
                           gttGCGCAGCAACCACCTCACCGTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGca                                                                                                              <  2:326245/142‑1 (MQ=255)
                             tGCGCAGCAACCACCTCAGCCTGCTCGTTGAGCAAAATATCTTCTNCGCCCGAGGTGCCAAGATCTATCCCGCTANACATATCGATCGTTCTTTNAAAAAANGCCCGGNATCGCTACCGATAACCGGGCCAACGGACTGcaca                                                                                                            <  2:538135/143‑1 (MQ=255)
                                gcagcaACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTCTCTTGGTTTTCACAACAAAAGAACGTAGGACTAGAGCTCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                                                                         >  2:190990/1‑143 (MQ=255)
                                         accTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAt                                                                                                <  1:595790/143‑1 (MQ=255)
                                           cTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAttt                                                                                              <  2:751165/143‑1 (MQ=255)
                                             cgcgcatcTCGTTGACCAAAATATGTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATGTCGATCGTTCCTTACAAAGGTGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGt                                                                                            <  2:599567/136‑1 (MQ=255)
                                             cACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGt                                                                                            <  1:108471/143‑1 (MQ=255)
                                                   gcgcGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTAGAACa                                                                                      <  2:450271/140‑1 (MQ=255)
                                                   gCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACa                                                                                      <  1:152649/143‑1 (MQ=255)
                                                       gTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTAACGATAACCGGGCCAACGGACTGCACAGTCAGCCGTTATTTGTCGAACAGATa                                                                                  >  2:390761/1‑143 (MQ=255)
                                                        ttGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAGCAGAtga                                                                                 >  1:508667/1‑141 (MQ=255)
                                                        ttGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATaa                                                                                 >  1:427607/1‑143 (MQ=255)
                                                          gAGCAATCTATCTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGAGCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATg                                                                               <  2:464787/143‑1 (MQ=255)
                                                          gAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATg                                                                               <  2:285042/143‑1 (MQ=255)
                                                           aGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAGAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATgg                                                                              >  2:396934/1‑143 (MQ=255)
                                                                 aTAACTTTTTCGCCCGAGGTGCCAAGCTCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTAcc                                                                        <  2:564307/143‑1 (MQ=255)
                                                                  tAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAAACGGGCACACTGACTGCACAGTTACCTGTTATTTGTCGAACCGATCATGCTTCActa                                                                       >  2:594865/1‑141 (MQ=255)
                                                                    aCTTCTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGa                                                                     <  1:891608/143‑1 (MQ=255)
                                                                     cTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCCGAt                                                                    >  1:614126/1‑143 (MQ=255)
                                                                      ttttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGAtt                                                                   >  1:326025/1‑143 (MQ=255)
                                                                      ttttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGAtg                                                                   >  1:227318/1‑142 (MQ=255)
                                                                      ttttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGACGTTATTTGTCGAACAGATAATGGTTTACCAGaat                                                                   >  2:201258/1‑141 (MQ=255)
                                                                            gCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCTGATTTTcctg                                                             >  1:153940/1‑141 (MQ=255)
                                                                             cccGAGGTGCCAAGATCTATCCAGATATACATATCGATCGTTCATTAAAAAAATGCCCGGTATCGCTACAGATAGCCGGGCCAACGGACTGCACAGTTAGCCGTTATGTGTCGCACAGATAAGGGTTTACAAGGTTTTCCAgc                                                            >  2:518870/1‑142 (MQ=255)
                                                                                    tGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCt                                                     >  1:836258/1‑143 (MQ=255)
                                                                                      ccAAGCTCTATCCCGATATACATATCGATCGTTCCTTAAAACAATGCCCGGTATCGCTACCGATACCCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGACTCTTCAGTTGTTCCTgg                                                   <  2:502591/143‑1 (MQ=255)
                                                                                        aaGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGACTTTCCAGTTGTcccagcca                                                 >  2:537100/1‑138 (MQ=255)
                                                                                            tCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAGCAGATAATAGTGTACCAGAGTTTCAAGTTGTTCCTCGCGACCa                                             >  1:262725/1‑143 (MQ=255)
                                                                                             cTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCAc                                            >  1:266059/1‑143 (MQ=255)
                                                                                               aTCCCGATCTACCTATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCActc                                          <  1:879517/143‑1 (MQ=255)
                                                                                                 cccGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTg                                        <  2:574921/143‑1 (MQ=255)
                                                                                                        taCAGATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCAcc                                 <  2:721787/143‑1 (MQ=255)
                                                                                                            tatCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGATATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGgag                             >  1:652026/1‑143 (MQ=255)
                                                                                                             atcgatcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTATCATGCACCgccga                            >  2:729539/1‑139 (MQ=255)
                                                                                                              tcgatcgTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCCCTCTGATGCACCGcagtc                           >  1:731924/1‑141 (MQ=255)
                                                                                                                    gTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAatga                     <  1:583580/143‑1 (MQ=255)
                                                                                                                       ccTTAAAAAAACGCCCGGTATCGCTACCGGTAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGt                  <  2:77138/143‑1 (MQ=255)
                                                                                                                        cTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGtt                 >  1:537085/1‑143 (MQ=255)
                                                                                                                         ttAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTc                <  2:280855/143‑1 (MQ=255)
                                                                                                                              aaaaTGCCGGGTAACGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGGTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAg           <  1:303371/143‑1 (MQ=255)
                                                                                                                              aaaaTGCCCGGTATCGCTACCGATAACCCGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACCGTTTATGGTTTACCAGGTTTTCCAGTTGTTCCTGGCGGCCGCGCGGATGGACCGGACACCGATGTTGTGCCTGAg           >  1:37869/1‑143 (MQ=255)
                                                                                                                                   gCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGCTTCCCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGACGCACCGGCGACAAATGCTGTTCCTGaccccat      >  2:913762/1‑137 (MQ=255)
                                                                                                                                    cccGGTTATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATAtt     >  1:552032/1‑143 (MQ=255)
                                                                                                                                      cggttATCGGTAGCGATACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGCGCATATTTg   >  1:778037/5‑143 (MQ=255)
                                                                                                                                      cGGTATCGCTACCGATATCCGGGCCAACGGTCTTTACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTGACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTg   <  1:390761/143‑1 (MQ=255)
                                                                                                                                       ggTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACACATGATCTTCCTGAGCATActtcg  >  1:554402/1‑141 (MQ=255)
                                                                                                                                       ggTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTgg  <  2:525803/143‑1 (MQ=255)
                                                                                                                                       ggTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTgg  <  2:326025/143‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                       
GCGGGCGCGAAACGGTCAGCTTTTCCGTTTGCGCAGCAACCACCTCACCCTGCTCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGAGACAAATGATGTTCCTGAGCATATTTGG  >  NC_000913/3729263‑3729540

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