Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,315 G→T 39.9% L11L (CTG→CTT narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,3150GT39.9% 146.5 / 125.6 135L11L (CTG→CTTnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (41/40);  new base T (29/25);  total (70/65)
Fisher's exact test for biased strand distribution p-value = 8.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

GATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGA  >  minE/835245‑835383
                                                                       |                                                                    
gATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTg                                                                      <  1:236991/71‑1 (MQ=255)
 aTGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGa                                                                     <  1:2242913/71‑1 (MQ=255)
    gCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtc                                                                  <  1:1585040/71‑1 (MQ=255)
      gAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATTGTGCTGAAtctc                                                                <  1:5349937/71‑1 (MQ=255)
       aaGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCg                                                               <  1:3902747/71‑1 (MQ=255)
         ggCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                                                             <  1:4381332/71‑1 (MQ=255)
            aaTGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAg                                                          >  1:4613436/1‑71 (MQ=255)
            aaTGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAg                                                          >  1:4002944/1‑71 (MQ=255)
            aaTGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGGGCTGAATCTCGATaa                                                           <  1:5613257/70‑1 (MQ=255)
             aTGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAg                                                          >  1:1670486/1‑70 (MQ=255)
                 tcAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTAt                                                     >  1:3005167/2‑71 (MQ=255)
                 tcAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTAt                                                     >  1:2954352/2‑71 (MQ=255)
                 ccAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                                                      <  1:3122629/70‑1 (MQ=255)
                  cAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATc                                                    >  1:556670/1‑71 (MQ=255)
                  cAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATc                                                    >  1:2118422/1‑71 (MQ=255)
                   aGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCg                                                   >  1:1446184/1‑71 (MQ=255)
                   aGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACaaat                                                         >  1:2580379/1‑63 (MQ=255)
                   aGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACaaat                                                         >  1:2366678/1‑63 (MQ=255)
                   aGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACaaat                                                         >  1:2155698/1‑63 (MQ=255)
                     gTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt                                                             <  1:4453539/59‑1 (MQ=255)
                       aCAGGAGAGCGTAAAATGAGAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg                                               >  1:2227744/1‑71 (MQ=255)
                       aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg                                               >  1:1495309/1‑71 (MQ=255)
                           gagaGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCcac                                           <  1:3211809/71‑1 (MQ=255)
                            agagCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct                                                                 >  1:3287573/1‑48 (MQ=255)
                              agCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACAcc                                        <  1:104920/71‑1 (MQ=255)
                                cGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGa                                                                     >  1:3072600/1‑40 (MQ=255)
                                cGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGa                                                                     >  1:2775739/1‑40 (MQ=255)
                                  aaaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGt                                     <  1:892720/69‑1 (MQ=255)
                                  aaaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGt                                     <  1:2601292/69‑1 (MQ=255)
                                   tAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGtt                                    <  1:829840/71‑1 (MQ=255)
                                   aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGc                                              >  1:458403/2‑61 (MQ=255)
                                   aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGc                                              >  1:1924013/2‑61 (MQ=255)
                                   aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGc                                              >  1:124944/2‑61 (MQ=255)
                                   aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGtt                                    <  1:4462691/70‑1 (MQ=255)
                                   aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGtt                                    <  1:3431545/70‑1 (MQ=255)
                                    aaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTAt                                                     >  1:3815755/1‑53 (MQ=255)
                                    aaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTc                                   >  1:2363345/1‑71 (MQ=255)
                                     aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGaa                                                                    >  1:3377056/1‑37 (MQ=255)
                                     aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGaa                                                                    >  1:3621258/1‑37 (MQ=255)
                                     aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGa                                                                     <  1:5326309/36‑1 (MQ=255)
                                     aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa                                                              >  1:4103550/1‑43 (MQ=255)
                                     aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa                                                              >  1:485090/1‑43 (MQ=255)
                                     aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATa                                                            <  1:2027004/45‑1 (MQ=255)
                                     aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATa                                                            <  1:5139248/45‑1 (MQ=255)
                                     aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATa                                                            <  1:4720835/45‑1 (MQ=255)
                                     aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGt                                                         <  1:4634673/48‑1 (MQ=255)
                                     aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGt                                                         <  1:3425474/48‑1 (MQ=255)
                                     aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                   <  1:4398764/70‑1 (MQ=255)
                                     aaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTaac                                                                   >  1:4297083/1‑37 (MQ=38)
                                     aaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTaac                                                                   >  1:1661258/1‑37 (MQ=38)
                                     aaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTaac                                                                   >  1:1654626/1‑37 (MQ=38)
                                     aaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTcg                                  >  1:2178191/1‑70 (MQ=255)
                                      aaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATgt                                                       <  1:2977072/49‑2 (MQ=255)
                                      aaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATgt                                                       <  1:1707439/49‑2 (MQ=255)
                                       aTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGt                                >  1:2942982/1‑71 (MQ=255)
                                        tGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                   <  1:2007373/67‑1 (MQ=255)
                                        tGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATg                                                        <  1:3807655/46‑1 (MQ=255)
                                        tGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATg                                                        <  1:1345615/46‑1 (MQ=255)
                                          aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTTCCACACCTGTTCAGTTAc                             >  1:1753273/1‑71 (MQ=255)
                                          aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAc                             >  1:2988737/1‑71 (MQ=255)
                                          aaaaTCCGTTCACAAGTCGGCATGGTGCTTAACCTCGaca                                                            <  1:5576613/40‑3 (MQ=255)
                                          aaaaTCCGTTCACAAGTCGGCATGGTGCTTAACCTCGaca                                                            <  1:5583323/40‑3 (MQ=255)
                                           aaaTCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGAcc                            >  1:3033300/1‑71 (MQ=255)
                                            aaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAg                                                          >  1:844935/1‑40 (MQ=255)
                                            aaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAg                                                          >  1:1941724/1‑40 (MQ=255)
                                            aaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAg                                                          >  1:4641487/1‑40 (MQ=255)
                                            aaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAg                                                          >  1:31596/1‑40 (MQ=255)
                                            aaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAcc                            <  1:5463800/70‑1 (MQ=255)
                                             aTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt                                                     >  1:4538003/1‑44 (MQ=255)
                                             aTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt                                                     >  1:2679041/1‑44 (MQ=255)
                                             aTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAg                                 >  1:5531578/1‑64 (MQ=255)
                                              ttCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGt                         <  1:1202460/71‑1 (MQ=255)
                                              tccGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGtg                               <  1:2521365/63‑2 (MQ=255)
                                               ccGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTg                          >  1:4929890/2‑69 (MQ=255)
                                                 gTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTaaa                      <  1:39445/71‑1 (MQ=255)
                                                 gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaa                       <  1:3824596/70‑1 (MQ=255)
                                                 gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaa                       <  1:2735648/70‑1 (MQ=255)
                                                  ttCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGtt                                    <  1:1244187/56‑1 (MQ=255)
                                                  ttCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                      <  1:2080557/70‑1 (MQ=255)
                                                   tCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGc                                              <  1:2813408/45‑1 (MQ=255)
                                                    cacaAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCca                                            >  1:5225096/1‑46 (MQ=255)
                                                    cacaAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAAc                   >  1:2548851/1‑71 (MQ=255)
                                                     acaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACg                  >  1:3183088/1‑71 (MQ=255)
                                                      caAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGtt                                    <  1:807893/52‑1 (MQ=255)
                                                       aaGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTaaaa                     >  1:5685662/1‑66 (MQ=255)
                                                        aGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCt               >  1:503056/1‑71 (MQ=255)
                                                          tCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTgg             <  1:5621635/71‑1 (MQ=255)
                                                            ggCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGAc           <  1:4918992/71‑1 (MQ=255)
                                                               aTGGTGCTGAATCTCTATAAGTGCATCGGCTGCCACAcc                                        <  1:1000773/39‑1 (MQ=255)
                                                                tGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGtt                                    <  1:1798066/42‑1 (MQ=255)
                                                                tGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTc                                   <  1:2512733/43‑1 (MQ=255)
                                                                tGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTc                                   <  1:836334/43‑1 (MQ=255)
                                                                tGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                   <  1:4865397/43‑1 (MQ=255)
                                                                  gTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTg              >  1:4608998/1‑62 (MQ=255)
                                                                  gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCg     <  1:696788/71‑1 (MQ=255)
                                                                    gCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTa                        >  1:921225/1‑50 (MQ=255)
                                                                    gCTGAATCTCGATAAGTGCATCGGCTTCCACACCTGTTCAGTTACCTGTaa                       <  1:757831/51‑1 (MQ=255)
                                                                    gCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGt    >  1:4773351/1‑70 (MQ=255)
                                                                     cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGt    >  1:4223772/1‑69 (MQ=255)
                                                                     cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg   <  1:4804319/70‑1 (MQ=255)
                                                                       |                                                                    
GATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGA  >  minE/835245‑835383

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: