Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,628 T→G 18.0% V362V (GTT→GTG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,6280TG18.0% 144.0 / 26.3 83V362V (GTT→GTGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (16/52);  new base G (6/9);  total (22/61)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGA  >  minE/832560‑832690
                                                                    |                                                              
cTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGtt                                                                <  1:432087/69‑1 (MQ=255)
cTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGtt                                                                <  1:1817463/69‑1 (MQ=255)
cTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGtt                                                                <  1:1215052/69‑1 (MQ=255)
   gTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGa                                                              >  1:538797/1‑68 (MQ=255)
   gTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGAt                                                             >  1:1142130/1‑69 (MQ=255)
   gTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGAt                                                             >  1:2289415/1‑69 (MQ=255)
   gTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGAt                                                             >  1:640874/1‑69 (MQ=255)
     gATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATgc                                                           <  1:193554/69‑1 (MQ=255)
     gATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATgc                                                           <  1:1849092/69‑1 (MQ=255)
     gATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATgc                                                           >  1:1447654/1‑69 (MQ=255)
     gATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATgc                                                           >  1:584735/1‑69 (MQ=255)
     gATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATgc                                                           >  1:634751/1‑69 (MQ=255)
     gATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATgc                                                           >  1:2243245/1‑69 (MQ=255)
     gATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATgc                                                           <  1:1238352/69‑1 (MQ=255)
     gATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATgc                                                           >  1:2371211/1‑69 (MQ=255)
     gATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATgc                                                           >  1:2566974/1‑69 (MQ=255)
     gATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATgc                                                           >  1:2687359/1‑69 (MQ=255)
     gATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATgc                                                           >  1:2669261/1‑69 (MQ=255)
     gATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTATCTGGTTGATgcg                                                          >  1:598906/1‑69 (MQ=255)
                            actacGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCgg                                    <  1:1803660/69‑1 (MQ=255)
                            actacGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCgg                                    <  1:1848740/69‑1 (MQ=255)
                            actacGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCgg                                    <  1:1740973/69‑1 (MQ=255)
                            actacGCAGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCgg                                    <  1:773562/69‑1 (MQ=255)
                              tacGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCgg                                    <  1:767428/67‑1 (MQ=255)
                                    gCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCgg                                    <  1:1855441/61‑1 (MQ=255)
                                    gCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGaaa                            <  1:874646/69‑1 (MQ=255)
                                       ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:2082442/69‑1 (MQ=255)
                                       ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:1260111/69‑1 (MQ=255)
                                       ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:3326041/69‑1 (MQ=255)
                                       ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:2108312/69‑1 (MQ=255)
                                       ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:2051810/69‑1 (MQ=255)
                                       ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:1298182/69‑1 (MQ=255)
                                       ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:1641068/69‑1 (MQ=255)
                                       ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:2658406/69‑1 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCTAGGAAAACAATCCGGAAt                                 <  1:2931530/54‑1 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:1054/54‑1 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:2718345/54‑1 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:2786182/54‑1 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:672036/54‑1 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:3269638/54‑1 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:360520/54‑1 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:501272/54‑1 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:475666/54‑1 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:934493/54‑1 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:1626137/54‑1 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 >  1:2682063/1‑54 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:1749417/54‑1 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 >  1:2202315/1‑54 (MQ=255)
                                              tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:1258142/54‑1 (MQ=255)
                                              tgctgcGCGCTGATGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 >  1:61823/1‑54 (MQ=255)
                                                  gcgcgcTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:2594031/69‑1 (MQ=255)
                                                  gcgcgcTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:142769/69‑1 (MQ=255)
                                                  gcgcgcTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:1356712/69‑1 (MQ=255)
                                                  gcgcgcTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:115349/69‑1 (MQ=255)
                                                  gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCtttga              <  1:1326384/69‑3 (MQ=255)
                                                  gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:614296/69‑1 (MQ=255)
                                                  gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:1372138/69‑1 (MQ=255)
                                                  gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:914634/69‑1 (MQ=255)
                                                  gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:779403/69‑1 (MQ=255)
                                                  gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:3331655/69‑1 (MQ=255)
                                                  gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:3166844/69‑1 (MQ=255)
                                                  gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:1220021/69‑1 (MQ=255)
                                                   cgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:1591980/68‑1 (MQ=255)
                                                      gctgctGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATAcc          >  1:604717/1‑69 (MQ=255)
                                                      gctgctGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATAcc          >  1:21499/1‑69 (MQ=255)
                                                      gctgctGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATAcc          >  1:472635/1‑69 (MQ=255)
                                                      gctgctGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATAcc          >  1:2913763/1‑69 (MQ=255)
                                                      gctgctGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATAcc          >  1:480987/1‑69 (MQ=255)
                                                      gctgctGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTATAATAcc          >  1:1113425/1‑69 (MQ=255)
                                                         gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:2572194/60‑1 (MQ=255)
                                                         gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:434449/60‑1 (MQ=255)
                                                         gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:3412648/60‑1 (MQ=255)
                                                         gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:2718572/60‑1 (MQ=255)
                                                         gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:2532449/60‑1 (MQ=255)
                                                         gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:3228026/60‑1 (MQ=255)
                                                         gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:2453296/60‑1 (MQ=255)
                                                         gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:1172207/60‑1 (MQ=255)
                                                         gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:3001812/60‑1 (MQ=255)
                                                         gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:3027089/60‑1 (MQ=255)
                                                               cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:64408/68‑1 (MQ=255)
                                                               cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:646898/68‑1 (MQ=255)
                                                               cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:496091/68‑1 (MQ=255)
                                                               cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:3217169/68‑1 (MQ=255)
                                                               cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:2227590/68‑1 (MQ=255)
                                                               cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:3387324/68‑1 (MQ=255)
                                                               cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:1735373/68‑1 (MQ=255)
                                                                    |                                                              
CTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGA  >  minE/832560‑832690

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: