Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 837,036 A→G 24.9% S73S (TCA→TCG 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,0360AG24.9% 97.5 / 33.2 62S73S (TCA→TCGnarJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base A (19/27);  new base G (7/8);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.74e-01

ACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCG  >  Exported/836968‑837105
                                                                    |                                                                     
acgacgATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCAc                                                                      <  1:1792229/70‑1 (MQ=255)
           tCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGttt                                                           >  1:1243331/1‑70 (MQ=255)
           tCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg                                                          <  1:3830050/71‑1 (MQ=255)
           tCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg                                                          <  1:2925923/71‑1 (MQ=255)
             cGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa                                                        <  1:4486640/71‑1 (MQ=255)
             cGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa                                                        <  1:2634176/71‑1 (MQ=255)
              gCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa                                                        <  1:3021745/70‑1 (MQ=255)
               cTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTg                                                                 <  1:3730991/60‑1 (MQ=255)
               cTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTg                                                                 <  1:3115401/60‑1 (MQ=255)
               cTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa                                                        <  1:679740/69‑1 (MQ=255)
               cTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACa                                                      <  1:5285151/71‑1 (MQ=255)
                  gATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtgt                                                   >  1:4743326/1‑71 (MQ=255)
                      cccAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCa                                                <  1:5129755/70‑1 (MQ=255)
                      cccAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCa                                                <  1:2142867/70‑1 (MQ=255)
                       ccAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACg                                              >  1:4304606/1‑71 (MQ=255)
                       ccAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACg                                              <  1:3982208/71‑1 (MQ=255)
                        cAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCAc                                               >  1:528943/1‑69 (MQ=255)
                        cAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCAc                                               >  1:1774486/1‑69 (MQ=255)
                        cAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACgg                                             >  1:4123467/1‑71 (MQ=255)
                        cAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACgg                                             >  1:5328383/1‑71 (MQ=255)
                          ggCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                           <  1:1142602/71‑1 (MQ=255)
                           gcgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                           >  1:1429854/1‑70 (MQ=255)
                           gcgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                           >  1:1968832/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)
                                   cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCAc                                                                      >  1:487399/1‑35 (MQ=255)
                                   cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACt                                                                     >  1:5232874/1‑36 (MQ=255)
                                   cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                   <  1:3180810/70‑1 (MQ=255)
                                   cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                   <  1:4655449/70‑1 (MQ=255)
                                   cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                   >  1:3275029/1‑70 (MQ=255)
                                   cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                   >  1:1954183/1‑70 (MQ=255)
                                     gCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg                                                          <  1:5411941/45‑1 (MQ=255)
                                       gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                   <  1:3267845/66‑1 (MQ=255)
                                        aaCTGTTCGACCGTGGCCGCGCCACGTCACTGTTg                                                                 <  1:4796508/35‑1 (MQ=255)
                                            gTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGt                                                             <  1:2861544/35‑1 (MQ=255)
                                             ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGa                                                         >  1:2390564/1‑38 (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)
                                                         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: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)
                                                                    aCTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                     <  1:900444/51‑1 (MQ=255)
                                                                    aCTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGgcg  >  1:4511865/1‑70 (MQ=255)
                                                                    |                                                                     
ACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCG  >  Exported/836968‑837105

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: