Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,858 G→C 7.6% L772L (CTG→CTC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,8580GC 7.6% 194.5 / 5.5 66L772L (CTG→CTCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (34/27);  new base C (2/3);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAA  >  minE/833790‑833925
                                                                    |                                                                   
gTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGa                                                                   >  1:2099251/1‑71 (MQ=255)
gTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGa                                                                   >  1:3922503/1‑71 (MQ=255)
   aaGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACtt                                                                >  1:2965366/1‑71 (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)
                             cAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGtctc                                          >  1:3766115/1‑67 (MQ=255)
                              aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTcc                                                              >  1:2038522/1‑46 (MQ=255)
                              aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTcc                                                              >  1:2844672/1‑46 (MQ=255)
                              aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTcc                                                              >  1:3195794/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: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:2738845/71‑1 (MQ=255)
                                         aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAttt                          <  1:2103152/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:2511607/1‑43 (MQ=255)
                                                            gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg                                   >  1:871411/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)
                                                                   tGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGAc             >  1:1817918/1‑58 (MQ=255)
                                                                    |                                                                   
GTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAA  >  minE/833790‑833925

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: