Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,036 A→G 10.5% 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
*minE837,0360AG10.5% 95.0 / 4.7 38S73S (TCA→TCGnarJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base A (19/15);  new base G (2/2);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

TAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCG  >  minE/836966‑837105
                                                                      |                                                                     
tAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCa                                                                       >  1:623575/1‑71 (MQ=255)
           gATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGt                                                             <  1:2362405/70‑1 (MQ=255)
             tCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGttt                                                           >  1:4065062/1‑70 (MQ=255)
             tCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGt                                                             >  1:4127202/1‑68 (MQ=255)
             tCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg                                                          <  1:80472/71‑1 (MQ=255)
              ccGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg                                                          >  1:590696/1‑70 (MQ=255)
               cGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa                                                        >  1:149827/1‑71 (MQ=255)
                gCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAAc                                                       <  1:2005756/71‑1 (MQ=255)
                 cTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACa                                                      >  1:3852910/1‑71 (MQ=255)
                    gATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt                                                     <  1:2311282/69‑1 (MQ=255)
                     aTGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtgtg                                                  <  1:422146/71‑1 (MQ=255)
                     aTGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtgtg                                                  <  1:1875212/71‑1 (MQ=255)
                       gCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCAc                                                                      >  1:950602/1‑49 (MQ=255)
                       gCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCAc                                                                      >  1:1495074/1‑49 (MQ=255)
                           aGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGc                                            >  1:2037232/1‑71 (MQ=255)
                           aGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGc                                            >  1:3799140/1‑71 (MQ=255)
                             gcgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                           <  1:3288558/70‑1 (MQ=255)
                             gcgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                           >  1:954371/1‑70 (MQ=255)
                                   tACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                           <  1:519187/64‑1 (MQ=255)
                                   tACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATccc                                     >  1:3704445/1‑70 (MQ=255)
                                   tACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcg                                    >  1:1823996/1‑71 (MQ=255)
                                     cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtg                                                    >  1:2546614/1‑53 (MQ=255)
                                     cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                   >  1:1744669/1‑70 (MQ=255)
                                     cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                   <  1:1710408/70‑1 (MQ=255)
                                        cGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                   <  1:3057701/67‑1 (MQ=255)
                                        cGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAcc                               >  1:160268/1‑71 (MQ=255)
                                        cGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAc                                >  1:2736397/1‑70 (MQ=255)
                                        cGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAc                                >  1:2800016/1‑70 (MQ=255)
                                           aCTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAc                                >  1:1340696/1‑67 (MQ=255)
                                           aCTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATCTACACGGCGAATCCCGCGAc                                >  1:528105/1‑67 (MQ=255)
                                               ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAAc                                                       <  1:1029069/40‑1 (MQ=255)
                                                     cGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgcg                                  <  1:404767/55‑1 (MQ=255)
                                                          ccGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtgg             <  1:3914957/71‑1 (MQ=255)
                                                           cgcaCCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGa            <  1:23359/71‑1 (MQ=255)
                                                              gcCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACct         <  1:2252290/71‑1 (MQ=255)
                                                              accACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACct         >  1:857270/2‑71 (MQ=255)
                                                                 aCGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc       >  1:46547/1‑70 (MQ=255)
                                                                  cGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGCACctgctg     <  1:3545037/71‑1 (MQ=255)
                                                                   gTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTgg    <  1:3787014/71‑1 (MQ=255)
                                                                      aCTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGgcg  <  1:3956713/70‑1 (MQ=255)
                                                                      aCTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGgcg  <  1:595993/70‑1 (MQ=255)
                                                                      |                                                                     
TAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCG  >  minE/836966‑837105

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: