Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,609,188 G→A 43.1% L76L (CTG→TTG)  nikE ← nickel transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,609,1880GA43.1% 11.3 / 29.3 44L76L (CTG→TTG) nikEnickel transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (17/8);  new base A (13/6);  total (30/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAAATACCATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTCGCCACGCCAGCTAATATTCCCCTGCGCGGGCGATTCTAAACCCACCAGCAACCGC  >  minE/2609126‑2609249
                                                              |                                                             
cTGAAATACCATTTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTc                                                         >  1:645122/1‑69 (MQ=255)
cTGAAATACCATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTc                                                         >  1:419966/1‑69 (MQ=255)
cTGAAATACCATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTc                                                         >  1:349043/1‑69 (MQ=255)
cTGAAATACCATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTc                                                         >  1:3419831/1‑69 (MQ=255)
cTGAAATACCATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTc                                                         >  1:3167416/1‑69 (MQ=255)
cTGAAATACCATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTc                                                         >  1:1777011/1‑69 (MQ=255)
cTGAAATACCATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTc                                                         >  1:2074092/1‑69 (MQ=255)
     aTACCATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTcgccac                                                    >  1:2028208/1‑69 (MQ=255)
       aCCATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTcgc                                                       <  1:571817/64‑1 (MQ=255)
       aCCATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTcgc                                                       <  1:619794/64‑1 (MQ=255)
       aCCATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTcgc                                                       <  1:3004084/64‑1 (MQ=255)
        ccATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTcgccacgcc                                                 >  1:2408594/1‑69 (MQ=255)
        ccATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTcgccacgcc                                                 <  1:349266/69‑1 (MQ=255)
        ccATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTcgccacgcc                                                 <  1:695482/69‑1 (MQ=255)
        ccATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTcgccacgcc                                                 >  1:2514322/1‑69 (MQ=255)
        ccATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTcgccacgcc                                                 >  1:250724/1‑69 (MQ=255)
        ccATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTcgccacgcc                                                 >  1:2305269/1‑69 (MQ=255)
            cTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTcgc                                                       <  1:323233/59‑1 (MQ=255)
             tGAATATCGCGGCGGAACGCTTTTCGCTGGGCGCGATTGAGTTTCGCCAGCGGTTCGCCACGCCAGCt                                             >  1:1468528/1‑68 (MQ=255)
             tGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTCGCCACGCCAGCt                                             >  1:876608/1‑68 (MQ=255)
                 tatCACGGCGGAACGCTTTACGTTGCGCGCGATTGAGCTTCGCCAACg                                                             >  1:114368/1‑48 (MQ=255)
                 tatCACGGCGGAACGCTTTACGTTGCGCGCGATTGAGCTTCGCCAACg                                                             >  1:559092/1‑48 (MQ=255)
                 tatCACGGCGGAACGCTTTACGTTGCGCGCGATTGAGCTTCGCCAACg                                                             >  1:430568/1‑48 (MQ=255)
                 tatCACGGCGGAACGCTTTACGTTGCGCGCGATTGAGCTTCGCCAACg                                                             >  1:1244141/1‑48 (MQ=255)
                 tatCACGGCGGAACGCTTTACGTTGCGCGCGATTGAGCTTCGCCAACg                                                             >  1:1389446/1‑48 (MQ=255)
                 tatCACGGCGGAACGCTTTACGTTGCGCGCGATTGAGCTTCGCCAACg                                                             >  1:1506691/1‑48 (MQ=255)
                 tatCACGGCGGAACGCTTTACGTTGCGCGCGATTGAGCTTCGCCAACg                                                             >  1:1555281/1‑48 (MQ=255)
                 tatCACGGCGGAACGCTTTACGTTGCGCGCGATTGAGCTTCGCCAACg                                                             >  1:1682117/1‑48 (MQ=255)
                 tatCACGGCGGAACGCTTTACGTTGCGCGCGATTGAGCTTCGCCAACg                                                             >  1:1702888/1‑48 (MQ=255)
                 tatCACGGCGGAACGCTTTACGTTGCGCGCGATTGAGCTTCGCCAACg                                                             >  1:1907399/1‑48 (MQ=255)
                 tatCACGGCGGAACGCTTTACGTTGCGCGCGATTGAGCTTCGCCAACg                                                             >  1:648420/1‑48 (MQ=255)
                 tatCACGGCGGAACGCTTTACGTTGCGCGCGATTGAGCTTCGCCAACg                                                             >  1:2043658/1‑48 (MQ=255)
                 tatCACGGCGGAACGCTTTACGTTGCGCGCGATTGAGCTTCGCCAACg                                                             >  1:730092/1‑48 (MQ=255)
                        gcggAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTcgccacgc                                                  >  1:32005/1‑52 (MQ=255)
                        gcggAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTCGCCACGCCAGCTAATATTCCCCt                                  >  1:2070171/1‑68 (MQ=255)
                                                  tGAGTTTCGCCAGCGGTTCGCCACGCCAGCTAATATTCCCCTGCGCGGGCGATTCTAAAcc               <  1:383514/61‑1 (MQ=255)
                                                   gAGTTTCGCCAGCGGGTCGCCACGCCAGCTAATATTCCCCTGCGCGGGCGATTCTAAAcc               <  1:1218795/60‑1 (MQ=255)
                                                        tCGCCAGCGGTTCGCCACGCCAGCTAATATTCCCCTGCGCGGGCGATTCTAAACCCACCAGCAACcgc  >  1:2597910/1‑68 (MQ=255)
                                                          gCCAACGGTTCGCCGCGCCAGCTAATATTTCCCTGTTAGGGCGACTCTAAccct              <  1:2320981/54‑2 (MQ=255)
                                                          gCCAACGGTTCGCCGCGCCAGCTAATATTTCCCTGTGAGGGCGACTCTAAccct              <  1:3417437/54‑2 (MQ=255)
                                                          gCCAACGGTTCGCCGCGCCAGCTAATATTTCCCTGTGAGGGCGACTCTAAccct              <  1:3279106/54‑2 (MQ=255)
                                                          gCCAACGGTTCGCCGCGCCAGCTAATATTTCCCTGTGAGGGCGACTCTAAccct              <  1:299/54‑2 (MQ=255)
                                                          gCCAACGGTTCGCCGCGCCAGCTAATATTTCCCTGTGAGGGCGACTCTAAccct              <  1:765125/54‑2 (MQ=255)
                                                          gCCAACGGTTCGCCGCGCCAGCTAATATTTCCCTGTGAGGGCGACTCTAAccct              <  1:2124894/54‑2 (MQ=255)
                                                              |                                                             
CTGAAATACCATCTGAATATCGCGGCGGAACGCTTTACGCTGGGCGCGATTGAGTTTCGCCAGCGGTTCGCCACGCCAGCTAATATTCCCCTGCGCGGGCGATTCTAAACCCACCAGCAACCGC  >  minE/2609126‑2609249

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: