Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,388 A→C 7.5% T282T (ACA→ACC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,3880AC 7.5% 203.9 / 8.7 80T282T (ACA→ACCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (34/40);  new base C (4/2);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAA  >  minE/832321‑832458
                                                                   |                                                                      
cgcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCa                                                                     <  1:568557/71‑1 (MQ=255)
 gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCacac                                                                       >  1:4665990/1‑68 (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)
          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)
                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:2640918/1‑70 (MQ=255)
                  aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGc                                                    >  1:2624974/1‑70 (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:2893078/47‑1 (MQ=255)
                                     ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaaa                                                        <  1:95249/47‑1 (MQ=255)
                                     ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt                                            <  1:1713914/59‑1 (MQ=255)
                                     ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt                                            <  1:617241/59‑1 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACTCTGAAATCGcc                                                   >  1:4834499/1‑50 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACTCTGAAATCGcc                                                   >  1:4828489/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:638277/51‑1 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                                  <  1:3783935/51‑1 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                                  >  1:5662090/1‑51 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                                  <  1:5438222/51‑1 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                                  <  1:3993389/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:1196310/1‑70 (MQ=255)
                                             ggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAcc                         >  1:4125620/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)
                                                            ccGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGAt         <  1:3196267/71‑1 (MQ=255)
                                                               tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACa                   <  1:894360/58‑1 (MQ=255)
                                                               tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACa                   <  1:3756818/58‑1 (MQ=255)
                                                               tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGa          <  1:3866823/67‑1 (MQ=255)
                                                                cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCtg                                    >  1:2355399/1‑40 (MQ=255)
                                                                cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCtg                                    >  1:2958692/1‑40 (MQ=255)
                                                                cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGc               >  1:2585835/1‑61 (MQ=255)
                                                                cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg     >  1:4275793/1‑71 (MQ=255)
                                                                cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg     >  1:2149322/1‑71 (MQ=255)
                                                                cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg     >  1:2155976/1‑71 (MQ=255)
                                                                 acacCAGACTACGCTGAAATCGCCAAACTGTGCGAt                                       >  1:4251374/1‑36 (MQ=255)
                                                                 acacCAGACTACGCTGAAATCGCCAAACTGTGCGAt                                       >  1:3032130/1‑36 (MQ=255)
                                                                 acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                     <  1:2837562/54‑1 (MQ=255)
                                                                 acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                     <  1:2830678/54‑1 (MQ=255)
                                                                 acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                 >  1:592428/1‑58 (MQ=255)
                                                                 acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGc    <  1:5383525/71‑1 (MQ=255)
                                                                   acCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACa                   <  1:2392515/54‑1 (MQ=255)
                                                                   acCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCaa  <  1:3964059/71‑1 (MQ=255)
                                                                   |                                                                      
CGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAA  >  minE/832321‑832458

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: