Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,558 A→G 6.2% intergenic (+1326/‑453) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,5580AG 6.2% 311.3 / 10.5 112intergenic (+1326/‑453)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (76/29);  new base G (4/3);  total (80/32)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 123 total aligned reads displayed.

GCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTA  >  minE/2572488‑2572627
                                                                      |                                                                     
gCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGaa                                                                       <  1:3579830/71‑1 (MQ=255)
 ccATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGaa                                                                       >  1:2520183/1‑70 (MQ=255)
      tATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGaa                                                                       >  1:3574421/1‑65 (MQ=255)
      tATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGaa                                                                       >  1:1919270/1‑65 (MQ=255)
      tATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGaa                                                                       >  1:421611/1‑65 (MQ=255)
      tATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAGTACt                                                                   >  1:571853/1‑69 (MQ=255)
      tATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAGTACt                                                                   >  1:1484245/1‑69 (MQ=255)
      tATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACt                                                                   >  1:2773236/1‑69 (MQ=255)
      tATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACt                                                                   >  1:3504144/1‑69 (MQ=255)
      tATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACt                                                                   >  1:4103727/1‑69 (MQ=255)
      tATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACt                                                                   >  1:552280/1‑69 (MQ=255)
      tATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACt                                                                   >  1:3381009/1‑69 (MQ=255)
      tATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACt                                                                   >  1:3340045/1‑69 (MQ=255)
      tATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACt                                                                   >  1:1169621/1‑69 (MQ=255)
      tATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTAATACCAGAAGCGACTGGAATACt                                                                   >  1:7436/1‑69 (MQ=255)
          tataAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGca                                                              >  1:545837/1‑70 (MQ=255)
          tataAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGca                                                              >  1:4080831/1‑70 (MQ=255)
           ataAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGcaca                                                            >  1:1998862/1‑71 (MQ=255)
                tgtgTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCAcatca                                                         >  1:2733851/1‑69 (MQ=255)
                   gtACAACGGGGATTGTTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATaaa                                                     >  1:3703292/1‑70 (MQ=255)
                   gtACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATaaa                                                     >  1:412276/1‑70 (MQ=255)
                    tACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATaaaa                                                    <  1:2127084/70‑1 (MQ=255)
                          ggggATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTatca                                             >  1:3310837/1‑71 (MQ=255)
                          ggggATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTatca                                             >  1:2953512/1‑71 (MQ=255)
                          ggggATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTatca                                             <  1:2920164/71‑1 (MQ=255)
                               ttattTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAAtt                                                 <  1:853292/62‑1 (MQ=255)
                               ttattTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAAtt                                                 <  1:645121/62‑1 (MQ=255)
                                 attTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTat                                               >  1:787284/1‑62 (MQ=255)
                                 attTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTat                                               >  1:3012190/1‑62 (MQ=255)
                                  tttttACATTGCATATTATTACCAGAAGCGACTGGaa                                                                       <  1:261891/37‑1 (MQ=255)
                                  tttttACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATgg                                     >  1:3474771/1‑71 (MQ=255)
                                    tttACATTGCATATTATTACCAGAAGCGACTGGAGTACTAAGCACATCATaa                                                      >  1:3557629/1‑52 (MQ=255)
                                    tttACATTGCATATTATTACCAGAAGCGACTGGAGTACTAAGCACATCATaa                                                      >  1:2384618/1‑52 (MQ=255)
                                    tttACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGta                                   >  1:3682822/1‑71 (MQ=255)
                                    tttACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGta                                   >  1:3400340/1‑71 (MQ=255)
                                       aCATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGtata                                 <  1:1130594/70‑1 (MQ=255)
                                           tGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTg                             >  1:960539/1‑70 (MQ=255)
                                            gaatatTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTAtcatca                                          >  1:2064276/3‑56 (MQ=255)
                                            gCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTAtcatca                                          >  1:1997283/1‑56 (MQ=255)
                                            gCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTAtcatca                                          >  1:2318437/1‑56 (MQ=255)
                                            gCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATACATGGTATATCTGaaa                          >  1:56067/1‑71 (MQ=255)
                                             cATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGaaa                          <  1:1928345/71‑1 (MQ=255)
                                              atatTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGaaa                          >  1:2356132/1‑70 (MQ=255)
                                              atatTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGaa                           >  1:1183985/1‑69 (MQ=255)
                                                attattACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTatcatc                                           <  1:1778906/51‑1 (MQ=255)
                                                 ttattaCCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTat                                               >  1:992190/1‑46 (MQ=255)
                                                  tattaCCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATcaca                                        >  1:677156/1‑52 (MQ=255)
                                                  tattaCCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAacacac                     >  1:1385989/1‑71 (MQ=255)
                                                  tattaCCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAacacac                     >  1:2642483/1‑71 (MQ=255)
                                                  tattaCCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAacacac                     >  1:3244285/1‑71 (MQ=255)
                                                  tattaCCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAacaca                      >  1:3857787/1‑70 (MQ=255)
                                                   attaCCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAacacac                     >  1:1606411/1‑70 (MQ=255)
                                                   attaCCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAacaca                      >  1:1748777/1‑69 (MQ=255)
                                                      aCCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGa                   >  1:2271919/1‑69 (MQ=255)
                                                      aCCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGa                   >  1:3541603/1‑69 (MQ=255)
                                                      aCCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGa                   >  1:4048192/1‑69 (MQ=255)
                                                      aCCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGa                   >  1:401572/1‑69 (MQ=255)
                                                      aCCAGAAGCGACTGGAATACTAAGAACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGa                   >  1:1550663/1‑69 (MQ=255)
                                                       ccAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAg                <  1:2544334/71‑1 (MQ=255)
                                                       ccAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAg                <  1:2090178/71‑1 (MQ=255)
                                                       ccAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAg                <  1:1378146/71‑1 (MQ=255)
                                                        cAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGa                   <  1:3806718/67‑1 (MQ=255)
                                                        cAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGa                   <  1:415985/67‑1 (MQ=255)
                                                         aGAAGCGACTGGAGTACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCa              <  1:3435548/71‑1 (MQ=255)
                                                         aGAAGCGACTGGAGTACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCa              <  1:3603101/71‑1 (MQ=255)
                                                         aGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCa              >  1:4097921/1‑71 (MQ=255)
                                                          gAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGa                            <  1:3877102/56‑1 (MQ=255)
                                                          gAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGa                            <  1:3343184/56‑1 (MQ=255)
                                                          gAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGct              >  1:1753123/1‑69 (MQ=255)
                                                          gAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCa              >  1:1096040/1‑70 (MQ=255)
                                                          gAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCa              >  1:701396/1‑70 (MQ=255)
                                                          gAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCa              >  1:2650140/1‑70 (MQ=255)
                                                          gAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCAt             <  1:3689813/71‑1 (MQ=255)
                                                          gAAGCGACTGGAATACTAAGAACATCATAAAAATTATCATCACATGGTATATCTGa                            <  1:1975077/56‑1 (MQ=255)
                                                           aaGCGACTGGAGTACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCa              <  1:1954542/69‑1 (MQ=255)
                                                           aaGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGa                            <  1:1235585/55‑1 (MQ=255)
                                                            aGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTatat                                >  1:2224499/1‑50 (MQ=255)
                                                            aGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTatat                                >  1:1394075/1‑50 (MQ=255)
                                                            aGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTatat                                >  1:2556686/1‑50 (MQ=255)
                                                            aGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGaaa                          >  1:2790714/1‑56 (MQ=255)
                                                            aGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGaaa                          >  1:3155539/1‑56 (MQ=255)
                                                            aGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATAc           >  1:3764308/1‑71 (MQ=255)
                                                             gCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGtat                                  >  1:3458593/1‑47 (MQ=255)
                                                                  tGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAACCACACGAAAGCATACCTaaa      >  1:2323171/1‑70 (MQ=255)
                                                                  tGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAg     >  1:1757028/1‑71 (MQ=255)
                                                                  tGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAg     >  1:3187121/1‑71 (MQ=255)
                                                                  tGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAg     >  1:891169/1‑71 (MQ=255)
                                                                  tGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAg     >  1:2869563/1‑71 (MQ=255)
                                                                    gAATACTAAGCACATCATAAAAATTATGTTCACATGTTATATCTGAAACACACGAAAGCATTCCTAAAg     <  1:1014412/69‑1 (MQ=255)
                                                                    gAATACTAAGCACATCATAAAAATTATCATCACATTGTATATCTGAAACACACGAAAGCATACCTAAAGt    <  1:2525824/70‑1 (MQ=255)
                                                                    gAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGt    <  1:1671510/70‑1 (MQ=255)
                                                                    gAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGt    <  1:3816326/70‑1 (MQ=255)
                                                                    gAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGt    <  1:234620/70‑1 (MQ=255)
                                                                    gAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGt    <  1:2475129/70‑1 (MQ=255)
                                                                    gAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGt    >  1:3108988/1‑70 (MQ=255)
                                                                    gAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGt    >  1:1122148/1‑70 (MQ=255)
                                                                    gAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGt    <  1:2681067/70‑1 (MQ=255)
                                                                     aaTACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTa  <  1:3169737/71‑1 (MQ=255)
                                                                     aaTACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTa  <  1:2947569/71‑1 (MQ=255)
                                                                     aaTACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTa  <  1:2895163/71‑1 (MQ=255)
                                                                      |                                                                     
GCCATCTATTTATAAATGTGTACAACGGGGATTATTTTTACATTGCATATTATTACCAGAAGCGACTGGAATACTAAGCACATCATAAAAATTATCATCACATGGTATATCTGAAACACACGAAAGCATACCTAAAGTTA  >  minE/2572488‑2572627

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