Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,086 T→G 19.5% R268R (CGT→CGG narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE836,0860TG19.5% 131.8 / 37.1 77R268R (CGT→CGGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (38/24);  new base G (7/8);  total (45/32)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

TTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATC  >  minE/836018‑836156
                                                                    |                                                                      
ttctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTa                                                                      >  1:5082527/1‑71 (MQ=255)
 tctGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTAt                                                                     >  1:3191624/1‑71 (MQ=255)
 tctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTAt                                                                     >  1:58117/1‑71 (MQ=255)
  ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGtt                                                                       <  1:5263314/68‑1 (MQ=255)
  ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGtt                                                                       <  1:675981/68‑1 (MQ=255)
  ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTa                                                                      <  1:5161340/69‑1 (MQ=255)
     cTATCCGCGTATTTAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCt                                                                   <  1:4448824/69‑1 (MQ=255)
        tCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCtt                                                                  <  1:2319512/67‑1 (MQ=255)
        tCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCtt                                                                  <  1:2269130/67‑1 (MQ=255)
          cgcgTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGt                                                             <  1:1879847/70‑1 (MQ=255)
           gcgTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGc                                                           >  1:2062906/1‑71 (MQ=255)
            cgTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCt                                                          >  1:208666/1‑71 (MQ=255)
               aTTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCt                                                          >  1:2992793/1‑68 (MQ=255)
               aTTGAAGCGGGTCAGACGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgtt                                                       >  1:2292953/1‑71 (MQ=255)
                ttGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                                      >  1:5677282/1‑71 (MQ=255)
                 tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                                      >  1:3607389/1‑70 (MQ=255)
                 tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                                      >  1:3838078/1‑70 (MQ=255)
                 tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                                      <  1:2321273/70‑1 (MQ=255)
                 tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                                      >  1:866778/1‑70 (MQ=255)
                 tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                                      >  1:4447281/1‑70 (MQ=255)
                  gAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGAGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTa                                                    >  1:269440/1‑71 (MQ=255)
                   aaGCGTGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTa                                                    <  1:3936445/70‑1 (MQ=255)
                    aGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACg                                                  <  1:4309061/71‑1 (MQ=255)
                     tccggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCt                                                          <  1:1297308/59‑1 (MQ=255)
                     tccggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCt                                                          <  1:4354187/59‑1 (MQ=255)
                     tccggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTACGa                                                 >  1:276075/4‑71 (MQ=255)
                     gCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                                 >  1:4667558/1‑71 (MQ=255)
                       cggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCtgc                                                        >  1:5670470/2‑61 (MQ=255)
                       cggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCtgc                                                        >  1:1297300/2‑61 (MQ=255)
                       cggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCtgc                                                        >  1:3915422/2‑61 (MQ=255)
                       cggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTACGac                                                <  1:812112/69‑2 (MQ=255)
                         gTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTAc                                                   <  1:2662269/65‑1 (MQ=255)
                          tCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTACGACGCCg                                            <  1:5551407/71‑1 (MQ=255)
                          tCAACCGACCGTGTGCTCAGAAACCTGCATGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTACGACGcc                                             >  1:2781159/1‑70 (MQ=255)
                           cAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGa                                           >  1:1801949/1‑71 (MQ=255)
                            aGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGc                                                               <  1:3612705/50‑1 (MQ=255)
                            aGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGc                                                               <  1:2213539/50‑1 (MQ=255)
                            aGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGAc                                          >  1:4856313/1‑71 (MQ=255)
                            aGCAGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGAc                                          >  1:4834711/1‑71 (MQ=255)
                             gCCGACCGTGTGATCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGcc                                             <  1:1669589/67‑1 (MQ=255)
                              ccGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTAc                                                   <  1:1369153/60‑1 (MQ=255)
                              ccGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTAt                                                                     >  1:3189075/1‑42 (MQ=25)
                               cGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGc                                        <  1:3257921/70‑1 (MQ=255)
                               cGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGc                                        >  1:125299/1‑70 (MQ=255)
                                   cGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                                 >  1:1920537/1‑57 (MQ=255)
                                      gtgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCt                                                                   <  1:5409495/36‑1 (MQ=255)
                                      gtgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCt                                                                   <  1:2584512/36‑1 (MQ=255)
                                      gtgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAAc                                >  1:1712571/1‑71 (MQ=255)
                                      gtgCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTAc                                                   <  1:1753677/52‑1 (MQ=255)
                                      gtgCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTAc                                                   <  1:4624859/52‑1 (MQ=255)
                                      gtgCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTAc                                                   <  1:1364331/52‑1 (MQ=255)
                                       tgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                                 >  1:5142558/1‑53 (MQ=255)
                                       tgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                                 >  1:1093989/1‑53 (MQ=255)
                                       tgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                                 >  1:1306591/1‑53 (MQ=255)
                                       tgCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTAc                                                   <  1:1447620/51‑1 (MQ=255)
                                         cTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTg                             <  1:5285365/71‑1 (MQ=255)
                                          tCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGt                                                             <  1:100427/38‑1 (MQ=255)
                                            aGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTa                                                    >  1:5124531/1‑45 (MQ=255)
                                            aGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGc                            <  1:447159/69‑1 (MQ=255)
                                               aaCCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGcagc                         >  1:1542849/1‑69 (MQ=255)
                                                aCCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGcagccag                      <  1:1057236/71‑1 (MQ=255)
                                                aCCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGcagccag                      <  1:4551761/71‑1 (MQ=255)
                                                  cTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGCCGCTATTGAACGTGCAGCCAGCa                    >  1:3687299/1‑71 (MQ=255)
                                                  cTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCa                    >  1:5308503/1‑71 (MQ=255)
                                                   tgtgTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCAc                   >  1:516212/1‑71 (MQ=255)
                                                    gtgtCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGcagc                         >  1:1822240/1‑64 (MQ=255)
                                                     tgtCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCAcc                  >  1:1109144/1‑70 (MQ=255)
                                                          gtcgTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGc                                        <  1:4479607/43‑1 (MQ=255)
                                                          gtcgTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGc                                        <  1:4986321/43‑1 (MQ=255)
                                                            cgTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaac            >  1:1731161/1‑69 (MQ=255)
                                                            cgTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaac            >  1:4635675/1‑69 (MQ=255)
                                                             gTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaac            >  1:4330242/1‑68 (MQ=255)
                                                              tATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaacgaga        <  1:5545138/71‑1 (MQ=255)
                                                               aTCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAAc                                >  1:4102071/1‑46 (MQ=255)
                                                               aTCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaacgaga        >  1:3000351/1‑70 (MQ=255)
                                                                tCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaacgaga        >  1:3121/1‑69 (MQ=255)
                                                                tCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaac            >  1:5270669/1‑65 (MQ=255)
                                                                tCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaac            >  1:84392/1‑65 (MQ=255)
                                                                  cGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgag               >  1:2697642/1‑60 (MQ=255)
                                                                    ttATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATc  >  1:3305631/1‑71 (MQ=255)
                                                                    ttATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATc  >  1:1709170/1‑71 (MQ=255)
                                                                    ttATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATc  >  1:1152525/1‑71 (MQ=255)
                                                                    |                                                                      
TTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATC  >  minE/836018‑836156

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: