Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,147,762 G→A intergenic (+57/+22) rpoD → / ← ygjF RNA polymerase sigma factor/G/U mismatch‑specific DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,147,7620GA99.6% 856.2 / NA 224intergenic (+57/+22)rpoD/ygjFRNA polymerase sigma factor/G/U mismatch‑specific DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (124/99);  total (124/100)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

GCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAAC  >  REL606/3147689‑3147834
                                                                         |                                                                        
 cTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCa                                                                       >  1:7581651/1‑76 (MQ=255)
 cTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCa                                                                       >  1:1595415/1‑76 (MQ=255)
 cTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCa                                                                       >  2:7937989/1‑76 (MQ=255)
 cTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCa                                                                       >  2:1884449/1‑76 (MQ=255)
 cTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCa                                                                       <  2:2382372/76‑1 (MQ=255)
    ccTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCt                                                                    >  1:1345059/1‑76 (MQ=255)
    ccTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCt                                                                    <  1:5609141/76‑1 (MQ=255)
    ccTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCt                                                                    <  1:351655/76‑1 (MQ=255)
     cTGGACGATTAATGGGTGGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCtt                                                                   <  2:7608406/76‑1 (MQ=38)
     cTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCtt                                                                   <  2:598227/76‑1 (MQ=255)
     cTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCtt                                                                   <  2:554836/76‑1 (MQ=255)
       ggACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCtttt                                                                 <  1:5797369/76‑1 (MQ=255)
       ggACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCtttt                                                                 <  2:3295889/76‑1 (MQ=255)
         aCGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCg                                                               >  2:3616563/1‑76 (MQ=255)
         aCGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCg                                                               <  1:812799/76‑1 (MQ=255)
           gATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGt                                                             <  2:4330651/76‑1 (MQ=255)
           gATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGt                                                             <  1:2566486/76‑1 (MQ=255)
            aTTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTg                                                            >  2:1987785/1‑76 (MQ=255)
            aTTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTg                                                            <  2:2195573/76‑1 (MQ=255)
            aTTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTg                                                            <  1:2216077/76‑1 (MQ=255)
             ttAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTgg                                                           >  2:6265684/1‑76 (MQ=255)
               aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCg                                                         >  2:6870347/1‑76 (MQ=255)
               aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCg                                                         >  2:129823/1‑76 (MQ=255)
               aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCg                                                         >  2:2234818/1‑76 (MQ=255)
               aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAAGCCCTCTGCACAAACACCACCTTTTCGGTGGCg                                                         >  2:6757966/1‑76 (MQ=255)
                aTCGGTAGGCCGGATGAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGt                                                        <  1:5384646/76‑1 (MQ=255)
                 tCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCATTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGtt                                                       <  1:6977498/76‑1 (MQ=255)
                 tCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGtt                                                       <  1:6126972/76‑1 (MQ=255)
                  cGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGttt                                                      >  2:8042164/1‑76 (MQ=255)
                   ggTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGtttt                                                     >  2:4028580/1‑76 (MQ=255)
                    gTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGttttt                                                    <  2:4454702/76‑1 (MQ=255)
                     tAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGtttttt                                                   >  2:507959/1‑76 (MQ=255)
                     tAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGtttttt                                                   >  2:4129279/1‑76 (MQ=255)
                     tAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGtttttt                                                   <  1:7116309/76‑1 (MQ=255)
                      aGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGTCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTa                                                  >  1:3960016/1‑76 (MQ=255)
                      aGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTa                                                  <  2:2746957/76‑1 (MQ=255)
                       ggCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTAt                                                 <  2:2424633/76‑1 (MQ=255)
                        gCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATc                                                >  2:2754724/1‑76 (MQ=255)
                         ccGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCg                                               <  2:3550955/76‑1 (MQ=255)
                         ccGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCg                                               >  1:4133828/1‑76 (MQ=255)
                          cGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGc                                              >  1:1188641/1‑76 (MQ=255)
                          cGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGc                                              >  2:1627794/1‑76 (MQ=255)
                          cGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGc                                              >  1:3183480/1‑76 (MQ=255)
                              tCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCAc                                          >  1:989422/1‑76 (MQ=255)
                               cAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACg                                         <  2:4515612/76‑1 (MQ=255)
                               cAGGCGTTAACCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGAGGCGTTCTTTACCGACCAcc                                         >  2:663992/1‑75 (MQ=255)
                                aGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGc                                        >  2:82057/1‑76 (MQ=255)
                                aGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGc                                        >  2:782302/1‑76 (MQ=255)
                                  gCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCAc                                      <  1:2099638/76‑1 (MQ=255)
                                  gCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCAc                                      >  1:6042070/1‑76 (MQ=255)
                                   cGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACt                                     >  2:495875/1‑76 (MQ=255)
                                   cGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACt                                     <  2:2903050/76‑1 (MQ=255)
                                   cGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACt                                     >  1:676350/1‑76 (MQ=255)
                                    gTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTa                                    <  1:4001703/76‑1 (MQ=255)
                                    gTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTa                                    >  2:7016821/1‑76 (MQ=255)
                                    gTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTa                                    >  1:5855929/1‑76 (MQ=255)
                                     ttACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTAc                                   <  2:8621/76‑1 (MQ=255)
                                     ttACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTAc                                   <  2:2898624/76‑1 (MQ=255)
                                     ttACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTAc                                   >  1:7085755/1‑76 (MQ=255)
                                      tACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTAcc                                  >  1:2398066/1‑76 (MQ=255)
                                      tACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTAcc                                  <  2:7253110/76‑1 (MQ=255)
                                      tACGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTAcc                                  >  2:4348115/1‑76 (MQ=255)
                                      tACGCCGCACCCGGCACTAGGCCCCCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTAcc                                  <  1:1090028/76‑1 (MQ=255)
                                       aCGCCGCACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCa                                 <  2:183708/76‑1 (MQ=255)
                                        cgccgcACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAg                                <  2:1687198/76‑1 (MQ=255)
                                            gcACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGcc                            >  1:7936678/1‑76 (MQ=255)
                                             cACCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCt                           <  1:6287924/76‑1 (MQ=255)
                                              ccccGGCACTAGGCCCTCTGCACAAACACCCCCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTg                          <  1:6298445/75‑1 (MQ=255)
                                              aCCCGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTg                          <  1:5000092/76‑1 (MQ=255)
                                               cccGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTgg                         <  2:1535161/76‑1 (MQ=255)
                                                ccGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGt                        >  2:2141876/1‑76 (MQ=255)
                                                ccGGCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGt                        >  1:5576188/1‑76 (MQ=255)
                                                  ggCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTcc                      <  1:475723/76‑1 (MQ=255)
                                                   gCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCa                     <  2:3730918/76‑1 (MQ=255)
                                                   gCACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCa                     >  2:1069941/1‑76 (MQ=255)
                                                    cACTAGTCCCTCTGCCAAACCGCCTCCCTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAg                    <  2:4748127/76‑1 (MQ=255)
                                                    cACTAGGCCCTCTGCACAAGCACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAg                    <  1:2007961/76‑1 (MQ=255)
                                                    cACTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCTCCTGGTCCAg                    <  2:4617474/76‑1 (MQ=255)
                                                     aCTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGc                   >  1:8065097/1‑76 (MQ=255)
                                                     aCTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGc                   <  1:5499530/76‑1 (MQ=255)
                                                      cTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCt                  >  1:7187680/1‑76 (MQ=255)
                                                      cTAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCt                  >  2:6431245/1‑76 (MQ=255)
                                                       tAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTc                 <  1:3647024/76‑1 (MQ=255)
                                                       tAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTc                 >  1:87670/1‑76 (MQ=255)
                                                       tAGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTc                 >  1:5177778/1‑76 (MQ=255)
                                                        aGGCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTcg                >  1:5443090/1‑76 (MQ=255)
                                                          gCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTcgcg              <  1:6411223/76‑1 (MQ=255)
                                                          gCCCTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTcgcg              <  2:3050572/76‑1 (MQ=255)
                                                           cccTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGa             >  2:5146721/1‑76 (MQ=255)
                                                           cccTCTGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGa             >  2:738344/1‑76 (MQ=255)
                                                              tctGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATAc          >  2:2508829/1‑76 (MQ=255)
                                                               ctGCACAAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACg         >  1:4530866/1‑76 (MQ=255)
                                                                  cacaAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCtt      >  2:3174980/1‑76 (MQ=255)
                                                                   acaAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTc     <  1:3940796/76‑1 (MQ=255)
                                                                   acaAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTc     <  2:5455159/76‑1 (MQ=255)
                                                                   acaAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTc     >  1:7433870/1‑76 (MQ=255)
                                                                    caAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCa    <  2:6027163/76‑1 (MQ=255)
                                                                    caAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCa    <  1:2943955/76‑1 (MQ=255)
                                                                    caAACACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCa    >  1:1967811/1‑76 (MQ=255)
                                                                      aaCACCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAAc  <  1:1780425/76‑1 (MQ=255)
                                                                         |                                                                        
GCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAAC  >  REL606/3147689‑3147834

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