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% 229.3 / 19.0 85intergenic (‑198/+206)ldrC/chaAtoxic polypeptide, small/calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (41/34);  new base A (4/6);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCACCAGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAG  >  minE/822153‑822288
                                                                     |                                                                  
tGCCACCAGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTc                                                                   >  1:2013062/1‑71 (MQ=255)
tGCCACCAGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTc                                                                   >  1:1787750/1‑71 (MQ=255)
        gggATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCa                                                           >  1:2202186/1‑71 (MQ=255)
        gggATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCa                                                           >  1:2846768/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:3341871/1‑70 (MQ=255)
                          aTTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTg                                          >  1:2404976/1‑70 (MQ=255)
                          aTTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGa                                         <  1:2977040/71‑1 (MQ=255)
                           ttGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGa                                         >  1:3083187/1‑70 (MQ=255)
                           ttGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGa                                         >  1:2850771/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:2494901/63‑1 (MQ=255)
                              cgcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATcc                                             <  1:3879605/63‑1 (MQ=255)
                               gcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATcttcactt                                                  >  1:3874831/1‑57 (MQ=255)
                               gcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGt                                    >  1:1570957/1‑71 (MQ=255)
                               gcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGt                                    <  1:2157772/71‑1 (MQ=255)
                               gcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGt                                    <  1:724007/71‑1 (MQ=255)
                               gcCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGt                                    <  1:1820328/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:818537/1‑70 (MQ=255)
                                    aCCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                >  1:308203/1‑70 (MQ=255)
                                    aCCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                >  1:3029493/1‑70 (MQ=255)
                                    aCCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                <  1:1311797/70‑1 (MQ=255)
                                    aCCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                                <  1:3863499/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:122685/44‑1 (MQ=38)
                                                          tgGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGt                                    <  1:2370879/44‑1 (MQ=38)
                                                            gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCaa                              >  1:3553029/1‑48 (MQ=38)
                                                            gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCtttt        >  1:741679/1‑70 (MQ=40)
                                                            gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc       <  1:2948501/71‑1 (MQ=39)
                                                            gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc       <  1:596503/71‑1 (MQ=39)
                                                            gATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc       <  1:677457/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:406724/67‑1 (MQ=40)
                                                                ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTc       <  1:1844206/67‑1 (MQ=40)
                                                                ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa     >  1:496374/1‑69 (MQ=40)
                                                                ctacGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa     >  1:1021518/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:695950/1‑50 (MQ=255)
                                                                    gTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTAc                    >  1:2047936/1‑50 (MQ=255)
                                                                     tCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa     >  1:1576839/1‑64 (MQ=40)
                                                                     tCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCaa     >  1:2428215/1‑64 (MQ=40)
                                                                     |                                                                  
TGCCACCAGGGATAAAGCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATGAAGAGTACATTCCGGCTTTTCAACAG  >  minE/822153‑822288

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: