Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,634 G→A 22.2% A364A (GCG→GCA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,6340GA22.2% 92.2 / 17.6 67A364A (GCG→GCAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (3/49);  new base A (6/9);  total (9/58)
Fisher's exact test for biased strand distribution p-value = 2.79e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGA  >  minE/832565‑832690
                                                                     |                                                        
gATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTATCTGGTTGATgcg                                                          >  1:598906/1‑69 (MQ=255)
                       actacGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCgg                                    <  1:1740973/69‑1 (MQ=255)
                       actacGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCgg                                    <  1:1803660/69‑1 (MQ=255)
                       actacGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCgg                                    <  1:1848740/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:1298182/69‑1 (MQ=255)
                                  ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:2051810/69‑1 (MQ=255)
                                  ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:1260111/69‑1 (MQ=255)
                                  ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:2082442/69‑1 (MQ=255)
                                  ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:3326041/69‑1 (MQ=255)
                                  ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:1641068/69‑1 (MQ=255)
                                  ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:2658406/69‑1 (MQ=255)
                                  ggTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACt                         <  1:2108312/69‑1 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCTAGGAAAACAATCCGGAAt                                 <  1:2931530/54‑1 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 >  1:2682063/1‑54 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:1054/54‑1 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:2718345/54‑1 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:672036/54‑1 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:2786182/54‑1 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:3269638/54‑1 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:501272/54‑1 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:360520/54‑1 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:475666/54‑1 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 >  1:2202315/1‑54 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:1749417/54‑1 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:934493/54‑1 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:1626137/54‑1 (MQ=255)
                                         tgctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 <  1:1258142/54‑1 (MQ=255)
                                         tgctgcGCGCTGATGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAAt                                 >  1:61823/1‑54 (MQ=255)
                                             gcgcgcTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:115349/69‑1 (MQ=255)
                                             gcgcgcTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:2594031/69‑1 (MQ=255)
                                             gcgcgcTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:142769/69‑1 (MQ=255)
                                             gcgcgcTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:1356712/69‑1 (MQ=255)
                                             gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCtttga              <  1:1326384/69‑3 (MQ=255)
                                             gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:779403/69‑1 (MQ=255)
                                             gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:3331655/69‑1 (MQ=255)
                                             gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:1220021/69‑1 (MQ=255)
                                             gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:3166844/69‑1 (MQ=255)
                                             gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:614296/69‑1 (MQ=255)
                                             gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:1372138/69‑1 (MQ=255)
                                             gcgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:914634/69‑1 (MQ=255)
                                              cgcgcTGCTGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa              <  1:1591980/68‑1 (MQ=255)
                                                 gctgctGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATAcc          >  1:604717/1‑69 (MQ=255)
                                                 gctgctGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATAcc          >  1:2913763/1‑69 (MQ=255)
                                                 gctgctGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATAcc          >  1:480987/1‑69 (MQ=255)
                                                 gctgctGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATAcc          >  1:21499/1‑69 (MQ=255)
                                                 gctgctGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATAcc          >  1:472635/1‑69 (MQ=255)
                                                 gctgctGATCTGGTGGATGCACTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTATAATAcc          >  1:1113425/1‑69 (MQ=255)
                                                    gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:3412648/60‑1 (MQ=255)
                                                    gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:1172207/60‑1 (MQ=255)
                                                    gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:2453296/60‑1 (MQ=255)
                                                    gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:2532449/60‑1 (MQ=255)
                                                    gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:2572194/60‑1 (MQ=255)
                                                    gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:2718572/60‑1 (MQ=255)
                                                    gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:3001812/60‑1 (MQ=255)
                                                    gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:3027089/60‑1 (MQ=255)
                                                    gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:3228026/60‑1 (MQ=255)
                                                    gctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                <  1:434449/60‑1 (MQ=255)
                                                          cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:64408/68‑1 (MQ=255)
                                                          cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:646898/68‑1 (MQ=255)
                                                          cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:3217169/68‑1 (MQ=255)
                                                          cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:3387324/68‑1 (MQ=255)
                                                          cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:496091/68‑1 (MQ=255)
                                                          cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:2227590/68‑1 (MQ=255)
                                                          cTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGa  <  1:1735373/68‑1 (MQ=255)
                                                                     |                                                        
GATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGA  >  minE/832565‑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: