Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,062 G→C 41.5% G82A (GGC→GCC)  narJ → molybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE837,0620GC41.5% 46.5 / 92.4 82G82A (GGC→GCC) narJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base G (20/28);  new base C (17/17);  total (37/45)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

CAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTA  >  minE/836992‑837132
                                                                      |                                                                      
cAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACgg                                                                        >  1:5328383/1‑71 (MQ=255)
cAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACgg                                                                        >  1:4123467/1‑71 (MQ=255)
  ggCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                                                      <  1:1142602/71‑1 (MQ=255)
   gcgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                                                      >  1:1968832/1‑70 (MQ=255)
   gcgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                                                      >  1:1429854/1‑70 (MQ=255)
      cAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATc                                                                  <  1:3038748/71‑1 (MQ=255)
      cAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATc                                                                  <  1:1273041/71‑1 (MQ=255)
         tACAGCGAACTGTTCGACCGTGGCCGCTCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcg                                                               >  1:5331585/1‑71 (MQ=255)
          aCAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                                              >  1:4351798/1‑71 (MQ=255)
           cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                                              <  1:3180810/70‑1 (MQ=255)
           cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                                              <  1:4655449/70‑1 (MQ=255)
           cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                                              >  1:1954183/1‑70 (MQ=255)
           cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                                              >  1:3275029/1‑70 (MQ=255)
               gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                                              <  1:3267845/66‑1 (MQ=255)
                     ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGa                                                            >  1:1313003/1‑62 (MQ=255)
                     ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGa                                                            >  1:3726861/1‑62 (MQ=255)
                        gACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                                                <  1:347441/71‑1 (MQ=255)
                        gACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAACGCGGCCAGGc                                                <  1:5569837/71‑1 (MQ=255)
                         aCCGCGGGCGCACCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGc                                                <  1:248931/70‑1 (MQ=255)
                          ccGCGGGCGCACCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGcaa                                              <  1:2033718/71‑3 (MQ=255)
                          ccGCGGGCGCACCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGcaa                                              <  1:4219983/71‑3 (MQ=255)
                            gTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg                                            <  1:3691372/71‑1 (MQ=255)
                               gCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggt                                          <  1:1477056/70‑1 (MQ=255)
                                gcgcaCCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggtgg                                        <  1:4929275/70‑1 (MQ=255)
                                gcgcaCCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggtgg                                        <  1:2844976/70‑1 (MQ=255)
                                 cgcgCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                                                      >  1:4555090/1‑40 (MQ=255)
                                     ccACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACct                                    >  1:5264352/1‑70 (MQ=255)
                                      cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtgttgg                                        <  1:5638980/65‑1 (MQ=255)
                                      cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggtgg                                        <  1:4664664/65‑1 (MQ=255)
                                      cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtgg                                           >  1:3892172/1‑62 (MQ=255)
                                      cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc                                  >  1:4472078/1‑71 (MQ=255)
                                      cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc                                  >  1:2444782/1‑71 (MQ=255)
                                      cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc                                  >  1:4174436/1‑71 (MQ=255)
                                      cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc                                  >  1:3472542/1‑71 (MQ=255)
                                      cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtgg                                        <  1:3683966/65‑1 (MQ=255)
                                       aCGTCCCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGcaa                                              >  1:3736215/1‑56 (MQ=255)
                                        cGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgctg                                >  1:3029830/1‑71 (MQ=255)
                                         gTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGc                                                <  1:5464509/54‑1 (MQ=255)
                                          tcgctgctgctgTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGgc                              <  1:1388791/68‑1 (MQ=255)
                                           cgctgctgctgTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGgcg                             >  1:4414028/3‑71 (MQ=255)
                                            gctgctgctgTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtg                                            >  1:3958707/2‑55 (MQ=255)
                                            gctgctgctgTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGgcg                             <  1:1021282/69‑1 (MQ=255)
                                            gctgctgctgTTCGAACATATTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGgcg                             <  1:3803419/69‑1 (MQ=255)
                                            aCTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                                                <  1:900444/51‑1 (MQ=255)
                                            aCTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGgcg                             >  1:4511865/1‑70 (MQ=255)
                                             cTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGgcg                             >  1:955244/1‑69 (MQ=255)
                                              tGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtgg                                        <  1:4333687/57‑1 (MQ=255)
                                              tGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtgg                                        <  1:352813/57‑1 (MQ=255)
                                              tGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAg                          >  1:2359348/1‑71 (MQ=255)
                                               gctGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAg                          >  1:2021444/3‑70 (MQ=255)
                                               gTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                                                <  1:973012/48‑1 (MQ=255)
                                               gTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                                                <  1:2335640/48‑1 (MQ=255)
                                                ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg                                            >  1:1548587/1‑51 (MQ=255)
                                                ctGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc                                  <  1:3437199/60‑1 (MQ=255)
                                                ctGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGt                         <  1:4773093/69‑1 (MQ=255)
                                                ctGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTa                        <  1:1249391/70‑1 (MQ=255)
                                                ctGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTa                        >  1:2197833/2‑71 (MQ=255)
                                                 tGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTgg                               >  1:2438912/1‑63 (MQ=255)
                                                   cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACg                      <  1:241880/70‑1 (MQ=255)
                                                   cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGa                     <  1:1022181/71‑1 (MQ=255)
                                                   cTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGcaa                                              <  1:1369561/46‑3 (MQ=16)
                                                   cTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTATGa                     <  1:229261/71‑1 (MQ=255)
                                                   cTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTATGa                     <  1:2408335/71‑1 (MQ=255)
                                                     gTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGag                    >  1:386382/1‑70 (MQ=255)
                                                     gTTTGAACATGTGCACGGCGAATCCCGCGAC‑‑CGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcag                 >  1:2686838/1‑71 (MQ=255)
                                                      tttGAACTTGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                                                <  1:1025207/41‑1 (MQ=255)
                                                      tttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                                                <  1:3737028/41‑1 (MQ=255)
                                                      tttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctg                                   <  1:2986800/54‑1 (MQ=255)
                                                      tttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGag                    <  1:2469044/69‑1 (MQ=255)
                                                       ttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGa                                              >  1:4370747/1‑42 (MQ=255)
                                                       ttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagca                  <  1:600224/70‑1 (MQ=255)
                                                        tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcag                 <  1:2064394/70‑1 (MQ=255)
                                                        cgAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggtg                                         >  1:5423722/2‑46 (MQ=16)
                                                         gAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGa                     <  1:2838369/65‑1 (MQ=255)
                                                           tcATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc                                  >  1:719290/2‑50 (MQ=255)
                                                           aCATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc                                  >  1:4525617/1‑50 (MQ=255)
                                                           aCATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc                                  >  1:456750/1‑50 (MQ=255)
                                                           aCATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctg                                   >  1:2669587/1‑49 (MQ=255)
                                                           aCATGTTCATGCCGAGTCCAGCGATCGCGGCCAGGCAATGGTGGACctgc                                  >  1:5371076/1‑50 (MQ=255)
                                                           aCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc                                  >  1:644956/1‑50 (MQ=255)
                                                           aCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc                                  >  1:848913/1‑50 (MQ=255)
                                                             aTGTGCACGGCGAATCCCGCGACCCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAgcagca               <  1:4636300/66‑1 (MQ=255)
                                                                tgCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTg        <  1:2243148/71‑1 (MQ=255)
                                                                  cACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGc       <  1:2103049/70‑1 (MQ=255)
                                                                     ggCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGt    >  1:5424827/1‑70 (MQ=255)
                                                                     ggCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGt    >  1:3496941/1‑70 (MQ=255)
                                                                     ggCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACATGCTGGCGCAGTACGAGCAGCACGGCTTGCAGt    >  1:5243615/1‑70 (MQ=255)
                                                                      gCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTa  <  1:4200296/71‑1 (MQ=255)
                                                                      |                                                                      
CAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTA  >  minE/836992‑837132

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: