Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,381 T→C 24.9% R33R (CGT→CGC narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,3810TC24.9% 152.8 / 53.9 91R33R (CGT→CGCnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base C (11/11);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

AATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTG  >  minE/835316‑835449
                                                                 |                                                                    
aacctcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGaa                                                                   <  1:2160995/66‑1 (MQ=255)
aacctcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAg                                                                  >  1:2376194/4‑70 (MQ=255)
aacctcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAg                                                                  >  1:1130916/4‑70 (MQ=255)
 aTCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGc                                                                <  1:70033/71‑1 (MQ=255)
 aTCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGc                                                                <  1:3043866/71‑1 (MQ=255)
 aTCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGc                                                                <  1:2687975/71‑1 (MQ=255)
  cctcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGca                                                               <  1:3542248/70‑2 (MQ=255)
    tcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGt                                                                      >  1:3842024/1‑62 (MQ=255)
     cGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg                                                             >  1:1826224/1‑70 (MQ=255)
      gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGa                                                           <  1:1793399/71‑1 (MQ=255)
       aTAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGaa                                                          <  1:3952560/71‑1 (MQ=255)
         aaGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                        >  1:2162018/1‑71 (MQ=255)
         aaGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                        >  1:2958253/1‑71 (MQ=255)
          aGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                        <  1:931629/70‑1 (MQ=255)
          aGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                        <  1:2416959/70‑1 (MQ=255)
          aGTGCATCGGCTGCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATg                                                      >  1:2513863/1‑71 (MQ=255)
           gTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAatat                                                       >  1:1832505/1‑70 (MQ=255)
            tGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGc                                                     <  1:2742999/71‑1 (MQ=255)
            tGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGc                                                     <  1:2472128/71‑1 (MQ=255)
              cATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGt                                                   >  1:3270565/1‑71 (MQ=255)
              cATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGt                                                   >  1:2047527/1‑71 (MQ=255)
               aTCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTg                                                  <  1:69819/71‑1 (MQ=255)
                  ggCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGtt                                               >  1:430549/1‑71 (MQ=255)
                  ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGt                                                <  1:2436538/70‑1 (MQ=255)
                   gCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGttt                                              <  1:109/71‑2 (MQ=255)
                   gCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGtt                                               >  1:2297446/1‑70 (MQ=255)
                   gCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc                                              <  1:2074764/71‑1 (MQ=255)
                    cTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGttt                                              <  1:3108489/70‑2 (MQ=255)
                     tGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaa                                            >  1:2398714/1‑71 (MQ=255)
                     tGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaa                                            >  1:2051205/1‑71 (MQ=255)
                     tGCCACACCTGCTCGATCACCTGCAAGAACGTCTGGACCAGCCGCGAAGGCATGGAGTACGCCTGGTTcaa                                            >  1:709105/1‑71 (MQ=255)
                       ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                           >  1:2892066/1‑70 (MQ=255)
                       ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                           >  1:1854520/1‑70 (MQ=255)
                       ccACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaa                                            <  1:3337043/69‑1 (MQ=255)
                          tacCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACg                                       <  1:2538245/70‑1 (MQ=255)
                          cacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATg                                                      >  1:3634736/1‑56 (MQ=255)
                          cacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg                                       <  1:2439287/71‑1 (MQ=255)
                          cacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg                                       <  1:2316976/71‑1 (MQ=255)
                           acCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGt                                      >  1:2166895/1‑71 (MQ=255)
                             cTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATg                                                             >  1:1677135/1‑46 (MQ=255)
                             cTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAat                                                         >  1:3027312/1‑50 (MQ=255)
                             cTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaaca                                          >  1:1893268/1‑65 (MQ=255)
                             cTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTgg                                    <  1:1327911/71‑1 (MQ=255)
                             cTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                     >  1:3401816/1‑70 (MQ=255)
                              tGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATg                                                      >  1:1440941/1‑52 (MQ=255)
                                ttgaGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaa                                  <  1:2134511/67‑1 (MQ=255)
                                ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaaa                                 >  1:2118523/1‑71 (MQ=255)
                                ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaa                                  <  1:243508/70‑1 (MQ=255)
                                ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGa                                   <  1:2919966/69‑1 (MQ=255)
                                ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTca                                             <  1:3307356/59‑1 (MQ=255)
                                ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                   <  1:3215110/69‑1 (MQ=255)
                                 tCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaac                                           >  1:699739/1‑60 (MQ=255)
                                 tCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaac                                           >  1:1884875/1‑60 (MQ=255)
                                  cAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGaa                                                          >  1:1486749/1‑44 (MQ=255)
                                  cAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAc                                <  1:3762977/70‑1 (MQ=255)
                                   aGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCa                              <  1:1729499/71‑1 (MQ=255)
                                    gTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAcc                               >  1:2873644/1‑69 (MQ=255)
                                    gTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCaa                             <  1:3256987/71‑1 (MQ=255)
                                     ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTca                                             <  1:1815241/54‑1 (MQ=255)
                                     ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAcc                               <  1:2143466/68‑1 (MQ=255)
                                     ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                            <  1:942218/71‑1 (MQ=255)
                                     tgaCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACg                                       <  1:3412602/58‑1 (MQ=255)
                                       aCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGtc                                     >  1:2752404/1‑59 (MQ=255)
                                       aCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaaaccaaacc                          >  1:690504/1‑71 (MQ=255)
                                       aCCTGTAAAAACGTCTGGACCGGGCGAGAAGGCATGGAGTACGCATGGTTTAACAACGtc                                     >  1:3463912/1‑59 (MQ=255)
                                       aCCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCaa                             <  1:1861070/68‑1 (MQ=255)
                                       aCCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGcc                          >  1:3361058/1‑71 (MQ=255)
                                        ccTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGcc                          <  1:2274587/70‑1 (MQ=255)
                                         cTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGtc                                     <  1:2574569/58‑2 (MQ=255)
                                          tGTAAAAACGTCTGGACCTGGCGCGAAGGCATGGAgtac                                                       <  1:1932449/39‑2 (MQ=16)
                                           ttAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCggg                       >  1:3332404/2‑70 (MQ=255)
                                            tAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTaacaac                                        <  1:1343900/52‑1 (MQ=255)
                                            tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa                                         >  1:2718493/1‑51 (MQ=255)
                                            tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGcc                     <  1:188282/71‑1 (MQ=255)
                                            tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGcc                     <  1:3938006/71‑1 (MQ=255)
                                            tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGc                      >  1:1180552/1‑70 (MQ=255)
                                              aaaaCGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCAAACCGGGcatt                   >  1:1561972/1‑68 (MQ=255)
                                              aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCa                    <  1:3446346/70‑1 (MQ=255)
                                                 aCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaac                                           >  1:1052481/1‑44 (MQ=255)
                                                 aCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaac                                           >  1:1236055/1‑44 (MQ=255)
                                                 aCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGc                >  1:1035147/1‑71 (MQ=255)
                                                 aCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGc                >  1:2327977/1‑71 (MQ=255)
                                                 aCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAAGGc                <  1:1573352/71‑1 (MQ=255)
                                                  cGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCt               <  1:1203280/71‑1 (MQ=255)
                                                     cTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCa                    >  1:2080787/1‑63 (MQ=255)
                                                        gACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc           <  1:2819519/69‑1 (MQ=255)
                                                             gCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac    >  1:3545396/1‑71 (MQ=255)
                                                              ccGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCa                              <  1:2673428/44‑1 (MQ=255)
                                                              ccGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgact   <  1:1795505/71‑1 (MQ=255)
                                                              ccGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgact   >  1:3544092/1‑71 (MQ=255)
                                                               cGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGactgactg  >  1:546169/1‑71 (MQ=255)
                                                               cGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGactgactg  >  1:3840697/1‑71 (MQ=255)
                                                                 tGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCtt              >  1:3829742/1‑57 (MQ=255)
                                                                 |                                                                    
AATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTG  >  minE/835316‑835449

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: