Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,145,726 G→A intergenic (+57/‑138) dnaG → / → rpoD DNA primase/RNA polymerase sigma factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,145,7260GA99.2% 959.6 / ‑1.7 234intergenic (+57/‑138)dnaG/rpoDDNA primase/RNA polymerase sigma factor
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (134/98);  total (134/100)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 242 total aligned reads displayed.

AGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAA  >  REL606/3145653‑3145798
                                                                          |                                                                        
aGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGc                                                                        <  2:5255766/76‑1 (MQ=255)
 gCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCa                                                                       >  2:23468/1‑76 (MQ=255)
  cTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAg                                                                      >  1:3654595/1‑76 (MQ=255)
   tGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAagcagc                                                                     >  2:2400860/1‑76 (MQ=255)
   tGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAagcagc                                                                     >  2:4330530/1‑76 (MQ=255)
    ggCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCg                                                                    <  1:2468674/76‑1 (MQ=255)
    ggCAAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCg                                                                    <  2:2232647/76‑1 (MQ=255)
     gCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCgg                                                                   <  1:5328996/76‑1 (MQ=255)
     gCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCgg                                                                   <  1:4921734/76‑1 (MQ=255)
      cGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGc                                                                  >  2:3751243/1‑76 (MQ=255)
       gAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCTGCACTGAGAAGCAGCGGCa                                                                 <  1:62953/76‑1 (MQ=255)
       gAAAAAGTGATTTAACGGCTAAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCa                                                                 <  1:4543009/76‑1 (MQ=255)
        aaaaaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCaa                                                                <  1:2618604/76‑1 (MQ=255)
          aaaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAat                                                              >  2:6680016/1‑76 (MQ=255)
           taGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAata                                                             >  2:3866981/2‑76 (MQ=255)
            aGTGATTTAAGGGCTAAAGGGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAatat                                                            <  2:5458138/76‑1 (MQ=255)
            aGTGATTTAACGGCTTAAGTGCCGAAGAGCGACCGAGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAatat                                                            <  1:4093456/76‑1 (MQ=255)
              tGATTAAACGGCTTAGGTGCGGAAGAGCGATCGGAAACCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAatatat                                                          <  2:2114596/76‑1 (MQ=255)
               gATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAAtatata                                                         >  1:696293/1‑76 (MQ=255)
               gATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAAtatata                                                         <  1:737694/76‑1 (MQ=255)
               gATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAAtatata                                                         <  1:4218325/76‑1 (MQ=255)
                aTTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCATTGAGAAGCAGCGGCAAATATATaa                                                        >  2:6844240/1‑76 (MQ=255)
                 tttAACGTCTAAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAg                                                       <  2:6427431/76‑1 (MQ=255)
                 tttAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAg                                                       <  1:5450898/76‑1 (MQ=255)
                 tttAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAg                                                       >  1:3587675/1‑76 (MQ=255)
                  ttAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGt                                                      >  1:3850608/1‑76 (MQ=255)
                   tAACGGATAAATTGCAGAAACGCGCTCGGTACTCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTa                                                     <  2:6812400/76‑1 (MQ=255)
                    aaCGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTAc                                                    >  1:2751892/1‑76 (MQ=255)
                     aCGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACg                                                   >  2:1875093/1‑76 (MQ=255)
                        gCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGccc                                                >  2:900577/1‑76 (MQ=255)
                         cTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCt                                               >  2:1555936/1‑76 (MQ=255)
                         cTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCt                                               >  1:2517298/1‑76 (MQ=255)
                          ttAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTc                                              >  2:1905724/1‑76 (MQ=255)
                          ttAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTc                                              >  1:3156064/1‑76 (MQ=255)
                          ttAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTc                                              >  2:1292321/1‑76 (MQ=255)
                          taaaGTCCGGAAGGGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTc                                              <  2:6447251/74‑1 (MQ=255)
                           tAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCg                                             >  2:3229957/1‑76 (MQ=255)
                           tAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCg                                             >  1:986654/1‑76 (MQ=255)
                           tAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCg                                             >  1:2159087/1‑76 (MQ=255)
                           tAAGTGCAGAGGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCg                                             <  1:815545/76‑1 (MQ=255)
                            aggTGCAGAAGAGCGTCGGGAATGCCCCCGCCGGCACACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTa                                           <  2:3898179/74‑1 (MQ=255)
                            aaGTGCGGAAGAGGAACGGGGACGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGt                                            <  1:1162199/76‑1 (MQ=255)
                            aaGTGCGAAAGAGGGATCGGAAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGt                                            <  2:297911/76‑1 (MQ=255)
                            aaGTGCCGAAGGCCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGt                                            <  2:4733920/76‑1 (MQ=255)
                            aaGTCCCGAAGAGGGACCGGAAAGCCCCCGCCAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGt                                            <  2:1233585/76‑1 (MQ=255)
                             aGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTa                                           >  2:6550324/1‑76 (MQ=255)
                             aGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTa                                           >  2:5188603/1‑76 (MQ=255)
                               tGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAAt                                         >  2:1715690/1‑76 (MQ=255)
                                 ccGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTa                                       >  2:294897/1‑76 (MQ=255)
                                 ccGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTa                                       >  1:3794121/1‑76 (MQ=255)
                                 ccGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTa                                       >  2:1348692/1‑76 (MQ=255)
                                   gAAGAGCGATGGGGACCCCCCCGCCAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATc                                     <  2:59180/76‑1 (MQ=255)
                                   gAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATc                                     >  1:5226278/1‑76 (MQ=255)
                                    aaGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCg                                    >  1:6608471/1‑76 (MQ=255)
                                      gagCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGtt                                  >  2:3595482/1‑76 (MQ=255)
                                      gagCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGtt                                  >  2:4630430/1‑76 (MQ=255)
                                       ggCGAGCGGAAACCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTg                                 <  2:5093544/75‑1 (MQ=255)
                                         cGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGc                               >  1:4697528/1‑76 (MQ=255)
                                          gATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCg                              >  2:5302721/1‑76 (MQ=255)
                                          gATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCg                              >  2:5155540/1‑76 (MQ=255)
                                          gATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCg                              >  1:908972/1‑76 (MQ=255)
                                             ggggAACCCCCCGCCAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTa                           <  2:6177199/75‑1 (MQ=255)
                                             cGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTa                           >  2:6231894/1‑76 (MQ=255)
                                              gggAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTaa                          >  2:5610310/1‑76 (MQ=255)
                                              gggAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTaa                          >  1:3500038/1‑76 (MQ=255)
                                                gAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaac                        >  2:330442/1‑76 (MQ=255)
                                                    gCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaacaac                     <  1:5055995/76‑1 (MQ=255)
                                                    gCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaacaac                     <  1:5013474/76‑1 (MQ=255)
                                                      ccccGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCg                   >  1:4332200/1‑76 (MQ=255)
                                                      ccccGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCg                   >  2:6044137/1‑76 (MQ=255)
                                                       cccGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGt                  >  2:4876281/1‑76 (MQ=255)
                                                        cggACAGCCGCACTGAGAAGCAGCGCCAAATATATAAGTGCGCCCTCGTAATTATCGTTGGCGGTAAACAACCGtt                 <  2:3689646/74‑1 (MQ=255)
                                                        ccGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGtt                 >  1:879582/1‑76 (MQ=255)
                                                          gccAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTgg               <  1:6025628/74‑1 (MQ=255)
                                                          gccAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTgg               <  2:4940395/74‑1 (MQ=255)
                                                          gACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTgg               <  2:886253/76‑1 (MQ=255)
                                                          gACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTgg               <  1:5966941/76‑1 (MQ=255)
                                                            cAGCCGCACTGAGAAGCAGCGGCAAATATATAAGAACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGAt             >  2:312524/1‑76 (MQ=255)
                                                             aGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGAtt            >  2:6705943/1‑76 (MQ=255)
                                                              gCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGAttt           <  1:1862497/76‑1 (MQ=255)
                                                               ccGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTc          >  2:2455799/1‑76 (MQ=255)
                                                                tgCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCa         >  2:352535/2‑76 (MQ=255)
                                                                cGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCa         <  1:2799703/76‑1 (MQ=255)
                                                                cGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCa         >  1:6405326/1‑76 (MQ=255)
                                                                cGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCa         >  2:6203499/1‑76 (MQ=255)
                                                                 gCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAg        <  2:4231183/76‑1 (MQ=255)
                                                                 gCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAg        <  1:1324664/76‑1 (MQ=255)
                                                                 gCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAg        <  2:1115185/76‑1 (MQ=255)
                                                                  cACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGc       >  1:6213851/1‑76 (MQ=255)
                                                                  cACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGc       >  2:3331298/1‑76 (MQ=255)
                                                                   aCTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCg      >  2:1798345/1‑76 (MQ=255)
                                                                   aCTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCg      >  1:6479872/1‑76 (MQ=255)
                                                                    cTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGt     >  2:1967199/1‑76 (MQ=255)
                                                                    cTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGt     >  1:6807559/1‑76 (MQ=255)
                                                                    cTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGt     >  2:5664617/1‑76 (MQ=255)
                                                                    cTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGt     >  2:3926254/1‑76 (MQ=255)
                                                                     tGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGtt    >  1:3194452/1‑76 (MQ=255)
                                                                     tGAGAAGCAGCGGCAAATATATAAGTACACCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGtt    >  1:4941186/1‑76 (MQ=255)
                                                                      gagaAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTa   >  2:3071978/1‑76 (MQ=255)
                                                                       agaaGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTaa  >  1:1410097/1‑76 (MQ=255)
                                                                          |                                                                        
AGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAA  >  REL606/3145653‑3145798

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