Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,071 T→C 25.8% C263C (TGT→TGC narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE836,0710TC25.8% 119.8 / 43.3 76C263C (TGT→TGCnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base C (9/10);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

AGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACC  >  minE/836001‑836140
                                                                      |                                                                     
aGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtgt                                                                       >  1:1209816/1‑71 (MQ=255)
      agaAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcggt                                                                 <  1:4389993/71‑1 (MQ=255)
      agaAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcgg                                                                  <  1:508103/70‑1 (MQ=255)
        aaatGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCgtgggtc                                                                >  1:777092/4‑70 (MQ=255)
         aatGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCgtgggtcgc                                                              <  1:1656730/69‑2 (MQ=255)
         aGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGGCGGTCGt                                                              <  1:523969/71‑1 (MQ=255)
          atGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCgtgggtcgca                                                             >  1:3327909/2‑69 (MQ=255)
           tGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTAt                                                            <  1:401398/71‑1 (MQ=255)
              aTCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATcct                                                         <  1:5139350/71‑2 (MQ=255)
              aTCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCg                                                         <  1:4443711/71‑1 (MQ=255)
              aTCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCg                                                         <  1:1523005/71‑1 (MQ=255)
              aTCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAAACTGTGTCGGTCGTATCCg                                                         <  1:566046/71‑1 (MQ=255)
                 ttctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTa                                                      >  1:5082527/1‑71 (MQ=255)
                  tctGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTAt                                                     >  1:3191624/1‑71 (MQ=255)
                  tctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTAt                                                     >  1:58117/1‑71 (MQ=255)
                   ctGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCgtgggtcgca                                                             <  1:2475869/62‑3 (MQ=255)
                   ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGtt                                                       <  1:5263314/68‑1 (MQ=255)
                   ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGtt                                                       <  1:675981/68‑1 (MQ=255)
                   ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTa                                                      <  1:5161340/69‑1 (MQ=255)
                      cTATCCGCGTATTTAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCt                                                   <  1:4448824/69‑1 (MQ=255)
                         tCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCtt                                                  <  1:2319512/67‑1 (MQ=255)
                         tCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCtt                                                  <  1:2269130/67‑1 (MQ=255)
                           cgcgTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGt                                             <  1:1879847/70‑1 (MQ=255)
                            gcgTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGc                                           >  1:2062906/1‑71 (MQ=255)
                             cgTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCt                                          >  1:208666/1‑71 (MQ=255)
                                aTTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCt                                          >  1:2992793/1‑68 (MQ=255)
                                aTTGAAGCGGGTCAGACGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgtt                                       >  1:2292953/1‑71 (MQ=255)
                                 ttGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                      >  1:5677282/1‑71 (MQ=255)
                                  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                      >  1:3838078/1‑70 (MQ=255)
                                  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                      >  1:4447281/1‑70 (MQ=255)
                                  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                      >  1:866778/1‑70 (MQ=255)
                                  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                      >  1:3607389/1‑70 (MQ=255)
                                  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                      <  1:2321273/70‑1 (MQ=255)
                                   gAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGAGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTa                                    >  1:269440/1‑71 (MQ=255)
                                    aaGCGTGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTa                                    <  1:3936445/70‑1 (MQ=255)
                                     aGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACg                                  <  1:4309061/71‑1 (MQ=255)
                                      tccggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCt                                          <  1:1297308/59‑1 (MQ=255)
                                      tccggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCt                                          <  1:4354187/59‑1 (MQ=255)
                                      tccggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTACGa                                 >  1:276075/4‑71 (MQ=255)
                                      gCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                 >  1:4667558/1‑71 (MQ=255)
                                        cggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCtgc                                        >  1:5670470/2‑61 (MQ=255)
                                        cggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCtgc                                        >  1:3915422/2‑61 (MQ=255)
                                        cggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCtgc                                        >  1:1297300/2‑61 (MQ=255)
                                        cggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTACGac                                <  1:812112/69‑2 (MQ=255)
                                          gTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTAc                                   <  1:2662269/65‑1 (MQ=255)
                                           tCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTACGACGCCg                            <  1:5551407/71‑1 (MQ=255)
                                           tCAACCGACCGTGTGCTCAGAAACCTGCATGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTACGACGcc                             >  1:2781159/1‑70 (MQ=255)
                                            cAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGa                           >  1:1801949/1‑71 (MQ=255)
                                             aGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGc                                               <  1:2213539/50‑1 (MQ=255)
                                             aGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGc                                               <  1:3612705/50‑1 (MQ=255)
                                             aGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGAc                          >  1:4856313/1‑71 (MQ=255)
                                             aGCAGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGAc                          >  1:4834711/1‑71 (MQ=255)
                                              gCCGACCGTGTGATCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGcc                             <  1:1669589/67‑1 (MQ=255)
                                               ccGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTAc                                   <  1:1369153/60‑1 (MQ=255)
                                               ccGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTAt                                                     >  1:3189075/1‑42 (MQ=25)
                                                cGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGc                        <  1:3257921/70‑1 (MQ=255)
                                                cGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGc                        >  1:125299/1‑70 (MQ=255)
                                                    cGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                 >  1:1920537/1‑57 (MQ=255)
                                                       gtgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCt                                                   <  1:2584512/36‑1 (MQ=255)
                                                       gtgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCt                                                   <  1:5409495/36‑1 (MQ=255)
                                                       gtgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAAc                >  1:1712571/1‑71 (MQ=255)
                                                       gtgCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTAc                                   <  1:1753677/52‑1 (MQ=255)
                                                       gtgCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTAc                                   <  1:4624859/52‑1 (MQ=255)
                                                       gtgCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTAc                                   <  1:1364331/52‑1 (MQ=255)
                                                        tgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                 >  1:5142558/1‑53 (MQ=255)
                                                        tgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                 >  1:1093989/1‑53 (MQ=255)
                                                        tgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                 >  1:1306591/1‑53 (MQ=255)
                                                        tgCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTAc                                   <  1:1447620/51‑1 (MQ=255)
                                                          cTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTg             <  1:5285365/71‑1 (MQ=255)
                                                           tCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGt                                             <  1:100427/38‑1 (MQ=255)
                                                             aGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTa                                    >  1:5124531/1‑45 (MQ=255)
                                                             aGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGc            <  1:447159/69‑1 (MQ=255)
                                                                aaCCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGcagc         >  1:1542849/1‑69 (MQ=255)
                                                                 aCCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGcagccag      <  1:4551761/71‑1 (MQ=255)
                                                                 aCCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGcagccag      <  1:1057236/71‑1 (MQ=255)
                                                                   cTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGCCGCTATTGAACGTGCAGCCAGCa    >  1:3687299/1‑71 (MQ=255)
                                                                   cTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCa    >  1:5308503/1‑71 (MQ=255)
                                                                    tgtgTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCAc   >  1:516212/1‑71 (MQ=255)
                                                                     gtgtCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGcagc         >  1:1822240/1‑64 (MQ=255)
                                                                      tgtCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCAcc  >  1:1109144/1‑70 (MQ=255)
                                                                      |                                                                     
AGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACC  >  minE/836001‑836140

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: