Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,330 C→T 41.5% C16C (TGC→TGT 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,3300CT41.5% 70.5 / 133.2 123C16C (TGC→TGTnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (30/42);  new base T (23/28);  total (53/70)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 123 total aligned reads displayed.

ACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGC  >  minE/835261‑835398
                                                                      |                                                                    
aCCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                                                                     >  1:233830/1‑71 (MQ=255)
 tcAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTAt                                                                    >  1:2954352/2‑71 (MQ=255)
 tcAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTAt                                                                    >  1:3005167/2‑71 (MQ=255)
 ccAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa                                                                     <  1:3122629/70‑1 (MQ=255)
  cAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATc                                                                   >  1:2118422/1‑71 (MQ=255)
  cAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATc                                                                   >  1:556670/1‑71 (MQ=255)
  cAGGTACAGGAGGCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATc                                                                   >  1:4734300/1‑71 (MQ=255)
   aGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCg                                                                  >  1:1446184/1‑71 (MQ=255)
       aCAGGAGAGCGTAAAATGAGAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg                                                              >  1:2227744/1‑71 (MQ=255)
       aCAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg                                                              >  1:1495309/1‑71 (MQ=255)
           gagaGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCcac                                                          <  1:3211809/71‑1 (MQ=255)
              agCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACAcc                                                       <  1:104920/71‑1 (MQ=255)
                  aaaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGt                                                    <  1:892720/69‑1 (MQ=255)
                  aaaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGt                                                    <  1:2601292/69‑1 (MQ=255)
                   aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGc                                                             >  1:458403/2‑61 (MQ=255)
                   aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGc                                                             >  1:124944/2‑61 (MQ=255)
                   aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGtt                                                   <  1:3431545/70‑1 (MQ=255)
                   aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGATTAACCTCGAGAAATGTATCGGCTGCCATACCTGtt                                                   <  1:515880/70‑1 (MQ=255)
                    aaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTAt                                                                    >  1:3815755/1‑53 (MQ=255)
                     aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                                  <  1:4398764/70‑1 (MQ=255)
                     aaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTcg                                                 >  1:2178191/1‑70 (MQ=255)
                        tGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                                  <  1:2007373/67‑1 (MQ=255)
                          aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTTCCACACCTGTTCAGTTAc                                            >  1:1753273/1‑71 (MQ=255)
                          aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTa                                             <  1:2111745/70‑1 (MQ=255)
                          aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAc                                            >  1:2988737/1‑71 (MQ=255)
                           aaaTCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGAcc                                           >  1:3033300/1‑71 (MQ=255)
                            aaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                                  <  1:643041/63‑1 (MQ=255)
                            aaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAcc                                           <  1:5463800/70‑1 (MQ=255)
                            aaTCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCt                                          <  1:3422085/71‑1 (MQ=255)
                             aTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt                                                                    >  1:2679041/1‑44 (MQ=255)
                             aTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt                                                                    >  1:4538003/1‑44 (MQ=255)
                             aTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt                                                                    >  1:1975393/1‑44 (MQ=255)
                             aTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAg                                                >  1:5531578/1‑64 (MQ=255)
                              ttCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGt                                        <  1:1202460/71‑1 (MQ=255)
                              tccGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCt                                          >  1:2955088/3‑69 (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)
                                  ttCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTTTaaaa                                    <  1:1391440/71‑1 (MQ=255)
                                  ttCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaaa                                    <  1:823694/71‑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)
                                          tCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTgg                            <  1:3266844/71‑1 (MQ=255)
                                          tCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTg                             <  1:4548543/70‑1 (MQ=255)
                                            ggCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGAc                          <  1:4918992/71‑1 (MQ=255)
                                               aTGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGAccg                        <  1:401755/70‑2 (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)
                                                    gCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTa                                       >  1:2454457/1‑50 (MQ=255)
                                                    gCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTa                                       >  1:921225/1‑50 (MQ=255)
                                                    gCTGAATCTCGATAAGTGCATCGGCTTCCACACCTGTTCAGTTACCTGTaa                                      <  1:757831/51‑1 (MQ=255)
                                                    gCTGAATCTCGATAAGTGCATCGGCTGCCACAGCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGt                   >  1:875915/1‑70 (MQ=255)
                                                    gCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGt                   >  1:4406792/1‑70 (MQ=255)
                                                    gCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGt                   >  1:4773351/1‑70 (MQ=255)
                                                     cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGt                   >  1:4223772/1‑69 (MQ=255)
                                                     cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                  <  1:4804319/70‑1 (MQ=255)
                                                     cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGa                 <  1:2550486/71‑1 (MQ=255)
                                                        aaTCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAAc                                  >  1:2485647/1‑51 (MQ=255)
                                                        aaTCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAAc                                  >  1:557873/1‑51 (MQ=255)
                                                        aaTCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAAc                                  >  1:290336/1‑51 (MQ=255)
                                                        aaTCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGaa                <  1:1816674/69‑1 (MQ=255)
                                                         acctcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAgg              <  1:1544852/68‑1 (MQ=255)
                                                         aTCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTTAAGGc             <  1:2791370/71‑1 (MQ=255)
                                                          tctcGATAAGTGCATCGGCTGCCACACCTGTTCAGt                                               <  1:2025118/36‑1 (MQ=255)
                                                          tctcGATAAGTGCATCGGCTGCCACACCTGTTCAGt                                               <  1:2429816/36‑1 (MQ=255)
                                                          tctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGt                                <  1:4109947/51‑1 (MQ=255)
                                                          tctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGt                                <  1:3492849/51‑1 (MQ=255)
                                                          cctcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGca            <  1:1450542/70‑2 (MQ=255)
                                                            tcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTa                                             >  1:690910/1‑36 (MQ=255)
                                                            tcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTgg                            >  1:511548/1‑53 (MQ=255)
                                                            tcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTaaa                                     <  1:1253319/44‑1 (MQ=16)
                                                            tcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGc             <  1:1833308/68‑1 (MQ=255)
                                                            tcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGc             <  1:1188328/68‑1 (MQ=255)
                                                            tcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGc             <  1:4446910/68‑1 (MQ=255)
                                                             cGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg                             >  1:4888898/1‑51 (MQ=255)
                                                              gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGaa                >  1:225986/1‑63 (MQ=255)
                                                              gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg         >  1:2555624/1‑70 (MQ=255)
                                                              gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg         <  1:757757/70‑1 (MQ=255)
                                                              gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg         <  1:2055536/70‑1 (MQ=255)
                                                              gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg         >  1:5387749/1‑70 (MQ=255)
                                                              gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGCACCAGCCGTGAAGGCGTgg         <  1:4789397/70‑1 (MQ=255)
                                                               aTAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAAc                                  <  1:3557967/44‑1 (MQ=255)
                                                               aTAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAAc                                  <  1:1283786/44‑1 (MQ=255)
                                                               aTAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAAc                                  <  1:2326283/44‑1 (MQ=255)
                                                               aTAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAAc                                  <  1:2184372/44‑1 (MQ=255)
                                                                  aGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata     <  1:5479613/70‑1 (MQ=255)
                                                                   gTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATg   <  1:1222241/71‑1 (MQ=255)
                                                                    tGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGc  <  1:2175775/71‑1 (MQ=255)
                                                                      |                                                                    
ACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGC  >  minE/835261‑835398

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: