Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,372 T→G 10.7% S610S (TCT→TCG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,3720TG10.7% 123.7 / 13.4 56S610S (TCT→TCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (20/29);  new base G (4/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATG  >  minE/833302‑833441
                                                                      |                                                                     
aaaCCGGCTGTCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct                                                                       >  1:3302308/1‑71 (MQ=255)
aaaCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct                                                                       <  1:1320267/71‑1 (MQ=255)
    cGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTAtt                                                                   <  1:152361/71‑1 (MQ=255)
      gCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTAttta                                                                 >  1:2680760/1‑70 (MQ=255)
          gCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                                              <  1:1040045/70‑1 (MQ=255)
           ttgCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACaa                                                            <  1:2131717/69‑1 (MQ=255)
               cGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATaa                                                            >  1:2773946/1‑67 (MQ=255)
                gCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACt                                                       >  1:775826/1‑71 (MQ=255)
                  tGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTc                                                      <  1:3227087/70‑1 (MQ=255)
                    gcctttGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccag                                                   >  1:3170020/4‑71 (MQ=255)
                    gcctttGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccag                                                   >  1:229897/4‑71 (MQ=255)
                     cGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagc                                                  <  1:254711/71‑1 (MQ=255)
                       tttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                                              <  1:613701/57‑1 (MQ=255)
                       tttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                                              <  1:2309997/57‑1 (MQ=255)
                         tGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTAt                                                                    >  1:597055/1‑49 (MQ=255)
                         tGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTcca                                                    >  1:3372621/1‑65 (MQ=255)
                         tGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt                                              >  1:1590333/1‑71 (MQ=255)
                           cccTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCa                                                         <  1:184800/58‑1 (MQ=255)
                            ccTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCta                                                               >  1:188770/1‑51 (MQ=255)
                             attGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTa                                                                     >  1:999774/2‑44 (MQ=255)
                              ttGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                                              <  1:832556/50‑1 (MQ=255)
                               tGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct                                                                       >  1:2879497/1‑40 (MQ=255)
                               tGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct                                                                       >  1:2874469/1‑40 (MQ=255)
                                gACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                                              >  1:22360/1‑48 (MQ=255)
                                gACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGtt                                        <  1:1186222/70‑1 (MQ=255)
                                gACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTa                                       >  1:3239879/1‑71 (MQ=255)
                                     gCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaaa                                  >  1:617901/1‑71 (MQ=255)
                                     gCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaaa                                  >  1:2589941/1‑71 (MQ=255)
                                        gCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt                                              <  1:1024171/56‑1 (MQ=255)
                                         cGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGt                              <  1:1574199/71‑1 (MQ=255)
                                         cGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGt                              <  1:3222984/71‑1 (MQ=255)
                                          gTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTc                             <  1:2144655/71‑1 (MQ=255)
                                             cGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt                                              >  1:685868/1‑51 (MQ=255)
                                             cGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt                                              >  1:2523068/1‑51 (MQ=255)
                                               gcgcGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggc                        <  1:3085662/71‑1 (MQ=255)
                                                cgcgCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccagcc                                                 <  1:529854/45‑1 (MQ=255)
                                                cgcgCCAAATGAACAGCACCTCGTTTTTCTACAATCATTCCAGCCAATGGCGCTATGaaaaag                               >  1:1769060/1‑60 (MQ=255)
                                                cgcgCCAAATGAACAGCACCTCGTTTTTCTACAATCATTCCAGCCAATGGCGCTATGaaaaag                               >  1:2353032/1‑60 (MQ=255)
                                                 gcgTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggcgg                      <  1:1408708/71‑1 (MQ=255)
                                                    tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaa                                   <  1:2789547/55‑1 (MQ=255)
                                                    tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggcgg                      >  1:4022660/1‑68 (MQ=255)
                                                    tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcg                          <  1:593061/64‑1 (MQ=255)
                                                    tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAag                   <  1:3654739/71‑1 (MQ=255)
                                                     cacaTGAACAGCACTTCTTATTTCTATAACCACTccagcc                                                 <  1:1225315/40‑1 (MQ=255)
                                                      acaTGAACAGCACTTCTTATTTCTATAACCACTcc                                                     >  1:1666908/1‑35 (MQ=255)
                                                      acaTGAACAGCACTTCTTATTTCTATAACCACTcc                                                     >  1:4035666/1‑35 (MQ=255)
                                                      acaTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAagag                 <  1:3111514/71‑1 (MQ=255)
                                                        aTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAaga                  <  1:3941187/68‑1 (MQ=255)
                                                         tGAACAGCACTTCTTATTTCTATAACCACTccagc                                                  <  1:2052576/35‑1 (MQ=255)
                                                         tGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcg                          <  1:4062324/59‑1 (MQ=255)
                                                         tGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGaa                    <  1:251319/65‑1 (MQ=255)
                                                         tGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTg              >  1:2453286/1‑71 (MQ=255)
                                                           aaCAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTa                                       >  1:1258885/1‑44 (MQ=255)
                                                            aCAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaa                                   <  1:2891241/47‑1 (MQ=255)
                                                              aGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTc         <  1:314178/71‑1 (MQ=255)
                                                               gCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTc         <  1:481431/70‑1 (MQ=255)
                                                               gCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCa        <  1:2750492/71‑1 (MQ=255)
                                                                cACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTg              >  1:2171872/1‑64 (MQ=255)
                                                                cACTTCATATTTCTATAACCACTCCAGCCAGTGGCGTTa                                       >  1:536322/1‑39 (MQ=255)
                                                                 aCTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCAcc      <  1:3387854/71‑1 (MQ=255)
                                                                    tcttATTTCTATAACCACTCCAGCCAGTGGCGTTATGaaa                                  >  1:2905927/1‑40 (MQ=255)
                                                                    tcttATTTCTATAACCACTCCAGCCAGTGGCGTTATGaaa                                  >  1:2813084/1‑40 (MQ=255)
                                                                     cttATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATg  >  1:635327/1‑71 (MQ=255)
                                                                     cttATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATg  >  1:837128/1‑71 (MQ=255)
                                                                     cttATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATg  >  1:1635795/1‑71 (MQ=255)
                                                                      ttATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAagag                 >  1:391556/1‑55 (MQ=255)
                                                                      |                                                                     
AAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATG  >  minE/833302‑833441

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: