Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,448 C→T 7.3% T302T (ACC→ACT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,4480CT 7.3% 247.3 / 6.9 82T302T (ACC→ACTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (40/36);  new base T (4/2);  total (44/38)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

TTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAG  >  minE/832378‑832517
                                                                      |                                                                     
ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                                                       >  1:3620646/1‑71 (MQ=255)
    cGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGAt                                                                    <  1:350066/70‑1 (MQ=255)
    cGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATg                                                                   <  1:708412/71‑1 (MQ=255)
     gTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGc                                                                  <  1:1945686/71‑1 (MQ=255)
        acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGc                                                                  >  1:1565246/1‑68 (MQ=255)
        acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                                                                >  1:1494829/1‑70 (MQ=255)
            cAGACTACGCTGCAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAGACAGGGCACCGATGCGGCAATg                                                           >  1:1818749/1‑71 (MQ=255)
            cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                                           <  1:2466817/71‑1 (MQ=255)
            cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                                           >  1:1469432/1‑71 (MQ=255)
            cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                                           >  1:384313/1‑71 (MQ=255)
             aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                          <  1:2809393/71‑1 (MQ=255)
              gACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                          >  1:3039408/1‑70 (MQ=255)
                 tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                                                         >  1:2164379/1‑68 (MQ=255)
                 tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                                                         >  1:3805538/1‑68 (MQ=255)
                 tACGCTGAAATCGCCAAACTGTGAGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                                                         >  1:1175749/1‑68 (MQ=255)
                   cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                                                      <  1:2326459/69‑1 (MQ=255)
                    gCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGc                                                   <  1:1074059/71‑1 (MQ=255)
                        aaaTCGCCTAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                <  1:1715310/70‑1 (MQ=255)
                        aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                <  1:2522704/70‑1 (MQ=255)
                        aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                <  1:2973687/70‑1 (MQ=255)
                        aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                <  1:3213067/70‑1 (MQ=255)
                        aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                <  1:966293/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)
                            cGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg                                                    >  1:1321583/1‑62 (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:1055514/71‑1 (MQ=255)
                                aaaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                       <  1:3310625/71‑1 (MQ=255)
                                 aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                                                                >  1:3767962/1‑45 (MQ=255)
                                 aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCAcct                                       >  1:1642185/1‑68 (MQ=255)
                                 aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACg                                        >  1:1401032/1‑69 (MQ=255)
                                 aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                       >  1:1569310/1‑70 (MQ=255)
                                 aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                       >  1:3090243/1‑70 (MQ=255)
                                 aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                       >  1:1703189/1‑70 (MQ=255)
                                    tgtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                          <  1:2401328/48‑1 (MQ=255)
                                    tgtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                          <  1:2405978/48‑1 (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)
                                         gATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                          <  1:1933014/43‑1 (MQ=255)
                                          aTCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGt                             >  1:3799497/1‑71 (MQ=255)
                                          aTCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGt                             >  1:3609478/1‑71 (MQ=255)
                                           tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                                                      <  1:3353623/45‑1 (MQ=255)
                                           tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                                                      <  1:1039792/45‑1 (MQ=255)
                                           tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                                                      <  1:3653872/45‑1 (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:765298/1‑71 (MQ=255)
                                                     gcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCt                  <  1:2313349/71‑1 (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:2512421/1‑71 (MQ=255)
                                                                 ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg      >  1:1755151/1‑71 (MQ=255)
                                                                 ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg      >  1:1667857/1‑71 (MQ=255)
                                                                 ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg      >  1:1619693/1‑71 (MQ=255)
                                                                 ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg      >  1:2616635/1‑71 (MQ=255)
                                                                  gCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc     >  1:1785149/1‑71 (MQ=255)
                                                                  gCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc     >  1:3395877/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:3907090/1‑71 (MQ=255)
                                                                     ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg  >  1:313832/1‑71 (MQ=255)
                                                                     ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg  >  1:2663934/1‑71 (MQ=255)
                                                                      cGATGCGGCAATGGCGCTGGCGATGGGCCACGTaa                                     >  1:1561248/1‑35 (MQ=255)
                                                                      cGATGCGGCAATGGCGCTGGCGATGGGCCACGTaa                                     >  1:350319/1‑35 (MQ=255)
                                                                      |                                                                     
TTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAG  >  minE/832378‑832517

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: