Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,379 T→C 6.6% V279V (GTT→GTC) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,3790TC 6.6% 259.8 / 6.4 90V279V (GTT→GTC) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (45/39);  new base C (4/2);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

CGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACC  >  minE/832310‑832448
                                                                     |                                                                     
cGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTg                                                                      >  1:824326/1‑71 (MQ=255)
 gTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTg                                                                      >  1:1725713/1‑70 (MQ=255)
 gTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTg                                                                      >  1:2942999/1‑70 (MQ=255)
 gTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGc                                                                     >  1:654801/1‑71 (MQ=255)
 gTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGc                                                                     >  1:2328699/1‑71 (MQ=255)
  tGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGcc                                                                    >  1:350784/1‑71 (MQ=255)
   gCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCg                                                                   >  1:2668483/1‑71 (MQ=255)
    ccGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                                                                  <  1:5044771/71‑1 (MQ=255)
    ccGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                                                                  <  1:2411531/71‑1 (MQ=255)
           cgcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCa                                                           <  1:568557/71‑1 (MQ=255)
            gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCacac                                                             >  1:4665990/1‑68 (MQ=255)
              gTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                                                                  <  1:4666541/61‑1 (MQ=255)
              gTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGAc                                                        <  1:962101/71‑1 (MQ=255)
                 cgccGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTa                                                      >  1:1199696/3‑70 (MQ=255)
                    cGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTActct                                                  >  1:3479817/1‑71 (MQ=255)
                    cGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTActct                                                  >  1:2756726/1‑71 (MQ=255)
                    cGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTActc                                                   >  1:259288/1‑68 (MQ=255)
                     ggATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                                                                  <  1:4998232/54‑1 (MQ=255)
                     ggATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                                                                  <  1:284573/54‑1 (MQ=255)
                     ggATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGc                                                   >  1:1864301/1‑69 (MQ=255)
                     ggACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTACTCTg                                                 <  1:5512562/71‑1 (MQ=255)
                     ggACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTACTCTg                                                 <  1:2691524/71‑1 (MQ=255)
                      gATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTg                                                 >  1:689125/1‑70 (MQ=255)
                       aTGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaa                                               <  1:4310725/71‑1 (MQ=255)
                        tGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaaa                                              <  1:2307398/71‑1 (MQ=255)
                        tGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaa                                               >  1:3174580/1‑70 (MQ=255)
                          cTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                                                                >  1:3486655/1‑51 (MQ=255)
                          cTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                                                                >  1:5423924/1‑51 (MQ=255)
                          cTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                                                                >  1:5139528/1‑51 (MQ=255)
                           tCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCg                                           >  1:3303168/1‑71 (MQ=255)
                            cACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTa                                                      <  1:3966909/59‑1 (MQ=255)
                            cACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTa                                                      <  1:5307474/59‑1 (MQ=255)
                             aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGa                                                         >  1:5292437/1‑55 (MQ=255)
                             aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGc                                          >  1:2624974/1‑70 (MQ=255)
                             aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGc                                          >  1:2640918/1‑70 (MQ=255)
                               ttcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGcc                                                                    <  1:2940867/42‑1 (MQ=255)
                                tcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa                                      >  1:2717180/1‑71 (MQ=255)
                                 cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa                                      <  1:2724908/70‑1 (MQ=255)
                                 cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                     >  1:3395995/1‑71 (MQ=255)
                                 cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                     >  1:1454583/1‑71 (MQ=255)
                                 cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                     >  1:270543/1‑71 (MQ=255)
                                 cttTACTGAAGTGCGTTACAAAGGGACAAAAACTGTTGCCGTCACACCAGACTACGCTGAAATc                                            >  1:3524381/1‑64 (MQ=255)
                                    tACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtg                                   <  1:2003354/70‑1 (MQ=255)
                                     aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaa                                       <  1:4275270/65‑1 (MQ=255)
                                     aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgtg                                 >  1:2901675/1‑71 (MQ=255)
                                     aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgtg                                 >  1:5422094/1‑71 (MQ=255)
                                     aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt                                  >  1:88103/1‑70 (MQ=255)
                                      cTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGc                                <  1:5286533/71‑1 (MQ=255)
                                        gAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGa                              >  1:5670754/1‑71 (MQ=255)
                                         aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                     >  1:5057104/1‑63 (MQ=255)
                                         aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCg                               <  1:1822553/69‑1 (MQ=255)
                                         aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCg                               <  1:1851614/69‑1 (MQ=255)
                                          aGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATc                            <  1:1794494/71‑1 (MQ=255)
                                          aGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATc                            >  1:1294859/1‑71 (MQ=255)
                                            tGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGa                              <  1:517313/67‑1 (MQ=255)
                                            tGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCtg                          >  1:2701369/1‑71 (MQ=255)
                                              cGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaaa                                              >  1:988673/1‑49 (MQ=255)
                                                ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGtaa                                              <  1:2958244/47‑1 (MQ=255)
                                                ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaaa                                              <  1:2893078/47‑1 (MQ=255)
                                                ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaaa                                              <  1:95249/47‑1 (MQ=255)
                                                ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt                                  <  1:617241/59‑1 (MQ=255)
                                                ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt                                  <  1:1713914/59‑1 (MQ=255)
                                                  aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACTCTGAAATCGcc                                         >  1:4828489/1‑50 (MQ=255)
                                                  aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACTCTGAAATCGcc                                         >  1:4834499/1‑50 (MQ=255)
                                                  aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATc                                            >  1:1235976/1‑47 (MQ=255)
                                                  aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGcc                                         >  1:2869033/1‑50 (MQ=255)
                                                  aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa                                      >  1:1711066/1‑53 (MQ=255)
                                                  aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                        <  1:3993389/51‑1 (MQ=255)
                                                  aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                        <  1:638277/51‑1 (MQ=255)
                                                  aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                        <  1:5438222/51‑1 (MQ=255)
                                                  aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                        <  1:3783935/51‑1 (MQ=255)
                                                  aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                        >  1:5662090/1‑51 (MQ=255)
                                                  aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggct                     <  1:2914552/70‑1 (MQ=255)
                                                    aaaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                   >  1:2249693/1‑70 (MQ=255)
                                                    aaaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                   >  1:4894099/1‑70 (MQ=255)
                                                     aaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                   <  1:4086787/69‑1 (MQ=255)
                                                      aGGGACCGAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAc                <  1:572442/71‑1 (MQ=255)
                                                        ggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAcc               >  1:4125620/1‑70 (MQ=255)
                                                        ggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAcc               >  1:1196310/1‑70 (MQ=255)
                                                         gACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGa             <  1:403498/71‑1 (MQ=255)
                                                         gACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGa             <  1:2448142/71‑1 (MQ=255)
                                                          aCCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaa            <  1:4646452/71‑1 (MQ=255)
                                                          aCCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaa            <  1:2999605/71‑1 (MQ=255)
                                                           ccAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa                                      <  1:63704/44‑1 (MQ=255)
                                                           ccAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa           <  1:4162029/71‑1 (MQ=255)
                                                             aaaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa           >  1:1259197/1‑69 (MQ=255)
                                                              aaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctggc                  <  1:2482576/61‑1 (MQ=255)
                                                              aaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCCCCGAAACa         >  1:5118343/1‑70 (MQ=255)
                                                               aaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                     >  1:2172988/1‑41 (MQ=255)
                                                                   gTTGCCGTCACACCAGACTACGCTGAAAACGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAc   >  1:4448240/1‑71 (MQ=255)
                                                                    ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg       <  1:669984/66‑1 (MQ=255)
                                                                    ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg       <  1:3690683/66‑1 (MQ=255)
                                                                    ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc  <  1:1910169/71‑1 (MQ=255)
                                                                     |                                                                     
CGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACC  >  minE/832310‑832448

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: