Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,469,879 A→G 13.1% K287K (AAA→AAG atpD → F1 sector of membrane‑bound ATP synthase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,469,8790AG13.1% 208.5 / 20.6 76K287K (AAA→AAGatpDF1 sector of membrane‑bound ATP synthase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (44/22);  new base G (5/5);  total (49/27)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TAGGTTATCAGCCGACCCTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTTGACTGACCCGTCTCCGGC  >  minE/2469815‑2469947
                                                                |                                                                    
tAGGTTATCAGCCGACCCTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGt                                                                >  1:227519/1‑71 (MQ=255)
  ggTTATCAGCCGACCCTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCaaa                                                                      <  1:3469018/63‑1 (MQ=255)
   gTTATCAGCCGACCCTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCt                                                             <  1:2340304/71‑1 (MQ=255)
     tATCAGCCGACCCTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTa                                                            >  1:3183917/1‑70 (MQ=255)
      aTCAGCCGACCCTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTAt                                                           <  1:1744570/70‑1 (MQ=255)
      aTCAGCCGACCCTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTAt                                                           >  1:2043747/1‑70 (MQ=255)
       tCAGCCGACCCTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCa                                                         >  1:3226352/1‑71 (MQ=255)
        cAGCCGACCCTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCAc                                                        >  1:3237604/1‑71 (MQ=255)
            cGACCCTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTcc                                                    >  1:3212824/1‑71 (MQ=255)
                 cTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTa                                                            <  1:2875385/58‑1 (MQ=255)
                 cTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTAt                                                           >  1:1567215/1‑59 (MQ=255)
                 cTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTa                                                 >  1:1963583/1‑69 (MQ=255)
                   ggCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATc                                                          >  1:2211391/1‑58 (MQ=255)
                   ggCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATc                                                          >  1:415046/1‑58 (MQ=255)
                   ggCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATc                                                          >  1:2466996/1‑58 (MQ=255)
                      ggAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAg                                          <  1:1257744/71‑1 (MQ=255)
                       gAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAg                                          <  1:1887403/70‑1 (MQ=255)
                        aaGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGta                                        <  1:861823/71‑1 (MQ=255)
                         agagATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGta                                        <  1:1806405/70‑1 (MQ=255)
                          gagaTGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCt                                                             >  1:1888344/1‑48 (MQ=255)
                          gagaTGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCt                                                             >  1:1956310/1‑48 (MQ=255)
                          gagaTGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCt                                                             >  1:829358/1‑48 (MQ=255)
                             aTGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACg                                    >  1:3367624/1‑70 (MQ=255)
                             aTGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACg                                    >  1:2949495/1‑70 (MQ=255)
                             aTGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACg                                    >  1:277513/1‑70 (MQ=255)
                                 gCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACa                                               >  1:98740/1‑55 (MQ=255)
                                 gCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACa                                               >  1:2833313/1‑55 (MQ=255)
                                 gCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACa                                               >  1:1205519/1‑55 (MQ=255)
                                 gCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACa                                               >  1:2234281/1‑55 (MQ=255)
                                 gCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACa                                               >  1:622869/1‑55 (MQ=255)
                                 gCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACa                                               >  1:3231985/1‑55 (MQ=255)
                                  cGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCt                               >  1:1155231/1‑70 (MQ=255)
                                  cGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTg                              >  1:2782294/1‑71 (MQ=255)
                                    ttCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCg                            <  1:3397186/71‑1 (MQ=255)
                                    ttCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCg                            <  1:3401195/71‑1 (MQ=255)
                                      cTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAg                                          <  1:2051383/55‑1 (MQ=255)
                                      cTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAg                                          <  1:1734919/55‑1 (MQ=255)
                                      cTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAg                                          <  1:1209989/55‑1 (MQ=255)
                                      cTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGa                          >  1:2044478/1‑71 (MQ=255)
                                      cTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGa                          >  1:1775969/1‑71 (MQ=255)
                                        gCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGta                                        >  1:3258571/1‑55 (MQ=255)
                                         cAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGAt                         <  1:602824/69‑1 (MQ=255)
                                            gAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGAtgactt                    <  1:3045456/71‑1 (MQ=255)
                                            gAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGAtgactt                    <  1:30089/71‑1 (MQ=255)
                                            gAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGAtgactt                    <  1:1441752/71‑1 (MQ=255)
                                                gTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGAtgacttgact                <  1:1575890/71‑1 (MQ=255)
                                                gTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGAtgacttga                  >  1:830049/1‑69 (MQ=255)
                                                 tATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGAtgacttgact                <  1:2384686/70‑1 (MQ=255)
                                                 tATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGAtgacttgact                <  1:2947086/70‑1 (MQ=255)
                                                  aTCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTtgactga              >  1:1691340/1‑71 (MQ=255)
                                                   tCACCTCCACCAAGAAAGGCTCGATCACCTCGATCCAGGCCGTATACGttc                                 <  1:2837128/51‑3 (MQ=255)
                                                   tCACCTCCACCAAGAAAGGCTCGATCACCTCGATCCAGGCCGTATACGttc                                 <  1:2320244/51‑3 (MQ=255)
                                                   tCACCTCCACCAAGAAAGGCTCGATCACCTCGATCCAGGCCGTATACGttc                                 <  1:1522286/51‑3 (MQ=255)
                                                   tCACCTCCACCAAGAAAGGCTCGATCACCTCGATCCAGGCCGTATACGttc                                 <  1:2273389/51‑3 (MQ=255)
                                                   tCACCTCCACCAAGAAAGGCTCGATCACCTCGATCCAGGCCGTATACGttc                                 <  1:357078/51‑3 (MQ=255)
                                                   tCACCTCCACCAAAACTGGTTCTATCACCTCCGTa                                                 <  1:2545540/35‑1 (MQ=255)
                                                   tCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTtgactgac             >  1:101238/1‑71 (MQ=255)
                                                    cACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGtat                                       >  1:1618831/1‑44 (MQ=255)
                                                    cACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGtat                                       >  1:3376782/1‑44 (MQ=255)
                                                    cACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGta                                        >  1:1589875/1‑43 (MQ=255)
                                                    cACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTTGACTGAcc            >  1:904085/1‑71 (MQ=255)
                                                    cACCTCAACCAAAACTGGTTCTATCACCTCCGTACAGGCAGtat                                       >  1:1853/1‑44 (MQ=255)
                                                     aCCTCCACCAAGAAAGGCTCGATCACCTCGATCCAGGCCGTATACGTTccc                               >  1:1656986/1‑50 (MQ=255)
                                                     aCCTCCACCAAGAAAGGCTCGATCACCTCGATCCAGGCCGTATACGTTccc                               >  1:1456628/1‑50 (MQ=255)
                                                     aCCTCCACCAAGAAAGGCTCGATCACCTCGATCCAGGCCGTATACGTTccc                               >  1:807223/1‑50 (MQ=255)
                                                     aCCTCCACCAAGAAAGGCTCGATCACCTCGATCCAGGCCGTATACGTTccc                               >  1:1711370/1‑50 (MQ=255)
                                                     aCCTCCACCAAGAAAGGCTCGATCACCTCGAACCAGGCCGTATACGttc                                 >  1:230767/1‑47 (MQ=255)
                                                     aCCTCCACCAAAACTGGTTCTATCACCTCCGTACa                                               >  1:1177211/1‑35 (MQ=255)
                                                     aCCTCCACCAAAACTGGTTCTATCACCTCCGTACa                                               >  1:902396/1‑35 (MQ=255)
                                                      ccTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTAc                                 >  1:847435/1‑48 (MQ=255)
                                                      ccTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTTGACTGACCCg          <  1:1487376/71‑1 (MQ=255)
                                                       cTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTTGACTGACCCGt         <  1:939294/71‑1 (MQ=255)
                                                          caccaAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTTGACTGACCCGtctc      >  1:521784/1‑71 (MQ=255)
                                                          caccaAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTTGACTGACCCGtctc      >  1:2021175/1‑71 (MQ=255)
                                                           accaAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTTGACTGACCCGTCTcc     >  1:1484170/1‑71 (MQ=255)
                                                           accaAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTTGACTGACCCGTCTcc     >  1:854831/1‑71 (MQ=255)
                                                           accaAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTTGACTGACCCGTCTcc     >  1:1149309/1‑71 (MQ=255)
                                                              aaaaCTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTTGACTGACCCGTCTCCgg   >  1:2002623/1‑70 (MQ=255)
                                                              aaaaCTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTTGACTGACCCGTCTCCGGc  <  1:493031/71‑1 (MQ=255)
                                                                aaCTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGAtgacttga                  >  1:56397/1‑53 (MQ=255)
                                                                |                                                                    
TAGGTTATCAGCCGACCCTGGCGGAAGAGATGGGCGTTCTGCAGGAACGTATCACCTCCACCAAAACTGGTTCTATCACCTCCGTACAGGCAGTATACGTACCTGCGGATGACTTGACTGACCCGTCTCCGGC  >  minE/2469815‑2469947

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: