Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,057 G→T 21.7% V80V (GTG→GTT 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,0570GT21.7% 106.8 / 15.9 46V80V (GTG→GTTnarJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base G (16/20);  new base T (5/5);  total (21/25)
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.39e-01

GCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCT  >  minE/836989‑837124
                                                                    |                                                                   
gCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCa                                                                   >  1:2484570/1‑71 (MQ=255)
     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:1256061/71‑1 (MQ=255)
                  gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcg                                                 <  1:2717555/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:2221230/1‑71 (MQ=255)
                                         cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc                          >  1:2225904/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:1634932/1‑48 (MQ=255)
                                                      cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg                                    >  1:1804337/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:1476639/1‑71 (MQ=255)
                                                                 tgtgCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt  <  1:2562968/71‑1 (MQ=255)
                                                                    |                                                                   
GCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCT  >  minE/836989‑837124

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: