Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,847 G→C 8.9% V769L (GTG→CTG)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,8470GC 8.9% 156.6 / 5.7 56V769L (GTG→CTG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (30/21);  new base C (2/3);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTG  >  minE/833777‑833917
                                                                      |                                                                      
gCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtgg                                                                        <  1:1579500/71‑1 (MQ=255)
   aCAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtggtg                                                                      <  1:3848923/70‑1 (MQ=255)
       ggcggcGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTAc                                                                 >  1:1042999/1‑71 (MQ=255)
       ggcggcGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTAc                                                                 >  1:353916/1‑71 (MQ=255)
       ggcggcGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTAc                                                                 >  1:93836/1‑71 (MQ=255)
             gTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGa                                                           >  1:3922503/1‑71 (MQ=255)
             gTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGa                                                           >  1:2099251/1‑71 (MQ=255)
              tGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtgg                                                                        <  1:218400/57‑1 (MQ=255)
              tGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGc                                                               >  1:1420810/1‑66 (MQ=255)
                aaGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACtt                                                        >  1:2965366/1‑71 (MQ=255)
                 aGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtggtg                                                                      >  1:1121658/1‑56 (MQ=255)
                 aGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtggtg                                                                      >  1:1397650/1‑56 (MQ=255)
                 aGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACtt                                                        <  1:662166/70‑1 (MQ=255)
                     ggAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTc                                                   <  1:1820759/71‑1 (MQ=255)
                          aagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTc                                               <  1:3367957/70‑1 (MQ=255)
                            gTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAg                                            >  1:170765/1‑71 (MQ=255)
                            gTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAg                                            <  1:1867046/71‑1 (MQ=255)
                               gACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCAc                                         <  1:353209/71‑1 (MQ=255)
                                   ggCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGt                                     >  1:2150872/1‑71 (MQ=255)
                                    gCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTg                                      <  1:3879027/69‑1 (MQ=255)
                                    gCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGtc                                    <  1:2423677/71‑1 (MQ=255)
                                     cAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGtct                                   <  1:841087/71‑1 (MQ=255)
                                       ggACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTg                                                             <  1:3636206/43‑1 (MQ=255)
                                       ggACAATGGTCTGGAAGGCAAGCTGGATCTGCTGGTTACGCTg                                                             <  1:1559453/43‑1 (MQ=255)
                                        gACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTg                                                 <  1:2608539/54‑1 (MQ=255)
                                        gACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTa                                >  1:1941801/1‑71 (MQ=255)
                                         aCAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACtt                                                        <  1:920669/46‑1 (MQ=255)
                                          cAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTa                                                                  >  1:230671/1‑35 (MQ=255)
                                          cAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTa                                                                  >  1:1091855/1‑35 (MQ=255)
                                          cAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGtctc                                  >  1:3766115/1‑67 (MQ=255)
                                           aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTcc                                                      >  1:2844672/1‑46 (MQ=255)
                                           aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTcc                                                      >  1:3195794/1‑46 (MQ=255)
                                           aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTcc                                                      >  1:2038522/1‑46 (MQ=255)
                                               gTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa                          >  1:517880/1‑70 (MQ=255)
                                               gTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa                          >  1:628468/1‑70 (MQ=255)
                                                tCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCg                                              <  1:1354309/49‑1 (MQ=255)
                                                 tttgAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGa                          <  1:1693588/65‑1 (MQ=255)
                                                 attgAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAGTACATGCCTGTTCTCCGatat                       <  1:78509/68‑1 (MQ=255)
                                                  tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcat                       >  1:3906359/1‑70 (MQ=255)
                                                  tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcat                       >  1:1794463/1‑70 (MQ=255)
                                                    gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCAc                                         >  1:562351/1‑50 (MQ=255)
                                                    gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca                     <  1:2543467/70‑1 (MQ=255)
                                                    gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcataat                    >  1:1805447/1‑71 (MQ=255)
                                                    gAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAg                                            >  1:1254039/1‑47 (MQ=255)
                                                    gAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAg                                            >  1:2319323/1‑47 (MQ=255)
                                                     aaGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCa                                          <  1:843662/48‑1 (MQ=255)
                                                     aaGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTc                             >  1:3540995/1‑61 (MQ=255)
                                                      aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAttt                  <  1:2103152/71‑1 (MQ=255)
                                                      aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAttt                  <  1:2738845/71‑1 (MQ=255)
                                                      aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAtt                   >  1:2705971/1‑70 (MQ=255)
                                                       ggCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAtttt                 <  1:3019309/71‑1 (MQ=255)
                                                        gCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGtctc                                  <  1:142174/53‑1 (MQ=255)
                                                        gCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg                           <  1:2905222/60‑1 (MQ=255)
                                                        gCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTg                >  1:465022/1‑71 (MQ=255)
                                                         cAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTa                                >  1:337403/1‑54 (MQ=255)
                                                           aGCTGGATCTGGTGGTTACCCTGGACTTCCGTCTGTCGAGCACCTGTCTCTa                                >  1:2885293/1‑52 (MQ=255)
                                                            gCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGa            >  1:996399/1‑71 (MQ=255)
                                                             cTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGAc           <  1:2875369/71‑1 (MQ=255)
                                                                  tCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTAt                               <  1:476899/46‑1 (MQ=255)
                                                                  tCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGa      >  1:1051907/1‑71 (MQ=255)
                                                                   cTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGAc     >  1:1257131/1‑71 (MQ=255)
                                                                    tggtggTTACGCTGGACTTCCGTCTGTCGAGCACCTGtctc                                  <  1:2172931/41‑1 (MQ=255)
                                                                    tggtggTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGAc           >  1:2466527/1‑64 (MQ=255)
                                                                      gtggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt   <  1:3694872/70‑1 (MQ=255)
                                                                      gtggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTg  >  1:1699652/1‑71 (MQ=255)
                                                                      |                                                                      
GCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTG  >  minE/833777‑833917

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: