Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,609,041 A→C 51.3% S125A (TCC→GCC)  nikE ← nickel transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,609,0410AC51.3% 12.2 / 89.0 82S125A (TCC→GCC) nikEnickel transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (19/21);  new base C (13/29);  total (32/50)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

AACGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGCGCG  >  minE/2608982‑2609103
                                                           |                                                              
aaCGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCa                                                        >  1:903532/1‑68 (MQ=255)
aaCGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCa                                                        >  1:364695/1‑68 (MQ=255)
 aCGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGt                                                      <  1:2948544/69‑1 (MQ=255)
 aCGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGt                                                      <  1:1987982/69‑1 (MQ=255)
 aCGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGt                                                      <  1:1510551/69‑1 (MQ=255)
 aCGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGt                                                      <  1:3338458/69‑1 (MQ=255)
 aCGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGt                                                      <  1:2794233/69‑1 (MQ=255)
           cGAGATCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                            <  1:2519016/69‑1 (MQ=255)
           cGAGATCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                            <  1:54582/69‑1 (MQ=255)
           cGAGATCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                            <  1:2195212/69‑1 (MQ=255)
           cGAGATCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                            <  1:2891543/69‑1 (MQ=255)
           cGAGATCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                            <  1:879539/69‑1 (MQ=255)
           cGAGATCGACCGCCTTCAGCATTTCACTGGCGCGCGCCAGTTGTTCGGCTTTTTTCAGCGATAGCAGGt                                            <  1:2109396/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:920148/1‑67 (MQ=255)
                           cAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCa                              >  1:40801/1‑67 (MQ=255)
                           cAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCa                              >  1:1681964/1‑67 (MQ=255)
                            aGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGt                                                      <  1:2099985/42‑1 (MQ=255)
                                  tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAgg                        <  1:2538344/66‑1 (MQ=255)
                                  tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAgg                        <  1:3306171/66‑1 (MQ=255)
                                  tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAgg                        <  1:1582782/66‑1 (MQ=255)
                                  tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAgg                        <  1:1787739/66‑1 (MQ=255)
                                  tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAgg                        <  1:16313/66‑1 (MQ=255)
                                  tCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAgg                        <  1:1749744/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:2791509/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:183580/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:2692281/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:1278998/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:3309871/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:2631107/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:1405535/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:933585/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:1507365/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:910291/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:452370/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:464259/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:514779/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:579620/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:708175/69‑1 (MQ=255)
                                                 aGTTGTTCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      <  1:828141/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:1639663/1‑69 (MQ=255)
                                                 aGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTg      >  1:1880358/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:254982/1‑46 (MQ=255)
                                                    tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                          >  1:1045604/1‑46 (MQ=255)
                                                    tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                          >  1:1927292/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:1021708/1‑46 (MQ=255)
                                                    tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                          >  1:1952121/1‑46 (MQ=255)
                                                    tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                          >  1:237696/1‑46 (MQ=255)
                                                    tgttCGGCTTTTTTCAGCGATAGCAGGTGGCGCATCGGTTCACGCa                          >  1:2013588/1‑46 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCa              >  1:2175930/1‑57 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCTCGGt          <  1:1538559/61‑1 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt          <  1:1722485/61‑1 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt          <  1:1687247/61‑1 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt          <  1:2938939/61‑1 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt          <  1:1575904/61‑1 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt          <  1:2398687/61‑1 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt          <  1:2870271/61‑1 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGt          <  1:1056239/61‑1 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGcgcg  >  1:2344185/1‑69 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGcgcg  >  1:714885/1‑69 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGcgcg  >  1:2348755/1‑69 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGcgcg  >  1:2899730/1‑69 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGcgcg  >  1:3243252/1‑69 (MQ=255)
                                                     gttCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGcgcg  >  1:3237035/1‑69 (MQ=255)
                                                           |                                                              
AACGCTGTCATCGAGATCGACCGCCTTCAGCATTTCGCTGGCGCGCGCCAGTTGTTCGGATTTTTTCAGTGATAGCAGGTGGCGCATCGGTTCACGCAGGATCTCGCGCACGGTTTTGCGCG  >  minE/2608982‑2609103

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: