Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,145,830 A→G intergenic (+161/‑34) 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,8300AG99.4% 684.8 / ‑4.1 176intergenic (+161/‑34)dnaG/rpoDDNA primase/RNA polymerase sigma factor
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (77/98);  total (77/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

Only 100 of 193 total aligned reads displayed.

AATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTC  >  REL606/3145757‑3145902
                                                                         |                                                                        
aaTTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCCACGcc                                                                        >  2:2155495/1‑76 (MQ=255)
aaTTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGcc                                                                        >  1:4102227/1‑76 (MQ=255)
aaTTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGcc                                                                        >  1:3330448/1‑76 (MQ=255)
 aTTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCCACGCCa                                                                       >  2:2061101/1‑76 (MQ=255)
 aTTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCa                                                                       >  1:1432641/1‑76 (MQ=255)
 aTTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAAACGCCCAACGccc                                                                       >  2:4858928/1‑75 (MQ=255)
  ttATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCaa                                                                      <  1:4440897/76‑1 (MQ=255)
  ttATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCaa                                                                      >  1:5263777/1‑76 (MQ=255)
  ttATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCaa                                                                      <  1:6057170/76‑1 (MQ=255)
  ttATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCaa                                                                      >  1:2440728/1‑76 (MQ=255)
   tATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTGAATCGCCCAAAGCCAAc                                                                     >  1:3128815/1‑76 (MQ=255)
   tATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAAc                                                                     <  1:1273400/76‑1 (MQ=255)
     tCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCt                                                                   >  1:3763594/1‑76 (MQ=255)
     tCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCt                                                                   >  1:5536758/1‑76 (MQ=255)
     tCGTTGGCGGTAAACAACCGTGGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCt                                                                   <  2:4254869/76‑1 (MQ=255)
      cGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTc                                                                  <  1:1095237/76‑1 (MQ=255)
      cGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTc                                                                  >  1:2197722/1‑76 (MQ=255)
       gTTGGCGTAAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGCCATCGGGTCAATCGCCCAACGCCAACCTCa                                                                 <  2:591615/76‑1 (MQ=255)
       gTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCa                                                                 <  2:2749286/76‑1 (MQ=255)
       gTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGGCAATCGCCCAACGCCAACCTCa                                                                 >  2:4193536/1‑76 (MQ=255)
        ttGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCCACGCCAACCTCAt                                                                >  1:1600447/1‑76 (MQ=255)
        ttGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCAt                                                                <  1:1092926/76‑1 (MQ=255)
         tGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATg                                                               <  1:196032/76‑1 (MQ=255)
         tGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATg                                                               <  1:3482419/76‑1 (MQ=255)
          ggCGGTAAACAACCGTTGGATTTCAGCGTTAACTGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGa                                                              <  1:757340/76‑1 (MQ=255)
          ggCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGa                                                              <  1:3700513/76‑1 (MQ=255)
            cGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGaaa                                                            >  1:122607/1‑76 (MQ=255)
            cGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGaaa                                                            <  1:4254002/76‑1 (MQ=255)
             ggTAAACAACCGTTGGATTTCAGCTTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAAt                                                           <  1:530774/76‑1 (MQ=255)
              gTAAACAGCCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATa                                                          <  2:49072/76‑1 (MQ=255)
              gTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATa                                                          >  1:5084496/1‑76 (MQ=255)
               tAAACAACCGTTGGATTTCAGCGTTAACGGCTGGAGGAGATCGGGGCCATCGCCCCACACCCAACACATGAAATaa                                                         >  2:4945919/1‑76 (MQ=255)
               tAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATaa                                                         <  2:1791645/76‑1 (MQ=255)
               tAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATaa                                                         <  2:3178052/76‑1 (MQ=255)
               tAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATaa                                                         <  2:3829271/76‑1 (MQ=255)
                aaaCAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAg                                                        >  2:539580/1‑76 (MQ=255)
                 aacaacCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAgt                                                       <  1:4662545/76‑1 (MQ=255)
                  acaacCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAgtg                                                      >  1:5156059/1‑76 (MQ=255)
                   caacCGTTGGATTTCAGCGTTAACGGCTGGAGGACATCGGGTGAATCGCCCAACGCCAACCTCATGAAATAAgtgt                                                     >  2:909837/1‑76 (MQ=255)
                   caacCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAgtgt                                                     <  2:3128815/76‑1 (MQ=255)
                   caacCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGACAATCGCCCAACGCCAACCTCATGAAATAAgggt                                                     >  1:3557684/1‑76 (MQ=255)
                    aacCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGtgtg                                                    <  1:2941390/76‑1 (MQ=255)
                     acCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTgg                                                   >  2:5344503/1‑76 (MQ=255)
                       cGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGAt                                                 >  1:2085013/1‑76 (MQ=255)
                       cGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGGCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGaa                                                 >  2:4729770/1‑75 (MQ=255)
                        gTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATa                                                <  2:778496/76‑1 (MQ=255)
                          tGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATAcc                                              >  1:3470072/1‑76 (MQ=255)
                          tGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATAcc                                              <  2:1358119/76‑1 (MQ=255)
                            gATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGt                                            <  2:3618815/76‑1 (MQ=255)
                            gATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGt                                            >  1:3638035/1‑76 (MQ=255)
                            gATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATAACGt                                            >  1:4835032/1‑76 (MQ=255)
                               ttCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACGGTCtt                                         <  2:3557684/76‑1 (MQ=255)
                               ttCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCtt                                         <  2:3579198/76‑1 (MQ=255)
                               ttCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCtt                                         <  1:4354304/76‑1 (MQ=255)
                                tCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTa                                        >  1:3196675/1‑76 (MQ=255)
                                tCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTa                                        <  1:3865693/76‑1 (MQ=255)
                                 cAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTAt                                       >  1:4743354/1‑76 (MQ=255)
                                   gCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATgg                                     <  1:1565729/76‑1 (MQ=255)
                                   gCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATgg                                     <  2:2592794/76‑1 (MQ=255)
                                    cGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATAACGTCTTAttgg                                    >  1:1633000/1‑73 (MQ=255)
                                     gTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAg                                   >  1:3038274/1‑76 (MQ=255)
                                     gTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAg                                   <  1:2142811/76‑1 (MQ=255)
                                      ttAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGc                                  <  1:862022/76‑1 (MQ=255)
                                      ttAACGGCTGAAGGACATCGGGCCAATCGCCCAACGCCAACCTCATAAAATAAGTGTGGATACCGTCTTTTGGAGc                                  >  1:267311/1‑76 (MQ=255)
                                      ttAACGGCTGAAAGACACGGGGTCAATCGGCCAACGCCAACCTCATGAAATAAGTGGGGAGACCCGCTTTTGGTGc                                  >  2:3302701/1‑76 (MQ=255)
                                       tAACGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCa                                 <  1:6261543/76‑1 (MQ=255)
                                        aaCGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCaa                                <  1:2541277/76‑1 (MQ=255)
                                         aCGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCaaa                               <  2:5475659/76‑1 (MQ=255)
                                         aCGGCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCaaa                               <  2:4361029/76‑1 (MQ=255)
                                           ggCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAAc                             <  1:4482323/76‑1 (MQ=255)
                                           ggCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAAc                             <  2:3588364/76‑1 (MQ=255)
                                           ggCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAAc                             <  1:6179355/76‑1 (MQ=255)
                                            gCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAAcc                            >  2:4619093/1‑76 (MQ=255)
                                            gCTGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCCTCTTATGGAGCAAAAcc                            >  2:1138978/1‑76 (MQ=255)
                                              tGAAGGACATCGGGTCAGTCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCg                          >  1:132667/1‑76 (MQ=255)
                                              tGAAGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCg                          >  2:398642/1‑76 (MQ=255)
                                                aaGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCa                        >  1:1811629/1‑76 (MQ=255)
                                                aaGGACAGTGGGTCAATCGCCCAACGCCAGCCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCa                        <  2:5614460/76‑1 (MQ=255)
                                                 aGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAg                       >  1:2556708/1‑76 (MQ=255)
                                                 aGGACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAg                       <  2:3328014/76‑1 (MQ=255)
                                                  ggACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGt                      >  1:57711/1‑76 (MQ=255)
                                                   gACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTc                     >  1:4756699/1‑76 (MQ=255)
                                                   gACATCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTc                     <  1:2049022/76‑1 (MQ=255)
                                                      aTCGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTcaca                  <  1:9703/76‑1 (MQ=255)
                                                        cGGGTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGGCACAGc                >  2:1609549/1‑76 (MQ=255)
                                                          ggTCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCCCAGTCACAGCTg              >  2:2140373/1‑76 (MQ=255)
                                                            tCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGaa            <  2:6281350/76‑1 (MQ=255)
                                                            tCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGaa            >  1:381630/1‑76 (MQ=255)
                                                            tCAATCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGaa            <  1:4193536/76‑1 (MQ=255)
                                                               aTCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCCGTCACAGCTGAAAct         >  1:2747108/1‑76 (MQ=255)
                                                               aTCGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAAct         <  1:1504932/76‑1 (MQ=255)
                                                                 cGCCCAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGGGCAAAAACCGCAGTCACAGCTGAAActtc       >  2:2800316/1‑76 (MQ=255)
                                                                   cccAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAActtctt     >  1:2038278/1‑76 (MQ=255)
                                                                   cccAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAActtctt     <  2:1966356/76‑1 (MQ=255)
                                                                   cccAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAActtctt     <  2:2555510/76‑1 (MQ=255)
                                                                    ccAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGGGCAAAACCCGCAGTCACAGCTGAAActtcttt    >  2:555271/1‑75 (MQ=255)
                                                                    ccAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTg    <  2:4785989/76‑1 (MQ=255)
                                                                     cAACGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGt   <  2:1169815/76‑1 (MQ=255)
                                                                      aaCGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTc  >  2:5206409/1‑76 (MQ=255)
                                                                      aaCGCCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTc  >  2:2946370/1‑76 (MQ=255)
                                                                         |                                                                        
AATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTC  >  REL606/3145757‑3145902

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