Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,062 G→C 9.3% 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,0620GC 9.3% 120.8 / 3.7 43G82A (GGC→GCC) narJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base G (22/17);  new base C (2/2);  total (24/19)
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.37e-01

AGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTA  >  minE/836993‑837132
                                                                     |                                                                      
aGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGc                                                                       >  1:3799140/1‑71 (MQ=255)
aGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGc                                                                       >  1:2037232/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)
          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:2800016/1‑70 (MQ=255)
             cGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAc                                                           >  1:2736397/1‑70 (MQ=255)
                aCTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAc                                                           >  1:1340696/1‑67 (MQ=255)
                aCTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATCTACACGGCGAATCCCGCGAc                                                           >  1:528105/1‑67 (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:595993/70‑1 (MQ=255)
                                           aCTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGgcg                             <  1:3956713/70‑1 (MQ=255)
                                              gTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAg                          >  1:1805778/1‑70 (MQ=255)
                                              gTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAg                          >  1:3776751/1‑70 (MQ=255)
                                               ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtgg                                           <  1:242214/52‑1 (MQ=255)
                                               ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtgg                                           <  1:2545857/52‑1 (MQ=255)
                                               ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGAcc                                     <  1:1417426/58‑1 (MQ=255)
                                               ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTa                        >  1:3485711/1‑71 (MQ=255)
                                                     tttGAACATGTGCACGGCGAATCCCGCGACCGCGGc                                                     <  1:362568/36‑1 (MQ=255)
                                                      ttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGa                                              >  1:3951436/1‑42 (MQ=255)
                                                       tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                                                >  1:2153105/1‑39 (MQ=255)
                                                       tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtg                                         >  1:1407091/1‑46 (MQ=255)
                                                       tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcag                 <  1:923048/70‑1 (MQ=255)
                                                          aCATGTGCACGGCGAATCCCGCGACCGCGTCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagca                  >  1:3913568/1‑66 (MQ=255)
                                                          aCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtgg                                           >  1:4076199/1‑41 (MQ=255)
                                                          aCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtgg                                           >  1:3273629/1‑41 (MQ=255)
                                                           cATGTGCACGGCGAATCCCGCGACCGCTGCCAGGc                                                <  1:1895560/35‑1 (MQ=255)
                                                            aTGTGCACGGCGAATCCCGCGACCTCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGc           <  1:3654186/71‑1 (MQ=255)
                                                             tgtgCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt          >  1:3275405/1‑71 (MQ=255)
                                                              gtgCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCtt         >  1:3032360/1‑71 (MQ=255)
                                                                     gCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTa  <  1:2712548/71‑1 (MQ=255)
                                                                     |                                                                      
AGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTA  >  minE/836993‑837132

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: