Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,608 T→C 26.1% S1022S (TCT→TCC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,6080TC26.1% 107.9 / 43.7 69S1022S (TCT→TCCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (33/18);  new base C (7/11);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 9.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

ATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACT  >  minE/834538‑834676
                                                                      |                                                                    
atatTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCt                                                                      <  1:3872957/71‑1 (MQ=255)
 tatTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCtt                                                                     <  1:1375015/71‑1 (MQ=255)
     cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTa                                                                    >  1:653439/1‑68 (MQ=255)
     cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTa                                                                    >  1:3752800/1‑68 (MQ=255)
     cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa                                                                  <  1:1228092/70‑1 (MQ=255)
      cggCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTAc                                                                   <  1:550323/67‑1 (MQ=255)
      aGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAAc                                                                >  1:1431713/1‑71 (MQ=255)
        gcgcaGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACg                                                               <  1:3136073/67‑1 (MQ=255)
          acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGcc                                                             <  1:710114/70‑1 (MQ=255)
          acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg                                                            >  1:3223269/1‑71 (MQ=255)
          acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg                                                            >  1:526679/1‑71 (MQ=255)
          acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg                                                            >  1:3656926/1‑71 (MQ=255)
                  gcAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACAcc                                                    <  1:1113689/71‑1 (MQ=255)
                   taaaaTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCggga                                                          <  1:3233394/63‑2 (MQ=255)
                   cAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTacac                                                     >  1:2292212/1‑69 (MQ=255)
                    aaaaTTATCTCCAGCCCGACCTGGTCTGGTCTGGATGATGAACACGTTTCTTACAACGCCGGTTACACCaa                                                  >  1:1759730/1‑71 (MQ=255)
                    aaaaTTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCtt                                                                     >  1:3949103/1‑52 (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:2390005/71‑1 (MQ=255)
                       aTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGt                                               <  1:2394697/71‑1 (MQ=255)
                        ttATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATa                                                        <  1:3339452/61‑1 (MQ=255)
                        ttATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATa                                                        <  1:2956570/61‑1 (MQ=255)
                          aTCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAAc                                                                >  1:1464103/1‑51 (MQ=255)
                           tctcCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCgg                                                           >  1:49419/1‑55 (MQ=255)
                               cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGa                                         <  1:3493738/69‑1 (MQ=255)
                               cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGc                                       >  1:1885374/1‑71 (MQ=255)
                                aGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCt                                                                      >  1:898653/1‑39 (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)
                                           gtctggtctgGAAGATGAACACGTTTCTTACAACGc                                                              <  1:2911476/36‑1 (MQ=255)
                                           gtccggtctgGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTgg                            <  1:3774248/70‑1 (MQ=255)
                                               ggtctggAAGATTAACACGTTTCTTACAACGCCGGTTa                                                        <  1:3135240/38‑1 (MQ=255)
                                                gtctggAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACg                      >  1:578677/1‑71 (MQ=255)
                                                gtctggAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACg                      >  1:1977973/1‑71 (MQ=255)
                                                 tctggAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACg                      <  1:3428785/70‑1 (MQ=255)
                                                     gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctct                  <  1:1625502/70‑1 (MQ=255)
                                                     gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctct                  <  1:1248991/70‑1 (MQ=255)
                                                     gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctct                  <  1:887036/70‑1 (MQ=255)
                                                       aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg               >  1:3331147/1‑71 (MQ=255)
                                                       aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg               >  1:3472485/1‑71 (MQ=255)
                                                       aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg               >  1:587615/1‑71 (MQ=255)
                                                       aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg               >  1:1127115/1‑71 (MQ=255)
                                                       aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg               >  1:1381597/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:1314246/1‑54 (MQ=255)
                                                             acacGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                          >  1:2481734/1‑54 (MQ=255)
                                                              cacGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctc                   >  1:189805/1‑60 (MQ=255)
                                                              cacGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGctgtc                 >  1:826260/1‑62 (MQ=255)
                                                               acGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTa                        >  1:1628715/1‑54 (MQ=255)
                                                               acGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTa                        >  1:2798653/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)
                                                                      |                                                                    
ATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACT  >  minE/834538‑834676

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: