Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,620,159 C→T 7.8% A194T (GCA→ACA)  zntA ← zinc, cobalt and lead efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,620,1590CT 7.8% 283.0 / 10.4 90A194T (GCA→ACA) zntAzinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base C (33/50);  new base T (5/2);  total (38/52)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

AGTCGTTCACCAATCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAACCGTAATGCC  >  minE/2620092‑2620229
                                                                   |                                                                      
aGTCGTTCACCAATCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGc                                                                        <  1:2949311/68‑1 (MQ=255)
aGTCGTTCACCAATCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCa                                                                       >  1:3212510/1‑69 (MQ=255)
aGTCGTTCACCAATCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCa                                                                       >  1:1948920/1‑69 (MQ=255)
  tCGTTCACCAATTAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCaa                                                                   <  1:2393611/71‑1 (MQ=255)
  tCGTTCACCAATCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCAc                                                                      >  1:2967426/1‑68 (MQ=255)
  tCGTTCACCAATCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCAc                                                                      >  1:1638634/1‑68 (MQ=255)
    gTTCACCAATCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGTACCaa                                                                   >  1:3316744/1‑69 (MQ=255)
    gTTCACCAATCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCaa                                                                   >  1:2369017/1‑69 (MQ=255)
    gTTCACAAATCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCaa                                                                   >  1:201812/1‑69 (MQ=255)
       cACCAATCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCg                                                              >  1:2583330/1‑71 (MQ=255)
        aCCAATCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTCCGCCAATAAACAGTGCACCAATAGCg                                                              <  1:2886677/70‑1 (MQ=255)
           aaTCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGTACCAATAGCGGCTa                                                          >  1:1023268/1‑71 (MQ=255)
             tCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCAcc                                                                     <  1:1734129/58‑1 (MQ=255)
             tCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACg                                                        <  1:177487/71‑1 (MQ=255)
              cAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCa                                                                       >  1:3052940/1‑55 (MQ=255)
              cAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTa                                                          >  1:37869/1‑68 (MQ=255)
               aaaaaCAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGc                                                                        <  1:3534269/53‑1 (MQ=255)
               aaaaaCAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGc                                                       <  1:2399430/70‑1 (MQ=255)
               aaaaaCAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGc                                                       >  1:3307465/1‑70 (MQ=255)
               aaaaaCAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGc                                                       <  1:1828731/70‑1 (MQ=255)
               aaaaaCAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCt                                                      <  1:280401/71‑1 (MQ=255)
               aaaaaCAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCt                                                      <  1:2897446/71‑1 (MQ=255)
                  aaCAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCAt                                                   >  1:2181295/1‑71 (MQ=255)
                   aCAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCAtt                                                  <  1:1649864/71‑1 (MQ=255)
                      gcagcaACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAg                                               <  1:3717744/71‑1 (MQ=255)
                       cagcaACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAgg                                              >  1:107298/1‑71 (MQ=255)
                         gcaACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGc                                                               <  1:3845221/52‑1 (MQ=255)
                          caACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCAcc                                                                     <  1:3733244/45‑1 (MQ=255)
                           aaCACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCTGCTACGCTCATTAAGGTTTc                                          <  1:101935/71‑1 (MQ=255)
                           aaCACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTc                                          <  1:3754322/71‑1 (MQ=255)
                           aaCACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTc                                          <  1:975994/71‑1 (MQ=255)
                           aaCACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTc                                          <  1:2603142/71‑1 (MQ=255)
                             cacCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATa                                                                 <  1:1364081/46‑1 (MQ=255)
                             cacCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCaa                                        <  1:2771954/71‑1 (MQ=255)
                               ccATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCaa                                                                   >  1:772096/1‑42 (MQ=255)
                               ccATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCAt                                                   <  1:2063498/58‑1 (MQ=255)
                                  tCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCaa                                                                   >  1:3253434/1‑39 (MQ=255)
                                   cGCAGCTTCAGCCGTTGCGCCAATAAACAGTGTACCAATAGCGGCTACGCTCATTAAGGTTTc                                          >  1:919201/1‑63 (MQ=255)
                                    gCAGCTTCAGCCGTTGCGCCAATAAACAGTGTACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCg                                   >  1:3969774/1‑69 (MQ=255)
                                    gCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGtt                                            <  1:2113965/60‑1 (MQ=255)
                                    gCAGCTTCAGCCGTTGCCCCAATAAACAGTGTACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCg                                   >  1:2460210/1‑69 (MQ=255)
                                      aGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGTCTACGCTCATTAAGGTTTCAATGGCGAAGt                               <  1:3803608/71‑1 (MQ=255)
                                      aGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGc                                    >  1:179590/1‑66 (MQ=255)
                                      aGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGt                               <  1:1308666/71‑1 (MQ=255)
                                      aGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGt                               <  1:2241279/71‑1 (MQ=255)
                                       gCTTCAGCCGTTGCGCCAATAAACAGTGTACCAATAGCGGCTACGCTCATTAAGGTTTCa                                         <  1:2838741/60‑1 (MQ=255)
                                        cTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGt                               <  1:3052074/69‑1 (MQ=255)
                                         ttCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGt                                             <  1:827236/54‑1 (MQ=255)
                                          tCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCaa                                        >  1:2650426/1‑58 (MQ=255)
                                          tCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCa                                         >  1:93150/1‑57 (MQ=255)
                                             gCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGcc                        <  1:1874567/71‑1 (MQ=255)
                                               cGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCaa                                        >  1:1513472/1‑53 (MQ=255)
                                                gTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGa                     <  1:1957528/71‑1 (MQ=255)
                                                 ttGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTc                                          <  1:3855958/49‑1 (MQ=255)
                                                 ttGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTc                                          <  1:3848412/49‑1 (MQ=255)
                                                 ttGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGAt                    <  1:784251/71‑1 (MQ=255)
                                                      ccAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTc                                          <  1:2923278/44‑1 (MQ=255)
                                                      ccAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGAt               <  1:1336855/71‑1 (MQ=255)
                                                       cAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCaa                                        >  1:458448/1‑45 (MQ=255)
                                                       cAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCaa                                        >  1:334615/1‑45 (MQ=255)
                                                       cAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCaa                                        >  1:3135896/1‑45 (MQ=255)
                                                       cAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCaa                                        >  1:979722/1‑45 (MQ=255)
                                                       cAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGAt               >  1:1743788/1‑70 (MQ=255)
                                                       cAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGAt               >  1:1049958/1‑70 (MQ=255)
                                                       cAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGAt               >  1:1890544/1‑70 (MQ=255)
                                                        aaTAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCg                       <  1:515382/61‑1 (MQ=255)
                                                        aaTAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATc              >  1:2792395/1‑70 (MQ=255)
                                                         aTAAACAGTGTACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTa                              <  1:1661552/53‑1 (MQ=255)
                                                         aTAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCaa            >  1:3647462/1‑71 (MQ=255)
                                                          tAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGt                                             <  1:2699195/37‑1 (MQ=255)
                                                          tAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCg                                   <  1:3007761/47‑1 (MQ=255)
                                                          tAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGAt               >  1:2946351/1‑67 (MQ=255)
                                                           aaaCAGCGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGa                >  1:830397/1‑65 (MQ=255)
                                                            aaCAGTGCACCAATAGCGGCTACGCTCATTAAGGtt                                            <  1:2930400/36‑1 (MQ=255)
                                                            aaCAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTc                                          <  1:2837455/38‑1 (MQ=255)
                                                            aaCAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGaa                                 <  1:13795/47‑1 (MQ=255)
                                                            aaCAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAAc           <  1:762568/69‑1 (MQ=255)
                                                            aaCAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAACCg         <  1:3260647/71‑1 (MQ=255)
                                                            aaCAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAACCg         >  1:1542861/1‑71 (MQ=255)
                                                             aCAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGaa                                 <  1:428274/46‑1 (MQ=255)
                                                             aCAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAACCGt        <  1:2192444/71‑1 (MQ=255)
                                                             aCAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAACCGt        <  1:1673181/71‑1 (MQ=255)
                                                             aCAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAACCGt        <  1:1217184/71‑1 (MQ=255)
                                                             aCAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAACCGt        <  1:2673820/71‑1 (MQ=255)
                                                              cAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATgg                                     >  1:674052/1‑41 (MQ=255)
                                                              cAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCg                                   >  1:2000989/1‑43 (MQ=255)
                                                               aGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTg                 <  1:4099098/60‑1 (MQ=255)
                                                               aGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAACCGTaa      <  1:3481106/71‑1 (MQ=255)
                                                                gTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTg                 <  1:886282/59‑1 (MQ=255)
                                                                gTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAACCGTAAt     <  1:3485398/71‑1 (MQ=255)
                                                                  gCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAACCGt        <  1:3326795/66‑1 (MQ=255)
                                                                  gCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAACCGTAATGc   >  1:2247398/1‑71 (MQ=255)
                                                                   cACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGaa                                 >  1:1734181/1‑40 (MQ=255)
                                                                   cACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAACCGTAATGcc  >  1:2018267/1‑71 (MQ=255)
                                                                   |                                                                      
AGTCGTTCACCAATCAAAAACAGCAGCAACACCATCGCAGCTTCAGCCGTTGCGCCAATAAACAGTGCACCAATAGCGGCTACGCTCATTAAGGTTTCAATGGCGAAGTAGCTGCCGGATTTGATCAACCGTAATGCC  >  minE/2620092‑2620229

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