Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 837,106 C→G 15.1% Q97E (CAG→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
*Exported837,1060CG15.1% 125.1 / 14.2 53Q97E (CAG→GAG) narJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base C (21/24);  new base G (3/5);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAG  >  Exported/837038‑837174
                                                                    |                                                                    
tGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAg                                                                    >  1:2359348/1‑71 (MQ=255)
 gctGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAg                                                                    >  1:2021444/3‑70 (MQ=255)
  ctGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGt                                                                   <  1:4773093/69‑1 (MQ=255)
  ctGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTa                                                                  <  1:1249391/70‑1 (MQ=255)
  ctGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTa                                                                  >  1:2197833/2‑71 (MQ=255)
     cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACg                                                                <  1:241880/70‑1 (MQ=255)
     cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGa                                                               <  1:1022181/71‑1 (MQ=255)
     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)
        tttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGag                                                              <  1:2469044/69‑1 (MQ=255)
         ttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagca                                                            <  1:600224/70‑1 (MQ=255)
          tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcag                                                           <  1:2064394/70‑1 (MQ=255)
           gAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGa                                                               <  1:2838369/65‑1 (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)
                              cccGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGc                                      >  1:5614802/1‑71 (MQ=255)
                                cgcgACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCcg                                    <  1:219901/71‑1 (MQ=255)
                                cgcgACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCcg                                    >  1:1229473/1‑71 (MQ=255)
                                 gcgACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCcgc                                   >  1:4320596/1‑71 (MQ=255)
                                 gcgACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCcg                                    >  1:2388511/1‑70 (MQ=255)
                                  cgACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCcgcg                                  <  1:1779828/71‑1 (MQ=255)
                                   gACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGa                                 <  1:2819522/71‑1 (MQ=255)
                                   gACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGa                                 <  1:284705/71‑1 (MQ=255)
                                     ccGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAAt                               <  1:1248093/71‑1 (MQ=255)
                                      cgcgGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAAtt                              >  1:2738719/1‑71 (MQ=255)
                                         ggCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGcc                           <  1:349926/71‑1 (MQ=255)
                                              ggCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGAcc                      <  1:3009003/71‑1 (MQ=255)
                                               gCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCcgcg                                  <  1:3701955/58‑1 (MQ=255)
                                               gCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCa                     >  1:1445113/1‑71 (MQ=255)
                                               gCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCa                     <  1:3257223/71‑1 (MQ=255)
                                                cGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCTAATTGCCGGACCAt                    <  1:2088455/71‑1 (MQ=255)
                                                 gATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATc                   >  1:4579692/1‑71 (MQ=255)
                                                   tggtggACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTaaa                                          >  1:3874702/1‑46 (MQ=255)
                                                        gACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACa                                        >  1:942909/1‑43 (MQ=255)
                                                        gACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACa                                        >  1:2935956/1‑43 (MQ=255)
                                                           ctgctgGCGCAGTACGAGCAGCACGGCTTGCAGTTaa                                           <  1:3325432/37‑1 (MQ=255)
                                                           ctgctgGCGCAGTACGAGCAGCACGGCTTGCAGTTaa                                           <  1:1507171/37‑1 (MQ=255)
                                                           ctgctgGCGCAGTACGAGCAGCACGGCTTGCAGTTAAAc                                         <  1:4240336/39‑1 (MQ=255)
                                                           ctgctgGCGCAGTACGAGCAGCACGGCTTGCAGTTAAAc                                         <  1:3137339/39‑1 (MQ=255)
                                                             gctgGCGCAGTACGAGCAGCACGGCTTGCAGTTAAAc                                         >  1:1851493/1‑37 (MQ=255)
                                                              ctgGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTg                 >  1:4333393/1‑60 (MQ=255)
                                                              ctgGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTg                 >  1:3206413/1‑60 (MQ=255)
                                                              ctgGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTg                 >  1:127286/1‑60 (MQ=255)
                                                                ggCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTg     >  1:225439/1‑70 (MQ=255)
                                                                 gcgcAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTgg    <  1:4604577/70‑1 (MQ=255)
                                                                 gcgcAGTACGAGCAGCACAGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTg     <  1:442613/69‑1 (MQ=255)
                                                                  cgcAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAg  >  1:4870630/1‑71 (MQ=255)
                                                                  cgcAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAg  >  1:5282529/1‑71 (MQ=255)
                                                                  cgcAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAg  >  1:821237/1‑71 (MQ=255)
                                                                    |                                                                    
TGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAG  >  Exported/837038‑837174

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