Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,469 G→A 6.8% A309A (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,4690GA 6.8% 266.2 / 7.3 88A309A (GCG→GCAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (47/35);  new base A (4/2);  total (51/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

AAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTG  >  minE/832402‑832535
                                                                   |                                                                  
aaaTCGCCTAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                                  <  1:1715310/70‑1 (MQ=255)
aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                                  <  1:3213067/70‑1 (MQ=255)
aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                                  <  1:966293/70‑1 (MQ=255)
aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                                  <  1:2522704/70‑1 (MQ=255)
aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                                  <  1:2973687/70‑1 (MQ=255)
aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                                  <  1:1942153/70‑1 (MQ=255)
 aaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATgg                                                                >  1:1098882/1‑71 (MQ=255)
  aTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggg                                                               >  1:2943493/1‑71 (MQ=255)
   tCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGc                                                              <  1:3717259/71‑1 (MQ=255)
     gCCAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGcc                                                             <  1:1701434/70‑1 (MQ=255)
     gCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCa                                                            >  1:3053525/1‑71 (MQ=255)
      ccAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCat                                                           <  1:418280/71‑2 (MQ=255)
       cAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCatg                                                          >  1:2124485/1‑69 (MQ=255)
        aaaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                                         <  1:3310625/71‑1 (MQ=255)
        aaaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                                         <  1:1055514/71‑1 (MQ=255)
         aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCAcct                                                         >  1:1642185/1‑68 (MQ=255)
         aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACg                                                          >  1:1401032/1‑69 (MQ=255)
         aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                                         >  1:1703189/1‑70 (MQ=255)
         aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                                         >  1:3090243/1‑70 (MQ=255)
         aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                                         >  1:1569310/1‑70 (MQ=255)
            tgtgCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGTGAtt                                                     >  1:2941500/1‑70 (MQ=255)
             gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgc                                                    >  1:1987441/1‑71 (MQ=255)
             gtgCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGc                                                              >  1:2269038/1‑61 (MQ=255)
              tgCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGTGAttt                                                    >  1:596674/1‑68 (MQ=255)
                cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc                                                 <  1:3913534/71‑1 (MQ=255)
                cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc                                                 <  1:2682659/71‑1 (MQ=255)
                  aTCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGt                                               >  1:3609478/1‑71 (MQ=255)
                  aTCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGt                                               >  1:3799497/1‑71 (MQ=255)
                   tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCg                                                >  1:2492217/1‑69 (MQ=255)
                   tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGt                                               >  1:1233680/1‑70 (MQ=255)
                     tgtgGCTGGCACCGAAACATGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAt                                                      >  1:2915741/1‑61 (MQ=255)
                       tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCg                                                <  1:2416428/65‑1 (MQ=255)
                       tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa                                             <  1:986865/68‑1 (MQ=255)
                       tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGACACGTAATGCTGCGTGa                                             <  1:2115255/68‑1 (MQ=255)
                       tggctgacACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAAtt                                          >  1:802/1‑71 (MQ=255)
                        ggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTc                                         >  1:2049401/1‑71 (MQ=255)
                           tggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCAc                                      >  1:2446409/1‑71 (MQ=255)
                             gcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCt                                    <  1:2313349/71‑1 (MQ=255)
                             gcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCt                                    >  1:765298/1‑71 (MQ=255)
                                ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa                                                                   <  1:2477854/37‑1 (MQ=255)
                                  gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACa                               <  1:843796/71‑1 (MQ=255)
                                  gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACa                               <  1:2194359/71‑1 (MQ=255)
                                   aaaCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACaa                              >  1:1302991/1‑71 (MQ=255)
                                   aaaCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACaa                              >  1:1622538/1‑71 (MQ=255)
                                    aaCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCt                                    >  1:2583947/1‑64 (MQ=255)
                                      cAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAAccc                           <  1:3215379/71‑1 (MQ=255)
                                      cAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAAcc                            >  1:2465597/1‑70 (MQ=255)
                                        gggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCaa                         <  1:918481/71‑1 (MQ=255)
                                         ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg                        >  1:2616635/1‑71 (MQ=255)
                                         ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg                        >  1:1667857/1‑71 (MQ=255)
                                         ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg                        >  1:1619693/1‑71 (MQ=255)
                                         ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg                        >  1:1755151/1‑71 (MQ=255)
                                         ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg                        >  1:2512421/1‑71 (MQ=255)
                                          gCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc                       >  1:3395877/1‑71 (MQ=255)
                                          gCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc                       >  1:1785149/1‑71 (MQ=255)
                                           cACCGATGCGGCAATGGCGCTTGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                      <  1:381264/71‑1 (MQ=255)
                                           cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg                        <  1:2637359/69‑1 (MQ=255)
                                           cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                      <  1:1731723/71‑1 (MQ=255)
                                             ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg                    >  1:2663934/1‑71 (MQ=255)
                                             ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg                    >  1:3907090/1‑71 (MQ=255)
                                             ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg                    >  1:313832/1‑71 (MQ=255)
                                              cGATGCGGCAATGGCGCTGGCGATGGGCCACGTaa                                                       >  1:350319/1‑35 (MQ=255)
                                              cGATGCGGCAATGGCGCTGGCGATGGGCCACGTaa                                                       >  1:1561248/1‑35 (MQ=255)
                                                aTGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGa                                 >  1:93241/1‑55 (MQ=255)
                                                aTGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGa                                 >  1:347961/1‑55 (MQ=255)
                                                aTGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt                 <  1:2626359/71‑1 (MQ=255)
                                                aTGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt                 <  1:2622479/71‑1 (MQ=255)
                                                 tGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAtt                <  1:1023727/71‑1 (MQ=255)
                                                  gCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc                                   <  1:182222/51‑1 (MQ=255)
                                                  gCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                                <  1:2418178/54‑1 (MQ=255)
                                                   cGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAttt               <  1:1676676/70‑1 (MQ=255)
                                                   cGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAttt               <  1:1953701/70‑1 (MQ=255)
                                                    ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAttt               <  1:547279/69‑1 (MQ=255)
                                                         tGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGAc        >  1:1980104/1‑71 (MQ=255)
                                                         tGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGAc        >  1:535717/1‑71 (MQ=255)
                                                         tGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGAc        >  1:317888/1‑71 (MQ=255)
                                                           gcgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg                        >  1:358105/1‑53 (MQ=255)
                                                           gcgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg                        >  1:1924440/1‑53 (MQ=255)
                                                           gcgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACt       <  1:530625/70‑1 (MQ=255)
                                                            cgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAt     <  1:1191480/71‑1 (MQ=255)
                                                             gcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAt     <  1:717870/70‑1 (MQ=255)
                                                              cTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAtgt   <  1:2270230/71‑1 (MQ=255)
                                                              cTGGCGATGGGCCACGTAATGATGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAtgt   >  1:245330/1‑71 (MQ=255)
                                                               tGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACaa                              >  1:2297936/1‑43 (MQ=255)
                                                               tGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAtgtg  >  1:3803149/1‑71 (MQ=255)
                                                               tGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTAACCGACTAtgtg  >  1:3011906/1‑71 (MQ=255)
                                                                ggCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGa                                 >  1:3743568/1‑39 (MQ=255)
                                                                ggCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGa                                 >  1:429052/1‑39 (MQ=255)
                                                                   gATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCa                     >  1:1639657/1‑48 (MQ=255)
                                                                   |                                                                  
AAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTG  >  minE/832402‑832535

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: