Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,620 T→A 28.6% G1026G (GGT→GGA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,6200TA28.6% 102.6 / 54.5 70G1026G (GGT→GGAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (28/22);  new base A (7/13);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

GCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAC  >  minE/834556‑834689
                                                                |                                                                     
gcAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACAcc                                                                 <  1:1113689/71‑1 (MQ=255)
 cAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTacac                                                                  >  1:2292212/1‑69 (MQ=255)
  aaaaTTATCTCCAGCCCGACCTGGTCTGGTCTGGATGATGAACACGTTTCTTACAACGCCGGTTACACCaa                                                               >  1:1759730/1‑71 (MQ=255)
   aaaTTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAAc                                                              >  1:3283398/1‑71 (MQ=255)
   aaaTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAAc                                                              >  1:3884871/1‑71 (MQ=255)
    aaTTATCTCCAGCCCGCCCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAAc                                                              >  1:2240582/1‑70 (MQ=255)
     aTTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTa                                                                     >  1:3101921/1‑62 (MQ=255)
     aTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGt                                                            <  1:2394697/71‑1 (MQ=255)
     aTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGt                                                            <  1:2390005/71‑1 (MQ=255)
      ttATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATa                                                                     <  1:3339452/61‑1 (MQ=255)
      ttATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATa                                                                     <  1:2956570/61‑1 (MQ=255)
             cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGa                                                      <  1:3493738/69‑1 (MQ=255)
             cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGc                                                    >  1:1885374/1‑71 (MQ=255)
                 ccGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGAt                                                >  1:2931444/1‑71 (MQ=255)
                  cGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATc                                               >  1:2420899/1‑71 (MQ=255)
                     ccTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                                            <  1:997730/71‑1 (MQ=255)
                      ctggtccggtCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCat                                                        <  1:2195101/58‑2 (MQ=255)
                         gtctggtctgGAAGATGAACACTTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                                            <  1:2066342/67‑1 (MQ=255)
                         gtccggtctgGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTgg                                         <  1:3774248/70‑1 (MQ=255)
                             ggtctggAAGATTAACACGTTTCTTACAACGCCGGTTa                                                                     <  1:3135240/38‑1 (MQ=255)
                              gtctggAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACg                                   >  1:1977973/1‑71 (MQ=255)
                              gtctggAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACg                                   >  1:578677/1‑71 (MQ=255)
                               tctggAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACg                                   <  1:3428785/70‑1 (MQ=255)
                                   gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctct                               <  1:1248991/70‑1 (MQ=255)
                                   gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctct                               <  1:887036/70‑1 (MQ=255)
                                   gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctct                               <  1:1625502/70‑1 (MQ=255)
                                     aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                            >  1:3331147/1‑71 (MQ=255)
                                     aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                            >  1:3472485/1‑71 (MQ=255)
                                     aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                            >  1:1127115/1‑71 (MQ=255)
                                     aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                            >  1:1381597/1‑71 (MQ=255)
                                     aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                            >  1:587615/1‑71 (MQ=255)
                                      gATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTgg                           <  1:893974/71‑1 (MQ=255)
                                      gATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTgg                           <  1:2667360/71‑1 (MQ=255)
                                       aTGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                                            <  1:897873/53‑1 (MQ=255)
                                        tGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtc                         >  1:2139387/1‑71 (MQ=255)
                                         gAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATg                                          >  1:667060/1‑53 (MQ=255)
                                          aaCACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtcgt                       <  1:552425/71‑1 (MQ=255)
                                          aaCACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtcgt                       <  1:1873440/71‑1 (MQ=255)
                                           tcacGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGctgt                               <  1:1475867/61‑3 (MQ=255)
                                           tcacGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGgacgc                       >  1:1666397/2‑69 (MQ=255)
                                           acacGTTTCTTACAACGCCGGTTACACCAACGTTCa                                                         >  1:3526942/1‑36 (MQ=255)
                                           acacGTTTCTTACAACGCCGGTTACACCAACGTTCACGAg                                                     >  1:1072750/1‑40 (MQ=255)
                                           acacGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                                       >  1:2481734/1‑54 (MQ=255)
                                           acacGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                                       >  1:1314246/1‑54 (MQ=255)
                                            cacGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctc                                >  1:189805/1‑60 (MQ=255)
                                            cacGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGctgtc                              >  1:826260/1‑62 (MQ=255)
                                             acGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTa                                     >  1:2798653/1‑54 (MQ=255)
                                             acGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTa                                     >  1:1628715/1‑54 (MQ=255)
                                             acGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTa                                     >  1:1998141/1‑54 (MQ=255)
                                             acGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGgacgc                       <  1:159881/68‑2 (MQ=255)
                                             acGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAg                    >  1:1962048/1‑71 (MQ=255)
                                              cGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTg                                                  <  1:1888435/40‑1 (MQ=255)
                                                  tCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCa                     >  1:3464236/1‑65 (MQ=255)
                                                  tCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCa                     >  1:3565001/1‑65 (MQ=255)
                                                  tCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACt               >  1:1259746/1‑71 (MQ=255)
                                                     tataaCGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCa                     <  1:2323948/59‑1 (MQ=255)
                                                      aCAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTa            <  1:2709250/70‑1 (MQ=255)
                                                       taaCGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTa            <  1:3553339/68‑1 (MQ=255)
                                                        aaCGCCGGTTACACCAACGTTCACGAGCTGATCCCAt                                           <  1:1515248/37‑1 (MQ=255)
                                                        aaCGCCGGTTACACCAACGTTCACGAGCTGATCCCAt                                           <  1:1747078/37‑1 (MQ=255)
                                                        aaCGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTAt           >  1:2042326/1‑69 (MQ=255)
                                                         aCGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCa         <  1:965805/70‑1 (MQ=255)
                                                         aCGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCAg        >  1:807224/1‑71 (MQ=255)
                                                          cgcgGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGgacgc                       <  1:281607/55‑2 (MQ=255)
                                                          cGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                 <  1:2240643/61‑1 (MQ=255)
                                                          cGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAgg       <  1:2041313/71‑1 (MQ=255)
                                                             cGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGAt     >  1:1650743/1‑70 (MQ=255)
                                                              ggTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTAt           <  1:2732990/63‑1 (MQ=255)
                                                              ggTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTAt           <  1:3104310/63‑1 (MQ=255)
                                                               gTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctc                                <  1:2077042/41‑1 (MQ=255)
                                                                ttACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAc  >  1:1316801/1‑70 (MQ=255)
                                                                |                                                                     
GCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAC  >  minE/834556‑834689

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: