Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,226 T→G 24.6% S228S (TCT→TCG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,2260TG24.6% 87.0 / 47.4 69S228S (TCT→TCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (22/30);  new base G (8/9);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCA  >  minE/832156‑832293
                                                                      |                                                                   
gCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                                                     >  1:1003353/1‑71 (MQ=255)
gCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                                                     >  1:3162483/1‑71 (MQ=255)
gCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                                                     >  1:3778550/1‑71 (MQ=255)
gCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                                                     >  1:510064/1‑71 (MQ=255)
 cTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtctc                                                                    <  1:1105305/71‑1 (MQ=255)
   aTCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGcc                                                                   <  1:2645165/70‑1 (MQ=255)
   aTCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGcc                                                                   >  1:3292544/1‑70 (MQ=255)
    tctgtcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGcc                                                                   <  1:3875853/69‑1 (MQ=255)
     ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                                                 >  1:1034909/1‑70 (MQ=255)
     ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                                                 >  1:2059903/1‑70 (MQ=255)
     ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                                                 >  1:1599583/1‑70 (MQ=255)
       gtcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCgat                                                                >  1:214126/2‑67 (MQ=255)
        tcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                                                               <  1:2128473/69‑1 (MQ=255)
        tcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAc                                                             >  1:2177040/1‑71 (MQ=255)
        tcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAc                                                             >  1:3680238/1‑71 (MQ=255)
          gCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtctc                                                                    <  1:3811769/62‑1 (MQ=255)
           cTGATTGGCTGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTg                                                          <  1:3590913/71‑1 (MQ=255)
           cTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTg                                                          <  1:1659438/71‑1 (MQ=255)
             gATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTggg                                                        >  1:1951126/1‑71 (MQ=255)
              aTTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTgggg                                                       >  1:3435915/1‑71 (MQ=255)
               ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGc                                                      >  1:2810613/1‑71 (MQ=255)
               ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGc                                                      <  1:3029009/71‑1 (MQ=255)
                ttgCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCg                                                     <  1:3865011/69‑1 (MQ=255)
                  gCGGTACTTGCTTAAGCTTCTACGACTGGTACAGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAg                                                   <  1:2912485/71‑1 (MQ=255)
                  gCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAg                                                   >  1:1652358/1‑71 (MQ=255)
                    ggTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCa                                                                >  1:2643338/1‑56 (MQ=255)
                     gTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                                                 >  1:3821945/1‑70 (MQ=255)
                      tACTTGCTTAAGCTTCTACGACTCGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTggg                                                        <  1:869607/62‑1 (MQ=255)
                        cTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACt                                             <  1:3917703/71‑1 (MQ=255)
                         ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaaa                                                              >  1:2709945/1‑53 (MQ=255)
                          tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                                                     >  1:247265/1‑45 (MQ=255)
                           gCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTg                                            <  1:700467/69‑1 (MQ=255)
                            cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTg                                            <  1:1064268/68‑1 (MQ=255)
                            cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACg                                         <  1:1482933/71‑1 (MQ=255)
                              tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTa                                       >  1:1842024/1‑71 (MQ=255)
                               aaGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                                                     <  1:824904/40‑1 (MQ=255)
                                 gCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtctc                                                                    >  1:824643/1‑39 (MQ=255)
                                 gCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                                                     <  1:2590855/38‑1 (MQ=255)
                                 gCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                   >  1:81729/1‑40 (MQ=255)
                                 gCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                   >  1:2858709/1‑40 (MQ=255)
                                  cTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCTCAAACCTGGGGCGAGCAAACTGACGTACCgg                                   <  1:375734/71‑1 (MQ=255)
                                  cTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTg                                                          <  1:1304482/48‑1 (MQ=255)
                                    tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGc                                                      <  1:1338106/50‑1 (MQ=255)
                                    tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCa                                                 <  1:3075661/55‑1 (MQ=255)
                                    tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGa                                  >  1:2231834/1‑70 (MQ=255)
                                    tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGa                                  <  1:3881838/70‑1 (MQ=255)
                                    tCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGaa                                 >  1:902386/1‑71 (MQ=255)
                                     cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                                                               <  1:3819244/40‑1 (MQ=255)
                                     cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                                                               <  1:2584082/40‑1 (MQ=255)
                                     cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                  <  1:877986/69‑1 (MQ=255)
                                     cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                >  1:1349517/1‑71 (MQ=255)
                                      tACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGa                                                    >  1:2636449/1‑50 (MQ=255)
                                      tACGACTGGTATTGCGACCTGCCGCCCGCATCGCCGATGACCTGGGGCg                                                     <  1:2179961/49‑1 (MQ=255)
                                       aCGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                >  1:3906156/1‑69 (MQ=255)
                                        cGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATctg                             >  1:2893140/1‑71 (MQ=255)
                                        cGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATct                              <  1:3163374/70‑1 (MQ=255)
                                        cGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                <  1:3273029/68‑1 (MQ=255)
                                           cTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCg                                                     <  1:2402401/44‑1 (MQ=255)
                                           cTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATctgctg                          <  1:2482226/71‑1 (MQ=255)
                                                attGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATctgc                            <  1:3019215/62‑1 (MQ=255)
                                                attGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGACTgg                     >  1:1164402/3‑71 (MQ=255)
                                                  tGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                  <  1:3892563/56‑1 (MQ=255)
                                                   gCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                  <  1:3828782/55‑1 (MQ=255)
                                                   gCGACTTGCCTCATGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTa                   >  1:1000755/1‑70 (MQ=255)
                                                    cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc                  <  1:1217644/70‑1 (MQ=255)
                                                    cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc                  <  1:3645032/70‑1 (MQ=255)
                                                    cGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGACTGGTat                  <  1:3768923/70‑2 (MQ=255)
                                                        ttGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTc             >  1:2335376/1‑71 (MQ=255)
                                                          gCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTgg                     <  1:2292016/61‑1 (MQ=255)
                                                          gCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTgg                     <  1:313955/61‑1 (MQ=255)
                                                           cctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCa           >  1:1639362/1‑70 (MQ=255)
                                                             tcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCa           <  1:2042996/68‑1 (MQ=255)
                                                             tcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCa           <  1:2338139/68‑1 (MQ=255)
                                                               ctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAc      <  1:2468480/71‑1 (MQ=255)
                                                                 gCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcat    >  1:2493073/1‑71 (MQ=255)
                                                                  cGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcat    >  1:2353916/1‑70 (MQ=255)
                                                                  cGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAca     <  1:3653117/69‑1 (MQ=255)
                                                                   gTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcatca  <  1:784003/71‑1 (MQ=255)
                                                                   gTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcatc   >  1:3764119/1‑70 (MQ=255)
                                                                    tctcCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTgg                     <  1:971659/51‑1 (MQ=255)
                                                                      |                                                                   
GCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCA  >  minE/832156‑832293

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: