Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,464 C→A 7.3% L308M (CTG→ATG)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,4640CA 7.3% 259.0 / 4.5 81L308M (CTG→ATG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (21/54);  new base A (3/3);  total (24/57)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

TACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAT  >  minE/832395‑832532
                                                                     |                                                                    
tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                                                                     >  1:1171134/1‑71 (MQ=255)
  cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg                                                                   >  1:4875849/1‑71 (MQ=255)
  cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg                                                                   <  1:2331287/71‑1 (MQ=255)
   gCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGc                                                                  >  1:4841478/1‑71 (MQ=255)
    cTGAATTCGCCAAAGTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCg                                                                 <  1:2510370/71‑1 (MQ=255)
     tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                                                                     <  1:471326/66‑1 (MQ=255)
     tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa                                                                >  1:4293305/1‑71 (MQ=255)
       aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa                                                                <  1:30606/69‑1 (MQ=255)
       aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                               <  1:2656352/70‑1 (MQ=255)
       aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATg                                                              <  1:4550437/71‑1 (MQ=255)
       aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATg                                                              >  1:4135535/1‑71 (MQ=255)
         aTCGCCAAACTGTGCGATCTGTGGCTGGCCCCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggt                                                            >  1:1041489/1‑70 (MQ=255)
         aTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggg                                                            >  1:1459041/1‑71 (MQ=255)
            gCCAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCa                                                         >  1:5162204/1‑71 (MQ=255)
             ccAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCat                                                        <  1:4956967/71‑2 (MQ=255)
              cAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCatg                                                       >  1:972398/1‑69 (MQ=255)
               aaaTTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGt                                                      <  1:4691403/71‑1 (MQ=255)
               aaaTTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGt                                                      <  1:3292348/71‑1 (MQ=255)
               aaaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                                      <  1:647448/71‑1 (MQ=255)
                       cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                                                    >  1:4667873/1‑49 (MQ=255)
                       cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                                                    >  1:1306486/1‑49 (MQ=255)
                       cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                                                    >  1:1951876/1‑49 (MQ=255)
                       cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                                                    >  1:711005/1‑49 (MQ=255)
                       cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGc                                                                  <  1:3870600/51‑1 (MQ=255)
                       cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGc                                                                  <  1:950805/51‑1 (MQ=255)
                          tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGt                                            >  1:1595232/1‑70 (MQ=255)
                           cTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa                                                                >  1:5161605/1‑49 (MQ=255)
                              tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                           <  1:4393558/67‑1 (MQ=255)
                              tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa                                          <  1:3415584/68‑1 (MQ=255)
                               ggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                           <  1:5047364/66‑1 (MQ=255)
                                gctggcACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATg                                                              >  1:3183114/1‑46 (MQ=255)
                                     cACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg                                                                   <  1:5106141/36‑1 (MQ=255)
                                      aCCGAAACAGGGCTCCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa                                          <  1:545220/60‑1 (MQ=255)
                                       ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                                      <  1:689026/47‑1 (MQ=255)
                                       ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                                      <  1:2098150/47‑1 (MQ=255)
                                       ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                                      <  1:3193421/47‑1 (MQ=255)
                                       ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc                                <  1:2653088/69‑1 (MQ=255)
                                        cGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGa                              >  1:991128/1‑70 (MQ=255)
                                        cGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                             <  1:2245183/71‑1 (MQ=255)
                                         gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtg                                                  >  1:1505257/1‑49 (MQ=255)
                                         gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                             >  1:2577609/1‑70 (MQ=255)
                                           aaCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACaa                           >  1:5637108/1‑70 (MQ=255)
                                            aCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa                                          <  1:4273547/54‑1 (MQ=255)
                                             cAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                             <  1:2510651/66‑1 (MQ=255)
                                             cAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAAccc                        <  1:3708554/71‑1 (MQ=255)
                                              aGGGCACCGATGCGTCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                             <  1:3677352/65‑1 (MQ=255)
                                                ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc                                              <  1:4930043/46‑1 (MQ=255)
                                                ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc                                              <  1:3481310/46‑1 (MQ=255)
                                                ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc                                              <  1:2955774/46‑1 (MQ=255)
                                                ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc                                              <  1:5364858/46‑1 (MQ=255)
                                                 gCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc                    <  1:903663/71‑1 (MQ=255)
                                                  cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                           <  1:3710010/47‑1 (MQ=255)
                                                  cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                           <  1:3703196/47‑1 (MQ=255)
                                                  cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                           <  1:249938/47‑1 (MQ=255)
                                                    ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                           <  1:116302/45‑1 (MQ=255)
                                                    ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                           <  1:1822337/45‑1 (MQ=255)
                                                    ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg                 <  1:1412590/71‑1 (MQ=255)
                                                    ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg                 <  1:860809/71‑1 (MQ=255)
                                                    ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGATGCGTGAATTCCACCTCGACAACCCAAGCCAg                 <  1:1297774/71‑1 (MQ=255)
                                                     cGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg                 >  1:2784953/1‑70 (MQ=255)
                                                     cGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGt                <  1:4327533/71‑1 (MQ=255)
                                                     cGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGt                <  1:5093007/71‑1 (MQ=255)
                                                      gATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc                                <  1:4033357/54‑1 (MQ=255)
                                                      gATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg                 <  1:5112555/69‑1 (MQ=255)
                                                       aTGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt              <  1:481577/71‑1 (MQ=255)
                                                       aTGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt              <  1:1572162/71‑1 (MQ=255)
                                                       aTGCGGCAATGGCGCTGGCGATGGGCCACCTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt              <  1:2553456/71‑1 (MQ=255)
                                                           ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc                                <  1:4976215/49‑1 (MQ=255)
                                                           ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc                                <  1:2450968/49‑1 (MQ=255)
                                                           ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCa          <  1:4598708/71‑1 (MQ=255)
                                                           ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCa          <  1:400012/71‑1 (MQ=255)
                                                             cAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAttt            >  1:4337798/1‑67 (MQ=255)
                                                             cAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCAAGTAttt            >  1:4649254/1‑67 (MQ=255)
                                                              aaTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                             <  1:868249/49‑1 (MQ=255)
                                                               aTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCaa                      <  1:3529823/55‑1 (MQ=255)
                                                               aTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                   <  1:805254/58‑1 (MQ=255)
                                                               aTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGa      <  1:1758198/71‑1 (MQ=255)
                                                                tGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACa                            <  1:2770755/48‑1 (MQ=255)
                                                                tGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACa                            <  1:1772831/48‑1 (MQ=255)
                                                                tGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACa                            <  1:160473/48‑1 (MQ=255)
                                                                tGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCTAc     >  1:4360562/1‑71 (MQ=255)
                                                                  gcgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACt    <  1:5683136/70‑1 (MQ=255)
                                                                  gcgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTa   <  1:2298415/71‑1 (MQ=255)
                                                                  gcgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTa   <  1:1635669/71‑1 (MQ=255)
                                                                  gcgcTGGCGATGGGCCAAGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCAc         >  1:4827672/1‑65 (MQ=255)
                                                                   cgcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAt  <  1:5069339/71‑1 (MQ=255)
                                                                     |                                                                    
TACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAT  >  minE/832395‑832532

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: