Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,762,428 C→T intergenic (‑16/+47) rpoS ← / ← nlpD RNA polymerase sigma factor/predicted outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,762,4280CT99.3% 1241.3 / ‑2.6 295intergenic (‑16/+47)rpoS/nlpDRNA polymerase sigma factor/predicted outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (155/138);  total (155/140)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 301 total aligned reads displayed.

TCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCG  >  REL606/2762355‑2762500
                                                                         |                                                                        
tctcATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGt                                                                        >  2:4507313/1‑76 (MQ=255)
 ctcATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCtgtg                                                                       >  2:3008847/1‑76 (MQ=255)
 ctcATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCtgtg                                                                       >  2:6657127/1‑76 (MQ=255)
   catcaAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGAt                                                                     >  1:2067022/1‑76 (MQ=255)
     tcaAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATcc                                                                   >  1:1332019/1‑76 (MQ=255)
      caAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATccc                                                                  <  1:3603987/76‑1 (MQ=255)
      caAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATccc                                                                  <  1:4561796/76‑1 (MQ=255)
      caAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATccc                                                                  >  2:6747170/1‑76 (MQ=255)
      caAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATccc                                                                  <  1:2469573/76‑1 (MQ=255)
       aaaTTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCt                                                                 >  2:6952629/1‑76 (MQ=255)
        aaTTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCtt                                                                <  2:6036880/76‑1 (MQ=255)
         aTTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTg                                                               <  1:5870146/76‑1 (MQ=255)
          ttCCGCATCTTCATTTAAATCATGACCTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGa                                                              <  1:5125660/76‑1 (MQ=255)
           tCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGAc                                                             <  2:7241937/76‑1 (MQ=255)
             cGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACgg                                                           >  1:1176565/1‑76 (MQ=255)
             cGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACgg                                                           <  2:5651892/76‑1 (MQ=255)
              gCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGa                                                          >  2:7021288/1‑76 (MQ=255)
               cATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGaa                                                         >  1:7277117/1‑76 (MQ=255)
               cATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGaa                                                         >  2:7659723/1‑76 (MQ=255)
                aTCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAAc                                                        >  1:1010754/1‑76 (MQ=255)
                aTCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAAc                                                        >  2:3598612/1‑76 (MQ=255)
                aTCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAAc                                                        >  2:6571428/1‑76 (MQ=255)
                aTCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAAc                                                        >  2:2171523/1‑76 (MQ=255)
                 tCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACa                                                       <  2:1084809/76‑1 (MQ=255)
                  cTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACAt                                                      <  1:6445701/76‑1 (MQ=255)
                  cTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACAt                                                      <  1:1924013/76‑1 (MQ=255)
                   ttCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACAtt                                                     <  2:3269712/76‑1 (MQ=255)
                   ttCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACAtt                                                     >  2:4483226/1‑76 (MQ=255)
                   ttCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACAtt                                                     >  2:2870056/1‑76 (MQ=255)
                    tCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTc                                                    >  2:1332444/1‑76 (MQ=255)
                     cATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCa                                                   >  1:3765420/1‑76 (MQ=255)
                     cATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCa                                                   <  2:434061/76‑1 (MQ=255)
                     cATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCa                                                   <  1:5049654/76‑1 (MQ=255)
                      aTTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCaa                                                  <  1:2645657/76‑1 (MQ=255)
                        ttAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGc                                                >  1:413204/1‑76 (MQ=255)
                        ttAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGc                                                >  2:6927975/1‑76 (MQ=255)
                         tAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCa                                               >  1:6746132/1‑76 (MQ=255)
                           aaTCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCa                                                   <  1:6224400/70‑1 (MQ=255)
                           aaTCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCa                                                   <  1:1372302/70‑1 (MQ=255)
                           aaTCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCa                                                   >  2:1372302/1‑70 (MQ=255)
                           aaTCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCaaa                                             >  2:7116112/1‑76 (MQ=255)
                           aaTCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCaaa                                             <  2:2288422/76‑1 (MQ=255)
                            aTCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCaaaa                                            <  1:615007/76‑1 (MQ=255)
                            aTCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCaaaa                                            <  1:950363/76‑1 (MQ=255)
                            aTCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCaaaa                                            >  1:2316460/1‑76 (MQ=255)
                              cATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGAAACATTCAAGCAAAAGc                                          >  2:5909562/1‑76 (MQ=255)
                               aTGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGcc                                         <  2:7802253/76‑1 (MQ=255)
                                 gAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTg                                       <  2:6794043/76‑1 (MQ=255)
                                 gAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTg                                       <  1:3709480/76‑1 (MQ=255)
                                 gAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTg                                       <  2:2024118/76‑1 (MQ=255)
                                   aCTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGt                                     <  1:4202131/76‑1 (MQ=255)
                                   aCTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGt                                     <  1:3844764/76‑1 (MQ=255)
                                     tttCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTc                                   >  2:4069046/1‑76 (MQ=255)
                                      ttCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTcc                                  >  2:6298290/1‑76 (MQ=255)
                                      ttCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTcc                                  <  2:4948965/76‑1 (MQ=255)
                                      ttCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTcc                                  >  1:1976550/1‑76 (MQ=255)
                                       tCAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTccg                                 <  1:3944566/76‑1 (MQ=255)
                                        cAGCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATCCAAGCAAAAGCCTGGTTccgc                                <  1:5600489/76‑1 (MQ=255)
                                          gCGTATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTccgccg                              <  1:1103278/76‑1 (MQ=255)
                                             tATTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAtt                           >  1:2846336/1‑76 (MQ=255)
                                              aTTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAttt                          <  2:4236478/76‑1 (MQ=255)
                                              aTTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAttt                          >  2:5352175/1‑76 (MQ=255)
                                              aTTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAttt                          <  1:3080731/76‑1 (MQ=255)
                                              aTTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAttt                          >  1:7388579/1‑76 (MQ=255)
                                              aTTCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAttt                          >  1:6686389/1‑76 (MQ=255)
                                               ttCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTa                         <  1:487272/76‑1 (MQ=255)
                                                tCTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTAt                        >  1:5390207/1‑76 (MQ=255)
                                                 cTGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATc                       >  1:7048591/1‑76 (MQ=255)
                                                  tGACTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCTTGGTTCCGCCGATTTATCg                      >  1:7866369/1‑76 (MQ=255)
                                                     cTCATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTg                   >  2:3701069/1‑76 (MQ=255)
                                                       cATAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATACAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCg                 >  1:474023/1‑76 (MQ=255)
                                                         tAAGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGc               >  2:462460/1‑76 (MQ=255)
                                                           aGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCaa             >  1:6302475/1‑76 (MQ=255)
                                                           aGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCaa             <  1:7183196/76‑1 (MQ=255)
                                                           aGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCaa             <  1:5300275/76‑1 (MQ=255)
                                                           aGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCaa             <  1:7599559/76‑1 (MQ=255)
                                                           aGGTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCaa             >  1:3983557/1‑76 (MQ=255)
                                                            ggTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCaaa            <  2:539036/76‑1 (MQ=255)
                                                             gTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAAt           <  2:6183013/76‑1 (MQ=255)
                                                             gTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAAt           >  1:7520206/1‑76 (MQ=255)
                                                             gTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAAt           <  1:4319732/76‑1 (MQ=255)
                                                             gTGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAACCTGGTTCCGCCGATTTATCGCTGCGGCaaaa           >  2:1439617/1‑75 (MQ=255)
                                                              tGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATa          <  2:387112/76‑1 (MQ=255)
                                                              tGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATa          <  1:4507313/76‑1 (MQ=255)
                                                              tGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATa          <  1:4634988/76‑1 (MQ=255)
                                                              tGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATa          <  1:5950278/76‑1 (MQ=255)
                                                              tGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATa          <  1:6491678/76‑1 (MQ=255)
                                                              tGGCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATa          <  1:7814713/76‑1 (MQ=255)
                                                                gCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTccgc                                >  2:124619/1‑52 (MQ=255)
                                                                gCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTccgc                                <  1:124619/52‑1 (MQ=255)
                                                                gCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAAc        <  2:5840077/76‑1 (MQ=255)
                                                                gCTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAAc        <  1:1563446/76‑1 (MQ=255)
                                                                 cTCCTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACg       <  1:5721999/76‑1 (MQ=255)
                                                                   ccTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCa     >  1:7722204/1‑76 (MQ=255)
                                                                   ccTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCa     >  1:847623/1‑76 (MQ=255)
                                                                   ccTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCa     <  2:4810209/76‑1 (MQ=255)
                                                                    cTACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAg    <  1:3343616/76‑1 (MQ=255)
                                                                     tACCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGc   <  2:5973274/76‑1 (MQ=255)
                                                                      aCCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCg  >  1:6573985/1‑76 (MQ=255)
                                                                      aCCTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCg  >  2:7976682/1‑76 (MQ=255)
                                                                         |                                                                        
TCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCG  >  REL606/2762355‑2762500

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