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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,3760TC 6.2% 278.0 / 5.7 96T278T (ACT→ACCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (49/41);  new base C (4/2);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACA  >  minE/832307‑832441
                                                                     |                                                                 
aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                                                                  >  1:3052503/1‑71 (MQ=255)
  aCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGtc                                                                <  1:3143133/71‑2 (MQ=255)
  aCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGt                                                                 <  1:164303/70‑1 (MQ=255)
   cGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTg                                                               >  1:824326/1‑71 (MQ=255)
    gTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTg                                                               >  1:2942999/1‑70 (MQ=255)
    gTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTg                                                               >  1:1725713/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)
                 gTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                                                                  >  1:482381/1‑54 (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:2691524/71‑1 (MQ=255)
                        ggACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTACTCTg                                          <  1:5512562/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:5139528/1‑51 (MQ=255)
                             cTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                                                         >  1:5423924/1‑51 (MQ=255)
                             cTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                                                         >  1:3486655/1‑51 (MQ=255)
                              tCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                                                                  >  1:4794281/1‑41 (MQ=255)
                              tCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                                                                  >  1:2005069/1‑41 (MQ=255)
                              tCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                                                                  >  1:3764434/1‑41 (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:270543/1‑71 (MQ=255)
                                    cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                              >  1:1454583/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:1294859/1‑71 (MQ=255)
                                             aGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATc                     <  1:1794494/71‑1 (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:95249/47‑1 (MQ=255)
                                                   ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaaa                                       <  1:2893078/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:5438222/51‑1 (MQ=255)
                                                     aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                 <  1:638277/51‑1 (MQ=255)
                                                     aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                 >  1:5662090/1‑51 (MQ=255)
                                                     aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                 <  1:3993389/51‑1 (MQ=255)
                                                     aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                 <  1:3783935/51‑1 (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:2999605/71‑1 (MQ=255)
                                                             aCCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaa     <  1:4646452/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)
                                                                     |                                                                 
AAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACA  >  minE/832307‑832441

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: