Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 822,224 G→T 11.4% 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.4% 277.9 / 13.1 87intergenic (‑200/+204)ldrC/chaAtoxic polypeptide, small/calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (41/36);  new base T (4/6);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

GGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTG  >  minE/822161‑822294
                                                               |                                                                      
gggATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCa                                                                 >  1:2846768/1‑71 (MQ=255)
gggATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCa                                                                 >  1:2202186/1‑71 (MQ=255)
 ggATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCa                                                                 <  1:1191939/70‑1 (MQ=255)
   aTAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCa                                                                 >  1:974432/1‑68 (MQ=255)
    tAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATctt                                                             <  1:2146569/71‑1 (MQ=255)
       gcccccGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATcttcac                                                          >  1:3273823/3‑71 (MQ=11)
       aGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATcttcac                                                          <  1:1666970/71‑1 (MQ=255)
        gCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATcttcact                                                         >  1:2919341/1‑71 (MQ=255)
        gCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATct                                                              <  1:3254999/66‑1 (MQ=255)
        cccccGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATcttcact                                                         >  1:2149067/2‑71 (MQ=255)
           ccGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGc                                                                       <  1:2942580/54‑1 (MQ=255)
           ccGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCa                                                                 <  1:1734175/60‑1 (MQ=11)
             gCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATc                                                    >  1:3324744/1‑71 (MQ=255)
             gCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATc                                                    >  1:1536436/1‑71 (MQ=255)
               aaCATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATccc                                                  >  1:2596323/1‑71 (MQ=255)
                 cATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCt                                                 >  1:2083464/1‑70 (MQ=255)
                  aTTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTg                                                >  1:2404976/1‑70 (MQ=255)
                  aTTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTg                                                >  1:3341871/1‑70 (MQ=255)
                  aTTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGa                                               <  1:2977040/71‑1 (MQ=255)
                   ttGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGa                                               >  1:2850771/1‑70 (MQ=255)
                   ttGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGa                                               >  1:3083187/1‑70 (MQ=255)
                   ttGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGa                                               >  1:210782/1‑70 (MQ=255)
                    tGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGaaa                                             >  1:1574881/1‑71 (MQ=255)
                     gcgcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAAt                                            >  1:1760770/1‑71 (MQ=255)
                      cgcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATcc                                                   <  1:3879605/63‑1 (MQ=255)
                      cgcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATcc                                                   <  1:2494901/63‑1 (MQ=255)
                       gcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATcttcactt                                                        >  1:3874831/1‑57 (MQ=255)
                       gcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGt                                          <  1:2157772/71‑1 (MQ=255)
                       gcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGt                                          >  1:1570957/1‑71 (MQ=255)
                       gcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGt                                          <  1:1820328/71‑1 (MQ=255)
                       gcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGt                                          <  1:724007/71‑1 (MQ=255)
                        ccTCGCCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGtt                                         >  1:134454/1‑71 (MQ=255)
                         cTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGc                                                                       <  1:2779871/40‑1 (MQ=255)
                          tCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATtttt                                        <  1:3547786/70‑1 (MQ=255)
                            aCCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTc                                                                  <  1:546351/42‑1 (MQ=255)
                            aCCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATcttc                                                            <  1:1004604/48‑1 (MQ=255)
                            aCCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                      >  1:3029493/1‑70 (MQ=255)
                            aCCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                      >  1:308203/1‑70 (MQ=255)
                            aCCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                      >  1:818537/1‑70 (MQ=255)
                            aCCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                      <  1:3863499/70‑1 (MQ=255)
                            aCCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                      <  1:1311797/70‑1 (MQ=255)
                              cGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATcttcac                                                          <  1:3536497/48‑1 (MQ=255)
                                gATAACGCCGGCTTGGTGTGGATACTACGTCGCaa                                                                     >  1:3817954/1‑35 (MQ=255)
                                 aTAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATct                                                              >  1:3427354/1‑41 (MQ=255)
                                  tAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATg                                  >  1:338974/1‑68 (MQ=255)
                                  tAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGaa                                <  1:2705853/70‑1 (MQ=255)
                                    aCGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATg                                  <  1:1584760/66‑1 (MQ=255)
                                    aCGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAag                               <  1:1304268/69‑1 (MQ=255)
                                         ggCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAagag                             <  1:2909769/66‑1 (MQ=40)
                                          gCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTTAAATGTTTGc                                      <  1:808797/56‑1 (MQ=255)
                                           cTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTc                      >  1:2463967/1‑71 (MQ=39)
                                             tGGTGTGGATACTACGTCGCAATTCCTCTTCACTTCATCCCTGAAATGTTTGCAATGaa                                >  1:3405652/1‑59 (MQ=255)
                                              ggTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACAtt                       >  1:121181/1‑67 (MQ=40)
                                                tgtgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAagag                             <  1:664239/59‑1 (MQ=255)
                                                  tgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGt                                          <  1:2370879/44‑1 (MQ=38)
                                                  tgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGt                                          <  1:122685/44‑1 (MQ=38)
                                                    gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCaa                                    >  1:3553029/1‑48 (MQ=38)
                                                    gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCtttt              >  1:741679/1‑70 (MQ=40)
                                                    gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc             <  1:677457/71‑1 (MQ=39)
                                                    gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc             <  1:596503/71‑1 (MQ=39)
                                                    gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc             <  1:2948501/71‑1 (MQ=39)
                                                      tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAagag                             <  1:3282921/53‑1 (MQ=255)
                                                      tactacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAag                               <  1:114633/51‑1 (MQ=255)
                                                       actacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa           >  1:3517161/1‑70 (MQ=40)
                                                        ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGaa                                >  1:1324150/1‑48 (MQ=255)
                                                        ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc             <  1:1844206/67‑1 (MQ=40)
                                                        ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc             <  1:406724/67‑1 (MQ=40)
                                                        ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa           >  1:1021518/1‑69 (MQ=40)
                                                        ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa           >  1:496374/1‑69 (MQ=40)
                                                        ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa           >  1:3558902/1‑69 (MQ=40)
                                                        ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa           >  1:3371847/1‑69 (MQ=40)
                                                         tacGTCGCAATTCATCTTCACTTCATCCCTGAAAt                                            <  1:1978852/35‑1 (MQ=40)
                                                         tacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                      <  1:3266395/41‑1 (MQ=38)
                                                         tacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAaga                              >  1:4060533/1‑49 (MQ=255)
                                                         tacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTa                           <  1:1121352/52‑1 (MQ=255)
                                                         tacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCtt                <  1:2674587/63‑1 (MQ=40)
                                                         tacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAg        >  1:362386/1‑71 (MQ=39)
                                                           cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAag                               >  1:1898399/1‑46 (MQ=255)
                                                           cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAg        <  1:1557856/69‑1 (MQ=40)
                                                           cGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCCTTTc             <  1:1594818/64‑1 (MQ=255)
                                                            gTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTAc                          >  1:1784858/1‑50 (MQ=255)
                                                            gTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTAc                          >  1:2047936/1‑50 (MQ=255)
                                                            gTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTAc                          >  1:695950/1‑50 (MQ=255)
                                                             tCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa           >  1:1576839/1‑64 (MQ=40)
                                                             tCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa           >  1:2428215/1‑64 (MQ=40)
                                                              cGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGc       <  1:1685831/67‑1 (MQ=39)
                                                              cGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGc       <  1:3551563/67‑1 (MQ=39)
                                                               gCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCt                 >  1:3474826/1‑56 (MQ=38)
                                                               gCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCt                 >  1:1882684/1‑56 (MQ=38)
                                                               gCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGtt   <  1:3427547/70‑1 (MQ=37)
                                                               gCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTg  >  1:2463625/1‑71 (MQ=35)
                                                               gCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTg  <  1:866019/71‑1 (MQ=35)
                                                               |                                                                      
GGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAGCTGTTG  >  minE/822161‑822294

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