Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 822,222 T→A 11.7% intergenic (‑198/+206) ldrC ← / ← chaA toxic polypeptide, small/calcium/sodium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE822,2220TA11.7% 205.2 / 17.8 77intergenic (‑198/+206)ldrC/chaAtoxic polypeptide, small/calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (40/28);  new base A (2/7);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 7.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

ATGCCACCAGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGT  >  minE/822152‑822292
                                                                      |                                                                      
aTGCCACCAGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGt                                                                        <  1:1786588/71‑1 (MQ=255)
     accaGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCaa                                                                   <  1:1319901/71‑1 (MQ=255)
     accaCGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCaa                                                                   <  1:960485/71‑1 (MQ=255)
      ccaGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAAt                                                                  >  1:1118681/1‑71 (MQ=255)
      ccaGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAAt                                                                  >  1:1890918/1‑71 (MQ=255)
       caGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGt                                                                        <  1:2079661/64‑1 (MQ=255)
       caGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGt                                                                        <  1:521496/64‑1 (MQ=255)
         gggATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCa                                                               >  1:3710635/1‑71 (MQ=255)
          ggATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCAt                                                              >  1:488825/1‑71 (MQ=255)
           gATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCg                                                                      <  1:147665/62‑1 (MQ=255)
               agcccccGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATcttca                                                         <  1:1171231/68‑1 (MQ=1)
                aGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATcttcac                                                        >  1:3693127/1‑71 (MQ=255)
                 gCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTAATcttcact                                                       >  1:797305/1‑71 (MQ=255)
                     cGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCAt                                                   <  1:1966961/71‑1 (MQ=255)
                      gCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCAGCTTCACTTCATc                                                  >  1:1616908/1‑71 (MQ=255)
                       cAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATcc                                                 <  1:2177149/71‑1 (MQ=255)
                       cAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATcc                                                 >  1:1440127/1‑71 (MQ=255)
                       cAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATcc                                                 >  1:3526780/1‑71 (MQ=255)
                       cAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATcc                                                 >  1:3464517/1‑71 (MQ=255)
                       cAACATTGAGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATcc                                                 >  1:2815525/1‑71 (MQ=255)
                          cATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATctt                                                           <  1:3083292/58‑1 (MQ=255)
                          cATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCt                                               >  1:785521/1‑70 (MQ=255)
                               cgcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAAt                                          <  1:1798398/70‑1 (MQ=255)
                               cgcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATg                                         <  1:3244186/71‑1 (MQ=255)
                               cgcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATg                                         <  1:1434139/71‑1 (MQ=255)
                               cgcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATg                                         <  1:1359373/71‑1 (MQ=255)
                                 ccTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGtt                                       >  1:426196/1‑71 (MQ=255)
                                  cTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATCttt                                      <  1:3455081/71‑1 (MQ=255)
                                   tCACCGGATAACGCCGGCTTGGTGTGGATACTACGt                                                                        <  1:2084305/36‑1 (MQ=255)
                                   tCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTc                                                                <  1:451728/44‑1 (MQ=255)
                                   tCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGtttt                                     >  1:1613180/1‑70 (MQ=255)
                                    cACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGTTTg                                     >  1:3678782/1‑70 (MQ=255)
                                     aCCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATc                                                             <  1:896309/45‑1 (MQ=255)
                                     aCCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                    >  1:1943507/1‑70 (MQ=255)
                                     aCCGGATAACGCCGGCTTGGTGTGGATACTAAGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                    >  1:3512277/1‑70 (MQ=255)
                                      ccGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCa                                                               >  1:1194245/1‑42 (MQ=255)
                                      ccGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGaaa                                           >  1:2671251/1‑62 (MQ=255)
                                      ccGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                    <  1:1976913/69‑1 (MQ=255)
                                       cGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATg                                         <  1:3217359/63‑1 (MQ=255)
                                       cGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCaa                                  >  1:3658768/1‑70 (MQ=255)
                                        ggATAACGCCGGCTTGGTGTGGATACTACGTCGCAAtt                                                                 <  1:1593523/38‑1 (MQ=255)
                                         gATAACGCCGGCTTGGTGTGGATACTACGACTCAAt                                                                  <  1:551162/36‑1 (MQ=14)
                                          aTAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATg                                >  1:2185348/1‑69 (MQ=255)
                                          aTAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATg                                >  1:834305/1‑69 (MQ=255)
                                           tAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATg                                >  1:1250754/1‑68 (MQ=255)
                                             aCGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAagag                           >  1:2056928/1‑71 (MQ=255)
                                              cGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCATGAAATGTTTGCAATGAagag                           >  1:1787643/1‑70 (MQ=255)
                                              cGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAagag                           >  1:3588394/1‑70 (MQ=255)
                                                ccGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGtt                                       >  1:151007/1‑56 (MQ=35)
                                                 cGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACa                       >  1:1808010/1‑71 (MQ=39)
                                                  ggCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACAt                      <  1:168253/71‑1 (MQ=39)
                                                    cTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCa                                   <  1:1692959/56‑1 (MQ=35)
                                                    cTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTc                    <  1:1373289/71‑1 (MQ=39)
                                                     ttGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTcc                   >  1:819932/1‑71 (MQ=39)
                                                       ggTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTcc                   <  1:458069/69‑1 (MQ=40)
                                                       ggTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCgg                 >  1:2280375/1‑71 (MQ=39)
                                                        gtgtGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGc                <  1:1834663/71‑1 (MQ=39)
                                                           tgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGaa                                            <  1:2664049/40‑1 (MQ=40)
                                                           tgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGaa                                            <  1:2365845/40‑1 (MQ=40)
                                                           tgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGt                          <  1:1337933/58‑1 (MQ=255)
                                                            ggATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                    <  1:2626367/47‑1 (MQ=38)
                                                             taTACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCtttt            <  1:1831327/69‑1 (MQ=255)
                                                             gATACTACGTCGCAATTCATCTTCACTTCATCCCt                                               <  1:3536035/35‑1 (MQ=40)
                                                             gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc           <  1:2201998/71‑1 (MQ=39)
                                                              aTACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCa          <  1:3503504/71‑1 (MQ=39)
                                                               tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa         <  1:582289/71‑1 (MQ=39)
                                                                 ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa       <  1:3813511/71‑1 (MQ=39)
                                                                  tacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                    >  1:1901168/1‑41 (MQ=38)
                                                                  tacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                    >  1:644106/1‑41 (MQ=38)
                                                                  tacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                    >  1:3789318/1‑41 (MQ=38)
                                                                  tacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                    >  1:2587123/1‑41 (MQ=38)
                                                                   acGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt             >  1:3703888/1‑63 (MQ=40)
                                                                   acGTCGCAATTCATCTTAACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCttt             >  1:153315/1‑63 (MQ=255)
                                                                    cGTCGCAATTCATCTTCACTTCATCCGTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCt    >  1:3787728/1‑71 (MQ=255)
                                                                    cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCgg                 >  1:2636052/1‑58 (MQ=255)
                                                                    cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCgg                 >  1:366524/1‑58 (MQ=255)
                                                                    cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa       >  1:2580893/1‑68 (MQ=40)
                                                                    cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACa       >  1:184691/1‑68 (MQ=40)
                                                                    cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAg      >  1:1772258/1‑69 (MQ=40)
                                                                    cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCt    <  1:3704765/71‑1 (MQ=39)
                                                                    cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCt    <  1:3229484/71‑1 (MQ=39)
                                                                      tCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTg   <  1:3484285/70‑1 (MQ=40)
                                                                      tCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGt  <  1:1046257/71‑1 (MQ=39)
                                                                      |                                                                      
ATGCCACCAGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGT  >  minE/822152‑822292

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: