Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,217 G→C 5.5% V1225V (GTG→GTC) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,2170GC 5.5% 286.1 / 4.9 90V1225V (GTG→GTC) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (39/46);  new base C (3/2);  total (42/48)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

ATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATG  >  minE/835147‑835285
                                                                      |                                                                    
aTGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAAGCGCGATGAGTTTGTTGTAGTg                                                                      >  1:71732/1‑71 (MQ=255)
aTGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTg                                                                      >  1:3152808/1‑71 (MQ=255)
    ccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTa                                                                  <  1:230824/71‑1 (MQ=255)
      aTCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGt                                                                   >  1:1865410/1‑68 (MQ=255)
       tCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCg                                                                    >  1:3939148/1‑66 (MQ=255)
        cTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGa                                                               >  1:1912363/1‑70 (MQ=255)
         tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTaa                                                                 <  1:1844840/67‑1 (MQ=255)
         tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGAt                                                              >  1:192605/1‑70 (MQ=255)
           gCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGAt                                                              <  1:2851865/68‑1 (MQ=255)
           gCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaa                                                           <  1:604381/71‑1 (MQ=255)
            cATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaag                                                          <  1:4072356/71‑1 (MQ=255)
             aTACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaaga                                                         <  1:182280/71‑1 (MQ=255)
                gggCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAc                                                       <  1:2679682/69‑1 (MQ=255)
                cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                                        >  1:601140/1‑69 (MQ=255)
                cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                                        >  1:3651660/1‑69 (MQ=255)
                cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                                                       <  1:3309420/70‑1 (MQ=255)
                cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                                                       <  1:2453992/70‑1 (MQ=255)
                cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                      <  1:1816317/71‑1 (MQ=255)
                cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                      >  1:1743056/1‑71 (MQ=255)
                cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                      >  1:2788969/1‑71 (MQ=255)
                 ggCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                                        <  1:232103/68‑1 (MQ=255)
                  gCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                                                    <  1:674326/71‑1 (MQ=255)
                  gCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGtt                                                    <  1:1072415/71‑1 (MQ=255)
                   cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                      >  1:2875070/1‑68 (MQ=255)
                   cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTg                                                   <  1:1492233/71‑1 (MQ=255)
                       aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                               >  1:1524225/1‑71 (MQ=255)
                       aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                               >  1:3964302/1‑71 (MQ=255)
                       aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                               >  1:1446072/1‑71 (MQ=255)
                        actactATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTgg                                              >  1:3977686/1‑71 (MQ=255)
                           actaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                             >  1:2601759/1‑69 (MQ=255)
                           actaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcg                                                      >  1:2101860/1‑59 (MQ=255)
                           actaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcg                                                      >  1:761130/1‑59 (MQ=255)
                            ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                             >  1:3308921/1‑68 (MQ=255)
                            ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                             >  1:436350/1‑68 (MQ=255)
                            ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                           <  1:3464078/70‑1 (MQ=255)
                            ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                           <  1:1400548/70‑1 (MQ=255)
                               tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATga                                                            <  1:3062234/50‑1 (MQ=255)
                               tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                                                       >  1:3293592/1‑55 (MQ=255)
                               cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATgaagaa                                                        >  1:1825414/2‑54 (MQ=255)
                                 gCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                                        >  1:3418545/1‑52 (MQ=255)
                                  cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg                                    <  1:3299478/71‑1 (MQ=255)
                                  cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg                                    <  1:834316/71‑1 (MQ=255)
                                    ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                      <  1:70176/51‑1 (MQ=255)
                                    ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg                                       <  1:2883082/66‑1 (MQ=255)
                                    ccGTAGGTTCTAACCGCGATGAGTTGTTGTAGTGCGTAAGATGAAGAACa                                                      <  1:3281758/50‑1 (MQ=255)
                                       tAGGTTCTAACCGCGATGAGTTTGTTGTATTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAgg                                >  1:4022729/1‑70 (MQ=255)
                                       tAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACt                                                >  1:2933549/1‑54 (MQ=255)
                                       tAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAgg                                >  1:233002/1‑70 (MQ=255)
                                          gTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGa                                   <  1:2801162/64‑1 (MQ=255)
                                            tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                          >  1:2495360/1‑71 (MQ=255)
                                            tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                          >  1:3131654/1‑71 (MQ=255)
                                            tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                          >  1:2742322/1‑71 (MQ=255)
                                             cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                         >  1:3203561/1‑55 (MQ=255)
                                             cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                          <  1:4099664/70‑1 (MQ=255)
                                             cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                         >  1:1286771/1‑71 (MQ=255)
                                             cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                         >  1:3929115/1‑71 (MQ=255)
                                              tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                                                    <  1:2065101/43‑1 (MQ=255)
                                              tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACt                                                <  1:2998402/47‑1 (MQ=255)
                                              tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGc                               <  1:1360076/64‑1 (MQ=255)
                                              tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGc                               <  1:1604771/64‑1 (MQ=255)
                                              tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATg                           <  1:3862420/68‑1 (MQ=255)
                                               aaCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                         >  1:1319084/1‑69 (MQ=255)
                                               aaCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                         >  1:1926929/1‑69 (MQ=255)
                                               aaCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                         >  1:3994764/1‑69 (MQ=255)
                                                aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                       >  1:716929/1‑70 (MQ=255)
                                                aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                       <  1:2917359/70‑1 (MQ=255)
                                                aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                       <  1:2785046/70‑1 (MQ=255)
                                                aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAg                      >  1:992980/1‑71 (MQ=255)
                                                 ccGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAgg                     <  1:756013/71‑1 (MQ=255)
                                                 ccGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAgg                     <  1:1514824/71‑1 (MQ=255)
                                                 ccGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAgg                     <  1:2236259/71‑1 (MQ=255)
                                                  cgcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg                                    <  1:3237832/55‑1 (MQ=255)
                                                   gcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGAt                                        <  1:1387209/50‑1 (MQ=255)
                                                   gcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGaa                                  >  1:2899431/1‑56 (MQ=255)
                                                   gcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                       >  1:3286444/1‑67 (MQ=255)
                                                   gcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                       >  1:3551702/1‑67 (MQ=255)
                                                    cgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTa                   >  1:3874902/1‑70 (MQ=255)
                                                      aTGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAgg                     <  1:792997/66‑1 (MQ=255)
                                                        gAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAAt                            >  1:4036455/1‑57 (MQ=255)
                                                         aGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGga              <  1:265130/70‑1 (MQ=255)
                                                           tttGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                           <  1:3866701/39‑1 (MQ=255)
                                                           tttGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAg                                          <  1:529778/40‑1 (MQ=255)
                                                           tttGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGagag           <  1:1979568/71‑1 (MQ=255)
                                                             tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCg         >  1:888426/1‑71 (MQ=255)
                                                                 gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa      <  1:1597182/70‑1 (MQ=255)
                                                                 gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa      <  1:4134512/70‑1 (MQ=255)
                                                                 gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa      <  1:4080190/70‑1 (MQ=255)
                                                                  tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaaaa    <  1:804020/71‑1 (MQ=255)
                                                                  tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa      <  1:3671209/69‑1 (MQ=255)
                                                                    gTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCa                              <  1:371431/43‑1 (MQ=255)
                                                                    gTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAAt   >  1:981673/1‑70 (MQ=255)
                                                                    gTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATg  <  1:323748/71‑1 (MQ=255)
                                                                      gCGTAAGATGAAGAACATTGACTGGTTAGATGGCg                                    <  1:3822226/35‑1 (MQ=255)
                                                                      |                                                                    
ATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATG  >  minE/835147‑835285

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: