Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,609,051 T→C 53.9% S121S (TCA→TCG nikE ← nickel transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,609,0510TC53.9% 11.9 / 63.7 78S121S (TCA→TCGnikEnickel transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (19/17);  new base C (13/29);  total (32/46)
Fisher's exact test for biased strand distribution p-value = 6.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

ACGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGCGCGGATTCACGGCGCT  >  minE/2608983‑2609116
                                                                    |                                                                 
aCGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGt                                                                   <  1:2794233/69‑1 (MQ=255)
aCGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGt                                                                   <  1:2948544/69‑1 (MQ=255)
aCGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGt                                                                   <  1:3338458/69‑1 (MQ=255)
aCGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGt                                                                   <  1:1987982/69‑1 (MQ=255)
aCGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGt                                                                   <  1:1510551/69‑1 (MQ=255)
          cGAGATCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                                         <  1:54582/69‑1 (MQ=255)
          cGAGATCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                                         <  1:2891543/69‑1 (MQ=255)
          cGAGATCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                                         <  1:2519016/69‑1 (MQ=255)
          cGAGATCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                                         <  1:2195212/69‑1 (MQ=255)
          cGAGATCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                                         <  1:2109396/69‑1 (MQ=255)
          cGAGATCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                                         <  1:879539/69‑1 (MQ=255)
           gagaTCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                                         <  1:1126324/68‑1 (MQ=255)
           gagaTCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                                         <  1:1631808/68‑1 (MQ=255)
                   ccGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                                         <  1:1069421/60‑1 (MQ=255)
                       cTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGt                                                                   <  1:2402276/46‑1 (MQ=255)
                          cAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCa                                           >  1:1526053/1‑67 (MQ=255)
                          cAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCa                                           >  1:40801/1‑67 (MQ=255)
                          cAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCa                                           >  1:920148/1‑67 (MQ=255)
                          cAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCa                                           >  1:1681964/1‑67 (MQ=255)
                           aGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGt                                                                   <  1:2099985/42‑1 (MQ=255)
                                 tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAgg                                     <  1:1582782/66‑1 (MQ=255)
                                 tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAgg                                     <  1:1749744/66‑1 (MQ=255)
                                 tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAgg                                     <  1:1787739/66‑1 (MQ=255)
                                 tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAgg                                     <  1:2538344/66‑1 (MQ=255)
                                 tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAgg                                     <  1:3306171/66‑1 (MQ=255)
                                 tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAgg                                     <  1:16313/66‑1 (MQ=255)
                                 tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGAtc                                  >  1:2761785/1‑69 (MQ=255)
                                 tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGAtc                                  >  1:2635612/1‑69 (MQ=255)
                                 tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGAt                                   >  1:155804/1‑68 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGTCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:1817470/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:2631107/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:2692281/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:2791509/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:933585/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:3309871/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:1278998/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:464259/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:1405535/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:828141/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:708175/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:183580/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:579620/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:1507365/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:514779/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:910291/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:452370/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGAGCTCGCGCACGGTTTTg                   <  1:2530828/69‑1 (MQ=255)
                                                aGTTGTTCGGCTTTTTTCAGCGACAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:3070536/69‑1 (MQ=255)
                                                aGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGtttt                    >  1:356861/1‑68 (MQ=255)
                                                aGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   >  1:1880358/1‑69 (MQ=255)
                                                aGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   >  1:1639663/1‑69 (MQ=255)
                                                 gttgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg                   <  1:2904514/68‑1 (MQ=255)
                                                   tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGc                                        >  1:2941500/1‑45 (MQ=255)
                                                   tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                                       >  1:949112/1‑46 (MQ=255)
                                                   tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                                       >  1:585741/1‑46 (MQ=255)
                                                   tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                                       >  1:1021708/1‑46 (MQ=255)
                                                   tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                                       >  1:1045604/1‑46 (MQ=255)
                                                   tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                                       >  1:1927292/1‑46 (MQ=255)
                                                   tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                                       >  1:1952121/1‑46 (MQ=255)
                                                   tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                                       >  1:2617168/1‑46 (MQ=255)
                                                   tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                                       >  1:2573225/1‑46 (MQ=255)
                                                   tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                                       >  1:256523/1‑46 (MQ=255)
                                                   tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                                       >  1:2013588/1‑46 (MQ=255)
                                                   tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                                       >  1:237696/1‑46 (MQ=255)
                                                   tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                                       >  1:254982/1‑46 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCa                           >  1:2175930/1‑57 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCTCGGt                       <  1:1538559/61‑1 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt                       <  1:2398687/61‑1 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt                       <  1:1056239/61‑1 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt                       <  1:2870271/61‑1 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt                       <  1:1575904/61‑1 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt                       <  1:2938939/61‑1 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt                       <  1:1687247/61‑1 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt                       <  1:1722485/61‑1 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGcgcg               >  1:3243252/1‑69 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGcgcg               >  1:2899730/1‑69 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGcgcg               >  1:714885/1‑69 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGcgcg               >  1:3237035/1‑69 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGcgcg               >  1:2344185/1‑69 (MQ=255)
                                                    gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGcgcg               >  1:2348755/1‑69 (MQ=255)
                                                             ttttCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGCGCGGATTCACgg      >  1:1440534/1‑69 (MQ=255)
                                                             ttttCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGCGCGGATTCACg       >  1:278671/1‑68 (MQ=255)
                                                                  aGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGCGCGGATTCACGGCGCt  <  1:559256/68‑1 (MQ=255)
                                                                  aGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGCGCGGATTCACGGCGCt  <  1:2148042/68‑1 (MQ=255)
                                                                  aGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGCGCGGATTCACGGCGCt  <  1:1529797/68‑1 (MQ=255)
                                                                    |                                                                 
ACGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGCGCGGATTCACGGCGCT  >  minE/2608983‑2609116

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: