Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,060 C→T 9.3% H81H (CAC→CAT 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,0600CT 9.3% 126.6 / 3.8 43H81H (CAC→CATnarJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base C (22/17);  new base T (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 = 9.27e-01

AGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTT  >  minE/836993‑837125
                                                                   |                                                                 
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)
                                                                   |                                                                 
AGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTT  >  minE/836993‑837125

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: