Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,868 C→A 13.4% L776M (CTG→ATG)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,8680CA13.4% 197.1 / 10.1 66L776M (CTG→ATG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base A (6/3);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGA  >  minE/833798‑833929
                                                                      |                                                             
ggAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTc                                                               <  1:1820759/71‑1 (MQ=255)
     aagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTc                                                           <  1:3367957/70‑1 (MQ=255)
       gTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAg                                                        <  1:1867046/71‑1 (MQ=255)
       gTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAg                                                        >  1:170765/1‑71 (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)
                   gACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTg                                                             <  1:2608539/54‑1 (MQ=255)
                   gACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTa                                            >  1:1941801/1‑71 (MQ=255)
                     cAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGtctc                                              >  1:3766115/1‑67 (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:1794463/1‑70 (MQ=255)
                             tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcat                                   >  1:3906359/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:2319323/1‑47 (MQ=255)
                               gAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAg                                                        >  1:1254039/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)
                                                  tggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg             <  1:2478699/71‑1 (MQ=255)
                                                   ggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg             >  1:2571116/1‑70 (MQ=255)
                                                   ggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg             >  1:2964215/1‑70 (MQ=255)
                                                    gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg                                       >  1:871411/1‑43 (MQ=255)
                                                    gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg                                       >  1:2511607/1‑43 (MQ=255)
                                                    gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGt            <  1:1800876/70‑1 (MQ=255)
                                                       aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTAt                                           >  1:260950/1‑36 (MQ=255)
                                                       aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTg              >  1:2002814/1‑65 (MQ=255)
                                                       aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTg              >  1:1969392/1‑65 (MQ=255)
                                                        cGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaa       <  1:1824569/71‑1 (MQ=255)
                                                         gCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAc                                     <  1:704334/40‑1 (MQ=255)
                                                         gCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaa      <  1:1861217/71‑1 (MQ=255)
                                                         gCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa      <  1:2948002/71‑1 (MQ=255)
                                                           tcgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa      >  1:3614129/3‑69 (MQ=255)
                                                           tcgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaa       >  1:1667589/3‑68 (MQ=255)
                                                           tGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGAc                 >  1:1817918/1‑58 (MQ=255)
                                                            cgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGACACCTGGTACGaaa      >  1:2271443/2‑68 (MQ=255)
                                                             gACTTCCGTCTGTCGAGCACGTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaaa  <  1:1768595/70‑2 (MQ=255)
                                                             gACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAga  >  1:1895780/1‑71 (MQ=255)
                                                                ttCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca                                 <  1:3573381/37‑1 (MQ=255)
                                                                ttCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca                                 <  1:3568828/37‑1 (MQ=255)
                                                                ttCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt               <  1:2947024/55‑1 (MQ=255)
                                                                ttCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt               <  1:1625578/55‑1 (MQ=255)
                                                                      |                                                             
GGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGA  >  minE/833798‑833929

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: