Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,066 A→G 42.9% E83E (GAA→GAG 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,0660AG42.9% 33.9 / 96.3 84E83E (GAA→GAGnarJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base A (20/28);  new base G (18/18);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

CAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTA  >  minE/836998‑837132
                                                                    |                                                                  
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:3275029/1‑70 (MQ=255)
     cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                                              <  1:4655449/70‑1 (MQ=255)
     cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                                              >  1:1954183/1‑70 (MQ=255)
         gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                                              <  1:3267845/66‑1 (MQ=255)
               ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGa                                                            >  1:3726861/1‑62 (MQ=255)
               ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGa                                                            >  1:1313003/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:4219983/71‑3 (MQ=255)
                    ccGCGGGCGCACCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGcaa                                              <  1:2033718/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)
                               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:3472542/1‑71 (MQ=255)
                                cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc                                  >  1:4472078/1‑71 (MQ=255)
                                cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc                                  >  1:2444782/1‑71 (MQ=255)
                                cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc                                  >  1:4174436/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:2408335/71‑1 (MQ=255)
                                             cTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTATGa                     <  1:229261/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:456750/1‑50 (MQ=255)
                                                     aCATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc                                  >  1:4525617/1‑50 (MQ=255)
                                                     aCATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctg                                   >  1:2669587/1‑49 (MQ=255)
                                                     aCATGTTCATGCCGAGTCCAGCGATCGCGGCCAGGCAATGGTGGACctgc                                  >  1:5371076/1‑50 (MQ=255)
                                                     aCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc                                  >  1:848913/1‑50 (MQ=255)
                                                     aCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc                                  >  1:644956/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)
                                                                ccGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTATGAAAAGGTCGGCCTGCAgct   <  1:4813993/69‑3 (MQ=255)
                                                                 cGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTATGAAAAGGTCGGCCTGCAgctg  >  1:1174080/1‑67 (MQ=255)
                                                                    |                                                                  
CAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTA  >  minE/836998‑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: