Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 822,224 G→T 11.5% intergenic (‑200/+204) 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,2240GT11.5% 243.4 / 10.4 78intergenic (‑200/+204)ldrC/chaAtoxic polypeptide, small/calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (40/29);  new base T (2/7);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 7.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

ACCAGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGT  >  minE/822157‑822292
                                                                   |                                                                    
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)
    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:3526780/1‑71 (MQ=255)
                  cAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATcc                                                 >  1:3464517/1‑71 (MQ=255)
                  cAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATcc                                                 <  1:2177149/71‑1 (MQ=255)
                  cAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATcc                                                 >  1:1440127/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:1359373/71‑1 (MQ=255)
                          cgcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATg                                         <  1:1434139/71‑1 (MQ=255)
                            ccTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGtt                                       >  1:426196/1‑71 (MQ=255)
                             cTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATCttt                                      <  1:3455081/71‑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:2587123/1‑41 (MQ=38)
                                                             tacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                    >  1:644106/1‑41 (MQ=38)
                                                             tacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                    >  1:1901168/1‑41 (MQ=38)
                                                             tacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                    >  1:3789318/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)
                                                                   |                                                                    
ACCAGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGT  >  minE/822157‑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: