Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,446 G→A 6.4% Q388Q (CAG→CAA 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,4460GA 6.4% 202.1 / 3.2 62Q388Q (CAG→CAAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base A (2/2);  total (29/33)
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.80e-01

CAGCAGACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGC  >  minE/836376‑836514
                                                                      |                                                                    
cagcagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAg                                                                      >  1:3973951/1‑71 (MQ=255)
cagcagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAg                                                                      <  1:460520/71‑1 (MQ=255)
 agcagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGt                                                                     <  1:3523021/71‑1 (MQ=255)
 agcagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGt                                                                     <  1:3938332/71‑1 (MQ=255)
   cagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTa                                                                    <  1:1149689/70‑1 (MQ=255)
       cgcgGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTgg                                                               >  1:2122209/1‑71 (MQ=255)
        gcgGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGc                                                              >  1:3698168/1‑71 (MQ=255)
          gggTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCg                                                             >  1:3883790/1‑70 (MQ=255)
          gggTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGa                                                            <  1:2927223/71‑1 (MQ=255)
                    ggTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctgctg                                                    <  1:1047054/69‑1 (MQ=255)
                    ggTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctgctg                                                    <  1:1042115/69‑1 (MQ=255)
                      tAGCAACGGCATTCTGCCAGACGTAGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATc                                                         >  1:2414971/1‑62 (MQ=255)
                        gCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctg                                                       <  1:3048575/62‑1 (MQ=255)
                          aaCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAAt                                                          >  1:3515285/1‑57 (MQ=255)
                          aaCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAAt                                                          >  1:3933056/1‑57 (MQ=255)
                            cgtcATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTg                                          <  1:3695159/68‑1 (MQ=255)
                             ggCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGa                                         <  1:3049400/71‑1 (MQ=255)
                             ggCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGa                                         <  1:2825019/71‑1 (MQ=255)
                               cATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATa                                       >  1:324108/1‑71 (MQ=255)
                                 ttCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGt                                                                     <  1:27138/39‑1 (MQ=255)
                                   cTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTg                                          >  1:2542316/1‑64 (MQ=255)
                                   cTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCa                                    >  1:1357481/1‑70 (MQ=255)
                                    tGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATAc                                      >  1:1358649/1‑67 (MQ=255)
                                    tGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCaaa                                  >  1:1501839/1‑71 (MQ=255)
                                      ccAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAAc                                 <  1:763241/70‑1 (MQ=255)
                                        aGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCt                                                                 >  1:1184506/1‑36 (MQ=255)
                                        aGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                           >  1:2182173/1‑58 (MQ=255)
                                        aGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                              <  1:3690589/71‑1 (MQ=255)
                                          aCGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAAtcttct                                                     <  1:1505454/46‑1 (MQ=255)
                                             tCGAAAGCCTGTGTATTCCGGTGCAATATCTCGCCaa                                                           <  1:759919/37‑1 (MQ=16)
                                             tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                           <  1:468029/37‑1 (MQ=21)
                                             tCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                           <  1:900946/53‑1 (MQ=255)
                                               gAAAGCCTGCGTATTCCGGTGCAATATCTCGCAAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTgcg                       >  1:305060/1‑71 (MQ=255)
                                               gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGaa                                                           >  1:1727636/1‑35 (MQ=255)
                                               gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctgctg                                                    >  1:2621622/1‑42 (MQ=255)
                                               gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctg                                                       <  1:573089/39‑1 (MQ=255)
                                               gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccgcc                                              <  1:599095/48‑1 (MQ=255)
                                               gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccgcc                                              <  1:594163/48‑1 (MQ=255)
                                               gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATAc                                      >  1:2178356/1‑56 (MQ=255)
                                               gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcg                       >  1:1760449/1‑71 (MQ=255)
                                                 aaGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccgcc                                              >  1:3962629/1‑46 (MQ=255)
                                                 aaGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgcg                     >  1:70665/1‑71 (MQ=255)
                                                  aGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgcgc                    <  1:2966477/71‑1 (MQ=255)
                                                  aGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgcg                     >  1:3622679/1‑70 (MQ=255)
                                                  aGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGATGATACCAAACCGGTACTgcgcg                     >  1:3739535/1‑70 (MQ=255)
                                                    ccTGCGTATTCCGGTACAGTATCTGGCGAATCTGGTGACCGCCGGTGATACCAAACCGGTACTGCGCGCAc                  <  1:60223/71‑1 (MQ=255)
                                                     cTGCGTATTACGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACt                 <  1:2159068/71‑1 (MQ=255)
                                                       gCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTg                >  1:1951945/1‑70 (MQ=255)
                                                        cGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCaaa                                  >  1:3933348/1‑51 (MQ=255)
                                                        cGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                              <  1:1218914/55‑1 (MQ=255)
                                                        cGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                              <  1:3096074/55‑1 (MQ=255)
                                                            ttCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgc                      >  1:840668/1‑59 (MQ=255)
                                                            ttCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACg           >  1:1156982/1‑70 (MQ=255)
                                                            ttCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACg           >  1:578007/1‑70 (MQ=255)
                                                              ccGGTACAGTATCTGGCGCATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTa         >  1:1809069/1‑70 (MQ=255)
                                                              ccGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATAc                                      >  1:175459/1‑41 (MQ=255)
                                                                ggTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGc      <  1:172290/71‑1 (MQ=255)
                                                                ggTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGc      <  1:412767/71‑1 (MQ=255)
                                                                ggTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGc      <  1:36380/71‑1 (MQ=255)
                                                                  tACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTg    <  1:632821/71‑1 (MQ=255)
                                                                    cAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGc  <  1:158314/71‑1 (MQ=255)
                                                                    cAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGc  <  1:153496/71‑1 (MQ=255)
                                                                    cAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGc  <  1:2347025/71‑1 (MQ=255)
                                                                     aGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGc  <  1:3795209/70‑1 (MQ=255)
                                                                      gTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTAt        <  1:3923861/63‑1 (MQ=255)
                                                                      |                                                                    
CAGCAGACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGC  >  minE/836376‑836514

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: