Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,060 C→T 20.0% 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,0600CT20.0% 103.5 / 19.8 50H81H (CAC→CATnarJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base C (16/24);  new base T (5/5);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

GGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCA  >  minE/836994‑837128
                                                                  |                                                                    
ggCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGc                                                                   >  1:850909/1‑70 (MQ=255)
ggCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                                                  <  1:953802/71‑1 (MQ=255)
ggCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                                                  <  1:3080495/71‑1 (MQ=255)
 gcgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGa                                                                 <  1:2542945/71‑1 (MQ=255)
  cgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGaa                                                                <  1:3502846/71‑1 (MQ=255)
        aCAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                                          >  1:859730/1‑71 (MQ=255)
          aGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGaa                                                                <  1:3395378/63‑1 (MQ=255)
            agAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAcc                                                      >  1:3444475/2‑71 (MQ=255)
             gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAcc                                                      >  1:1583246/1‑70 (MQ=255)
             gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcg                                                     <  1:2717555/71‑1 (MQ=255)
             gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcg                                                     <  1:1256061/71‑1 (MQ=255)
             gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcg                                                     >  1:2089670/1‑71 (MQ=255)
                   ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGcc                                                <  1:1099489/70‑1 (MQ=255)
                   ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGc                                                 >  1:1050852/1‑69 (MQ=255)
                   ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCa                                               <  1:1531943/71‑1 (MQ=255)
                     cGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGa                                                                 <  1:3029579/51‑1 (MQ=255)
                       aCCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                                            >  1:2398564/1‑70 (MQ=255)
                          gTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg                                        <  1:2108400/71‑1 (MQ=255)
                            gggcgcaCCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggt                                      <  1:874864/68‑1 (MQ=255)
                               cgcaCCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGa                                   <  1:2391791/71‑1 (MQ=255)
                                  accACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACct                                >  1:3603164/2‑71 (MQ=255)
                                  accACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACct                                >  1:1722484/2‑71 (MQ=255)
                                   ccACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGc                                            <  1:3107988/58‑1 (MQ=255)
                                    cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc                              >  1:1785592/1‑71 (MQ=255)
                                    cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc                              >  1:2225904/1‑71 (MQ=255)
                                    cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc                              >  1:2221230/1‑71 (MQ=255)
                                    cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc                              <  1:1465211/71‑1 (MQ=255)
                                      cGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtg                                     <  1:1637635/62‑1 (MQ=255)
                                          gctgctgctgTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggt                                      >  1:1731759/2‑57 (MQ=255)
                                          aCTGTTGCTGTTTGAACATGTGCACGGCTAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgctg                            <  1:2037465/67‑1 (MQ=255)
                                            tGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggt                                      <  1:2800032/55‑1 (MQ=255)
                                            tGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGAcc                                 >  1:1267926/1‑60 (MQ=255)
                                              ctGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTa                    <  1:29106/70‑1 (MQ=255)
                                                 cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg                                        >  1:1804337/1‑48 (MQ=255)
                                                 cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg                                        >  1:1634932/1‑48 (MQ=255)
                                                 cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGa                 <  1:3321884/71‑1 (MQ=255)
                                                 cTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggtgg                                    <  1:537710/52‑1 (MQ=255)
                                                   gTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGag                >  1:3108591/1‑70 (MQ=255)
                                                    tttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGag                <  1:3216655/69‑1 (MQ=255)
                                                    ttcgAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTATGa                 >  1:3425150/4‑68 (MQ=255)
                                                     ttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGag                <  1:2047565/68‑1 (MQ=255)
                                                      tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcag             >  1:2668388/1‑70 (MQ=255)
                                                      tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGAGCAGTACAagcag             >  1:326620/1‑70 (MQ=255)
                                                       gAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcagc            <  1:2181591/70‑1 (MQ=255)
                                                          cATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACgg        <  1:930769/71‑1 (MQ=255)
                                                            tgtgCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt      <  1:2562968/71‑1 (MQ=255)
                                                            tgtgCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt      >  1:1476639/1‑71 (MQ=255)
                                                                 aCGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGa                                   <  1:277858/37‑1 (MQ=255)
                                                                 aCGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGa                                   <  1:1546293/37‑1 (MQ=255)
                                                                 aCGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCa  <  1:2165574/70‑1 (MQ=255)
                                                                  |                                                                    
GGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCA  >  minE/836994‑837128

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: