Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,007 T→A 15.9% V1155V (GTT→GTA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,0070TA15.9% 236.9 / 32.0 94V1155V (GTT→GTAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (29/50);  new base A (7/8);  total (36/58)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACC  >  minE/834938‑835076
                                                                     |                                                                     
tAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTc                                                                      <  1:1020665/71‑1 (MQ=255)
tAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTc                                                                      <  1:2633304/71‑1 (MQ=255)
  aCGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTcc                                                                     <  1:1988737/70‑1 (MQ=255)
   cGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCgg                                                                   <  1:422534/71‑1 (MQ=255)
            tGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGa                                                          >  1:120265/1‑71 (MQ=255)
              aaGTCTTCAACAGCAACGGTGCACTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGc                                                                  <  1:2488119/61‑1 (MQ=255)
               aGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCgg                                                                   <  1:2297830/59‑1 (MQ=255)
               aGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCa                                                       <  1:992511/71‑1 (MQ=255)
               aGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCa                                                       <  1:2098452/71‑1 (MQ=255)
               aGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCa                                                       <  1:3534834/71‑1 (MQ=255)
               aGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCa                                                       <  1:3099332/71‑1 (MQ=255)
                gtattCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCat                                                      <  1:2917968/68‑1 (MQ=255)
                gTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCat                                                      <  1:1371996/71‑1 (MQ=255)
                gTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCat                                                      >  1:3940907/1‑71 (MQ=255)
                   ttCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTcc                                                                     >  1:2716712/1‑53 (MQ=255)
                   ttCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCatgat                                                   >  1:2106827/1‑71 (MQ=255)
                   ttCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCat                                                      <  1:1572509/68‑1 (MQ=255)
                   ttCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCat                                                      <  1:1945452/68‑1 (MQ=255)
                   atCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCatga                                                    <  1:755669/69‑1 (MQ=255)
                    tCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCatgatg                                                  >  1:2692629/1‑71 (MQ=255)
                     cAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCatgatg                                                  <  1:2265125/70‑1 (MQ=255)
                     cAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCatgatg                                                  <  1:2167323/70‑1 (MQ=255)
                     cAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCatgatg                                                  <  1:1467150/70‑1 (MQ=255)
                     cAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCatga                                                    <  1:1150338/68‑1 (MQ=255)
                       aCAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTAc                                               <  1:2522119/71‑1 (MQ=255)
                        cAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTAcc                                              <  1:465427/71‑1 (MQ=255)
                        cAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTAc                                               >  1:1048619/1‑70 (MQ=255)
                        cAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTAc                                               >  1:3352417/1‑70 (MQ=255)
                           cAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCAc                                            <  1:1064500/70‑1 (MQ=255)
                            aaCGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCAcgc                                          >  1:2375329/1‑71 (MQ=255)
                            aaCGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCAcgc                                          <  1:2197837/71‑1 (MQ=255)
                             aCGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCAcgcg                                         <  1:2548474/71‑1 (MQ=255)
                             aCGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCAcgcg                                         <  1:2606534/71‑1 (MQ=255)
                              cGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCAcgcg                                         <  1:999822/70‑1 (MQ=255)
                              cGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCAcgcc                                         >  1:2651998/1‑69 (MQ=255)
                                 tGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAgg                                     <  1:1451599/71‑1 (MQ=255)
                                  gCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAgg                                     >  1:2691125/1‑70 (MQ=255)
                                   ctctGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCTCAGGa                                    <  1:2769790/70‑1 (MQ=255)
                                   ctctGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGa                                    <  1:3419235/70‑1 (MQ=255)
                                     ctGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACg                                 >  1:2245694/1‑71 (MQ=255)
                                     ctGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCatg                                                     <  1:1404635/51‑1 (MQ=255)
                                      tGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACg                                 >  1:2452496/1‑70 (MQ=255)
                                        aCTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTa                               <  1:3343759/70‑1 (MQ=255)
                                        aCTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTAt                              <  1:2340125/71‑1 (MQ=255)
                                        aCTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTAt                              <  1:3157499/71‑1 (MQ=255)
                                        aCTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCAt                              >  1:2705619/1‑71 (MQ=255)
                                         cTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAgg                                                               >  1:3412902/1‑37 (MQ=255)
                                         cTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATc                             <  1:3439265/71‑1 (MQ=255)
                                         cTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCAtt                             <  1:2391046/71‑2 (MQ=255)
                                         cTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCAtt                             <  1:1094271/71‑2 (MQ=255)
                                         cTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCAt                              >  1:1391483/1‑70 (MQ=255)
                                          tGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCg                            >  1:3326240/1‑71 (MQ=255)
                                          tGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCg                            >  1:1554806/1‑71 (MQ=255)
                                            cccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                          <  1:3282058/71‑1 (MQ=255)
                                            cccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                          <  1:3189789/71‑1 (MQ=255)
                                            cccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                          <  1:3185670/71‑1 (MQ=255)
                                             ccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                         <  1:2400344/71‑1 (MQ=255)
                                             ccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                         <  1:1291355/71‑1 (MQ=255)
                                              cGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCatg                                                     <  1:3121471/42‑1 (MQ=255)
                                               gTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAc                       >  1:243385/1‑71 (MQ=255)
                                                tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCatgatg                                                  <  1:2912819/43‑1 (MQ=255)
                                                tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTa                                                <  1:71734/45‑1 (MQ=255)
                                                tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCAcgcg                                         <  1:1470713/52‑1 (MQ=255)
                                                tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcc                      >  1:218630/1‑71 (MQ=255)
                                                tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcc                      >  1:3059265/1‑71 (MQ=255)
                                                cgCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGaa                                   >  1:3400297/2‑58 (MQ=255)
                                                  cGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctg                    <  1:104930/71‑1 (MQ=255)
                                                  cGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctg                    >  1:3925061/1‑71 (MQ=255)
                                                     ttGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCAcgcg                                         <  1:1695276/47‑1 (MQ=255)
                                                     ttGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGt                                <  1:2749344/56‑1 (MQ=255)
                                                     ttGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGt                                <  1:1402295/56‑1 (MQ=255)
                                                     ttGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTa                               >  1:1742048/1‑57 (MQ=255)
                                                     ttGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctgcc                  <  1:456491/70‑1 (MQ=255)
                                                       gTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCa                                       >  1:2524787/1‑47 (MQ=255)
                                                        tCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGt              <  1:2091218/71‑1 (MQ=255)
                                                        tCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACgca                               <  1:3557089/54‑3 (MQ=16)
                                                          agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc            >  1:1489476/1‑71 (MQ=255)
                                                          agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc            >  1:3273817/1‑71 (MQ=255)
                                                          agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc            >  1:3278299/1‑71 (MQ=255)
                                                          agccaacGTGTACCGCCGGGCATGACCATGATGTATCAcgccc                                        >  1:50181/1‑41 (MQ=12)
                                                           ggcagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                         <  1:1995826/55‑1 (MQ=255)
                                                           gccaacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACgca                               <  1:1784727/51‑3 (MQ=16)
                                                            ccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCgg          <  1:835364/71‑1 (MQ=255)
                                                            ccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCg           >  1:525714/1‑70 (MQ=255)
                                                             cagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa         <  1:3403729/71‑1 (MQ=255)
                                                             cagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa         <  1:3941251/71‑1 (MQ=255)
                                                              agcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAAc                                  >  1:2446594/1‑45 (MQ=255)
                                                              agcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCGTGGTTCGGa         >  1:3309073/1‑70 (MQ=255)
                                                              agcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa         >  1:2686197/1‑70 (MQ=255)
                                                               gcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGaaa       >  1:1025423/1‑71 (MQ=255)
                                                               gcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGaaa       >  1:452458/1‑71 (MQ=255)
                                                               gcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGaa        >  1:2592330/1‑70 (MQ=255)
                                                               acGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAAc                                  >  1:341739/2‑44 (MQ=12)
                                                                  tgtTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCg                            <  1:1959701/47‑1 (MQ=255)
                                                                    ttCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAcc  <  1:3672553/71‑1 (MQ=255)
                                                                     tCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAcc  >  1:1647102/1‑70 (MQ=255)
                                                                     |                                                                     
TAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACC  >  minE/834938‑835076

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: