Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,367 G→A 48.5% G275G (GGG→GGA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,3670GA48.5% 4.2 / 93.4 105G275G (GGG→GGAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base A (34/17);  total (60/45)
Fisher's exact test for biased strand distribution p-value = 7.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

AAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACT  >  minE/832307‑832414
                                                            |                                               
aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGAGCAAAAc                                         <  1:6115/69‑1 (MQ=255)
aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAAc                                         <  1:2088056/69‑1 (MQ=255)
aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAAc                                         <  1:2062863/69‑1 (MQ=255)
aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAAc                                         <  1:909060/69‑1 (MQ=255)
aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAAc                                         <  1:1243684/69‑1 (MQ=255)
aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAAc                                         <  1:840286/69‑1 (MQ=255)
aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAAc                                         <  1:686738/69‑1 (MQ=255)
aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAAc                                         <  1:524416/69‑1 (MQ=255)
     tGCCGCAGACGCGTACCCCGGATGCTCTCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTg                                    >  1:3305942/1‑69 (MQ=255)
     tGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTTCAAAGGGACCAAAACTGtt                                     >  1:2454683/1‑68 (MQ=255)
     tGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTg                                    >  1:2033709/1‑69 (MQ=255)
     tGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTg                                    >  1:3124900/1‑69 (MQ=255)
         gCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                                <  1:977691/69‑1 (MQ=255)
         gCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                                <  1:1932340/69‑1 (MQ=255)
         gCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                                <  1:2475726/69‑1 (MQ=255)
         gCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                                <  1:2532552/69‑1 (MQ=255)
         gCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                                <  1:750612/69‑1 (MQ=255)
         gCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                                <  1:1444193/69‑1 (MQ=255)
         acAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                                <  1:827331/68‑1 (MQ=255)
           aGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                              <  1:1194807/69‑1 (MQ=255)
           aGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                              <  1:349534/69‑1 (MQ=255)
           aGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                              <  1:2028946/69‑1 (MQ=255)
           aGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGCTGCCGTca                              <  1:2056573/69‑1 (MQ=255)
            gACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTcac                             <  1:1599003/69‑1 (MQ=255)
            gACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTcac                             <  1:1707574/69‑1 (MQ=255)
               gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:2673479/1‑66 (MQ=255)
               gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:2750783/1‑66 (MQ=255)
               gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:365146/1‑66 (MQ=255)
               gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:2796466/1‑66 (MQ=255)
               gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:2119523/1‑66 (MQ=255)
               gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:2888164/1‑66 (MQ=255)
               gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:854136/1‑66 (MQ=255)
               gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:2902869/1‑66 (MQ=255)
               gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:3252924/1‑66 (MQ=255)
               gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:2739729/1‑66 (MQ=255)
               gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:1344615/1‑66 (MQ=255)
               gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:1423912/1‑66 (MQ=255)
               gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:1516648/1‑66 (MQ=255)
               gcgTACCCCGGATGATCACTTCTTTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTcac                             >  1:1441239/1‑66 (MQ=255)
                           tGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGa              >  1:1292526/1‑69 (MQ=255)
                                   tcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTa                    >  1:1419748/1‑55 (MQ=255)
                                   tcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa      >  1:1768132/1‑69 (MQ=255)
                                   tcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa      >  1:2089889/1‑69 (MQ=255)
                                   tcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa      >  1:33358/1‑69 (MQ=255)
                                   tcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa      >  1:1467892/1‑69 (MQ=255)
                                   tcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa      >  1:2539714/1‑69 (MQ=255)
                                   tcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa      >  1:1782922/1‑69 (MQ=255)
                                   tcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa      >  1:1409941/1‑69 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaa     >  1:379062/1‑69 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaa     >  1:1782694/1‑69 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaa     >  1:755842/1‑69 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaa     >  1:1656404/1‑69 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGt                                >  1:3305869/1‑42 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACAcc                          >  1:791005/1‑48 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:398957/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:1357882/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:416831/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:477813/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:75388/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:788153/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:2141559/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:2238312/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:2348460/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:9678/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:2556660/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:2810869/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:2650034/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:270589/1‑49 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGAACCAAAACTCTTGCCGTCACACCa                         >  1:3169715/1‑49 (MQ=255)
                                    cttTACTAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:2680167/1‑48 (MQ=255)
                                    cattACTGAAGTGCGTTACAAAGGAACCAAAACTGTTGCCGTCACACCa                         >  1:1206365/3‑49 (MQ=255)
                                       tACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACt  >  1:2926200/1‑69 (MQ=255)
                                       tACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACt  >  1:1495762/1‑69 (MQ=255)
                                       tACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACt  >  1:1717963/1‑69 (MQ=255)
                                       tACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACt  >  1:1663495/1‑69 (MQ=255)
                                       tACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACt  >  1:851308/1‑69 (MQ=255)
                                       tACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACt  >  1:2885438/1‑69 (MQ=255)
                                       tACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACt  >  1:2946457/1‑69 (MQ=255)
                                                         aGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACt  <  1:385741/51‑1 (MQ=255)
                                                         aGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACt  <  1:1799617/51‑1 (MQ=255)
                                                         aGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACt  <  1:2936338/51‑1 (MQ=255)
                                                         aGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACt  <  1:1384482/51‑1 (MQ=255)
                                                         aGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACt  <  1:1086417/51‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:878576/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:3176624/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:614783/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:60994/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:598878/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:5479/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:3227062/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:1605735/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:2982893/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:1609619/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:2828730/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:2811246/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:2780923/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:2760470/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:1002081/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt           <  1:1994566/42‑1 (MQ=255)
                                                         aGGAACCAAAACTGTTGACGTCACACCAGACTACGCTGAAAt           <  1:2357290/42‑1 (MQ=255)
                                                            |                                               
AAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACT  >  minE/832307‑832414

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: