Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,647 A→G 19.0% P1035P (CCA→CCG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,6470AG19.0% 165.2 / 33.5 79P1035P (CCA→CCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (36/28);  new base G (6/9);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

CCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCC  >  minE/834577‑834717
                                                                      |                                                                      
ccTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                                                                        <  1:997730/71‑1 (MQ=255)
    gtctggtctgGAAGATGAACACTTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                                                                        <  1:2066342/67‑1 (MQ=255)
    gtccggtctgGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTgg                                                                     <  1:3774248/70‑1 (MQ=255)
         gtctggAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACg                                                               >  1:578677/1‑71 (MQ=255)
         gtctggAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACg                                                               >  1:1977973/1‑71 (MQ=255)
          tctggAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACg                                                               <  1:3428785/70‑1 (MQ=255)
              gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctct                                                           <  1:887036/70‑1 (MQ=255)
              gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctct                                                           <  1:1248991/70‑1 (MQ=255)
              gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctct                                                           <  1:1625502/70‑1 (MQ=255)
                aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                                                        >  1:1127115/1‑71 (MQ=255)
                aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                                                        >  1:587615/1‑71 (MQ=255)
                aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                                                        >  1:3472485/1‑71 (MQ=255)
                aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                                                        >  1:3331147/1‑71 (MQ=255)
                aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                                                        >  1:1381597/1‑71 (MQ=255)
                 gATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTgg                                                       <  1:2667360/71‑1 (MQ=255)
                 gATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTgg                                                       <  1:893974/71‑1 (MQ=255)
                  aTGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                                                                        <  1:897873/53‑1 (MQ=255)
                   tGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtc                                                     >  1:2139387/1‑71 (MQ=255)
                    gAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATg                                                                      >  1:667060/1‑53 (MQ=255)
                     aaCACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtcgt                                                   <  1:552425/71‑1 (MQ=255)
                     aaCACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtcgt                                                   <  1:1873440/71‑1 (MQ=255)
                      tcacGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGctgt                                                           <  1:1475867/61‑3 (MQ=255)
                      tcacGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGgacgc                                                   >  1:1666397/2‑69 (MQ=255)
                      acacGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                                                                   >  1:2481734/1‑54 (MQ=255)
                      acacGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                                                                   >  1:1314246/1‑54 (MQ=255)
                       cacGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctc                                                            >  1:189805/1‑60 (MQ=255)
                       cacGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGctgtc                                                          >  1:826260/1‑62 (MQ=255)
                        acGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTa                                                                 >  1:1998141/1‑54 (MQ=255)
                        acGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTa                                                                 >  1:1628715/1‑54 (MQ=255)
                        acGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTa                                                                 >  1:2798653/1‑54 (MQ=255)
                        acGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGgacgc                                                   <  1:159881/68‑2 (MQ=255)
                        acGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAg                                                >  1:1962048/1‑71 (MQ=255)
                             tCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCa                                                 >  1:3565001/1‑65 (MQ=255)
                             tCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCa                                                 >  1:3464236/1‑65 (MQ=255)
                             tCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACt                                           >  1:1259746/1‑71 (MQ=255)
                                tataaCGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCa                                                 <  1:2323948/59‑1 (MQ=255)
                                 aCAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTa                                        <  1:2709250/70‑1 (MQ=255)
                                  taaCGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTa                                        <  1:3553339/68‑1 (MQ=255)
                                   aaCGCCGGTTACACCAACGTTCACGAGCTGATCCCAt                                                                       <  1:1747078/37‑1 (MQ=255)
                                   aaCGCCGGTTACACCAACGTTCACGAGCTGATCCCAt                                                                       <  1:1515248/37‑1 (MQ=255)
                                   aaCGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTAt                                       >  1:2042326/1‑69 (MQ=255)
                                    aCGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCa                                     <  1:965805/70‑1 (MQ=255)
                                    aCGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCAg                                    >  1:807224/1‑71 (MQ=255)
                                     cgcgGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGgacgc                                                   <  1:281607/55‑2 (MQ=255)
                                     cGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                             <  1:2240643/61‑1 (MQ=255)
                                     cGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAgg                                   <  1:2041313/71‑1 (MQ=255)
                                        cGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGAt                                 >  1:1650743/1‑70 (MQ=255)
                                         ggTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTAt                                       <  1:2732990/63‑1 (MQ=255)
                                         ggTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTAt                                       <  1:3104310/63‑1 (MQ=255)
                                          gTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctc                                                            <  1:2077042/41‑1 (MQ=255)
                                           ttACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAc                              >  1:1316801/1‑70 (MQ=255)
                                            tACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCAGGATCatcc                            <  1:3462352/71‑4 (MQ=255)
                                            tACACCAACGTTCACGAGCTTATCCCATGGCGTACGCTCTCTGGTCGTCa                                                 >  1:655461/1‑50 (MQ=255)
                                            tACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACt                                           >  1:2049253/1‑56 (MQ=255)
                                            tACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACt                                           >  1:3221082/1‑56 (MQ=255)
                                             acacCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCa                            <  1:3838051/70‑1 (MQ=255)
                                               acCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAcc                             <  1:7733/67‑1 (MQ=255)
                                               acCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTg                         <  1:3906863/71‑1 (MQ=255)
                                                 cAACGTTCACGAGCTGATCCCATGGCGTACGctctc                                                          >  1:2888050/1‑36 (MQ=255)
                                                 cAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgt                                                      >  1:1363833/1‑40 (MQ=255)
                                                 cAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGa                       >  1:1683362/1‑71 (MQ=255)
                                                 cAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGa                       >  1:107192/1‑71 (MQ=255)
                                                     gTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGc                    >  1:2720589/1‑70 (MQ=255)
                                                     gTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCg                   >  1:3904524/1‑71 (MQ=255)
                                                       tCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGt                  >  1:1241453/1‑70 (MQ=255)
                                                       tCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGt                  >  1:2369895/1‑70 (MQ=255)
                                                        cACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAg                           >  1:2554934/1‑60 (MQ=255)
                                                          cGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAtt              <  1:269066/71‑1 (MQ=255)
                                                          cGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAtt              <  1:2117088/71‑1 (MQ=255)
                                                          cGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAtt              <  1:202207/71‑1 (MQ=255)
                                                           gAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAttt             <  1:205397/71‑1 (MQ=255)
                                                            aGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCa                                              <  1:2546258/37‑1 (MQ=255)
                                                              cTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAttt             <  1:906536/68‑1 (MQ=255)
                                                                gATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCg           <  1:519447/68‑1 (MQ=255)
                                                                 aTCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTg                         <  1:2906606/53‑1 (MQ=255)
                                                                 aTCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGt         <  1:1306503/69‑1 (MQ=255)
                                                                 aTCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGa       <  1:3213661/71‑1 (MQ=255)
                                                                  tCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGa                       >  1:1831850/1‑54 (MQ=255)
                                                                  tCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGa                       >  1:1478975/1‑54 (MQ=255)
                                                                  tCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGa       >  1:3933976/1‑70 (MQ=255)
                                                                    ccATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAg    >  1:3103255/1‑71 (MQ=255)
                                                                     cATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGc   >  1:3604711/1‑71 (MQ=255)
                                                                      aTGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGcc  >  1:1294340/1‑71 (MQ=255)
                                                                      |                                                                      
CCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCC  >  minE/834577‑834717

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: