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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,4580AC 6.7% 236.1 / 9.5 90M306L (ATG→CTG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (33/51);  new base C (3/3);  total (36/54)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

ACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGA  >  minE/832388‑832528
                                                                      |                                                                      
acCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCaa                                                                        <  1:3964059/71‑1 (MQ=255)
  cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                                                      >  1:5291713/1‑71 (MQ=255)
  cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                                                      >  1:4531944/1‑71 (MQ=255)
   aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                                     >  1:1044184/1‑71 (MQ=255)
   aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                                     <  1:4082993/71‑1 (MQ=255)
   aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                                     <  1:3169226/71‑1 (MQ=255)
    gACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                                                                    <  1:2195511/71‑1 (MQ=255)
     aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                                                                   >  1:2937726/1‑71 (MQ=255)
     aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                                                                   >  1:4369683/1‑71 (MQ=255)
     aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                                                                   >  1:3295111/1‑71 (MQ=255)
     aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                                                                   >  1:3504627/1‑71 (MQ=255)
     aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                                                                    >  1:5640309/1‑70 (MQ=255)
       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)
            tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                                     <  1:2759535/62‑1 (MQ=255)
            tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                                                                 <  1:471326/66‑1 (MQ=255)
            tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa                                                            >  1:4293305/1‑71 (MQ=255)
              aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                                                                   >  1:4746651/1‑62 (MQ=255)
              aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa                                                            <  1:30606/69‑1 (MQ=255)
              aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                           <  1:2656352/70‑1 (MQ=255)
              aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATg                                                          >  1:4135535/1‑71 (MQ=255)
              aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATg                                                          <  1:4550437/71‑1 (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)
                         cTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                                     >  1:726344/1‑49 (MQ=255)
                         cTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                                     >  1:234138/1‑49 (MQ=255)
                              cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                                                >  1:1951876/1‑49 (MQ=255)
                              cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                                                >  1:1306486/1‑49 (MQ=255)
                              cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                                                >  1:4667873/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)
                                  cTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                                                                    >  1:2328187/1‑41 (MQ=255)
                                  cTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                                                                    >  1:4393556/1‑41 (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:3193421/47‑1 (MQ=255)
                                              ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                                  <  1:2098150/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:3481310/46‑1 (MQ=255)
                                                       ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc                                          <  1:5364858/46‑1 (MQ=255)
                                                       ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc                                          <  1:4930043/46‑1 (MQ=255)
                                                       ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc                                          <  1:2955774/46‑1 (MQ=255)
                                                        gCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc                <  1:903663/71‑1 (MQ=255)
                                                         cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                       <  1:249938/47‑1 (MQ=255)
                                                         cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                       <  1:3710010/47‑1 (MQ=255)
                                                         cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                       <  1:3703196/47‑1 (MQ=255)
                                                           ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                       <  1:1822337/45‑1 (MQ=255)
                                                           ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                       <  1:116302/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:5093007/71‑1 (MQ=255)
                                                            cGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGt            <  1:4327533/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:2450968/49‑1 (MQ=255)
                                                                  ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc                            <  1:4976215/49‑1 (MQ=255)
                                                                  ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCa      <  1:400012/71‑1 (MQ=255)
                                                                  ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCa      <  1:4598708/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)
                                                                      |                                                                      
ACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGA  >  minE/832388‑832528

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: