Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,852 T→G 6.1% V770V (GTT→GTG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,8520TG 6.1% 170.3 / 4.9 66V770V (GTT→GTGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (36/26);  new base G (2/2);  total (38/28)
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.83e-01

GGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGT  >  minE/833784‑833919
                                                                    |                                                                   
ggcggcGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTAc                                                                   >  1:1042999/1‑71 (MQ=255)
ggcggcGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTAc                                                                   >  1:93836/1‑71 (MQ=255)
ggcggcGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTAc                                                                   >  1:353916/1‑71 (MQ=255)
      gTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGa                                                             >  1:2099251/1‑71 (MQ=255)
      gTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGa                                                             >  1:3922503/1‑71 (MQ=255)
       tGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGc                                                                 >  1:1420810/1‑66 (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: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)
                                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:2038522/1‑46 (MQ=255)
                                    aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTcc                                                        >  1:3195794/1‑46 (MQ=255)
                                    aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTcc                                                        >  1:2844672/1‑46 (MQ=255)
                                        gTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa                            >  1:628468/1‑70 (MQ=255)
                                        gTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa                            >  1:517880/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:2964215/1‑70 (MQ=255)
                                                                 ggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg   >  1:2571116/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)
                                                                    |                                                                   
GGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGT  >  minE/833784‑833919

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: