Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,635 C→T 20.5% H1031H (CAC→CAT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,6350CT20.5% 152.5 / 32.6 73H1031H (CAC→CATnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (35/23);  new base T (6/9);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATT  >  minE/834569‑834705
                                                                  |                                                                      
cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGa                                                                      <  1:3493738/69‑1 (MQ=255)
cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGc                                                                    >  1:1885374/1‑71 (MQ=255)
    ccGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGAt                                                                >  1:2931444/1‑71 (MQ=255)
     cGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATc                                                               >  1:2420899/1‑71 (MQ=255)
        ccTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                                                            <  1:997730/71‑1 (MQ=255)
            gtctggtctgGAAGATGAACACTTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                                                            <  1:2066342/67‑1 (MQ=255)
            gtccggtctgGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTgg                                                         <  1:3774248/70‑1 (MQ=255)
                 gtctggAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACg                                                   >  1:1977973/1‑71 (MQ=255)
                 gtctggAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACg                                                   >  1:578677/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:1381597/1‑71 (MQ=255)
                        aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                                            >  1:3472485/1‑71 (MQ=255)
                        aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                                            >  1:587615/1‑71 (MQ=255)
                        aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                                            >  1:3331147/1‑71 (MQ=255)
                         gATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTgg                                           <  1:893974/71‑1 (MQ=255)
                         gATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTgg                                           <  1:2667360/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:1873440/71‑1 (MQ=255)
                             aaCACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtcgt                                       <  1:552425/71‑1 (MQ=255)
                              tcacGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGctgt                                               <  1:1475867/61‑3 (MQ=255)
                              tcacGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGgacgc                                       >  1:1666397/2‑69 (MQ=255)
                              acacGTTTCTTACAACGCCGGTTACACCAACGTTCACGAg                                                                     >  1:1072750/1‑40 (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:2798653/1‑54 (MQ=255)
                                acGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTa                                                     >  1:1628715/1‑54 (MQ=255)
                                acGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGgacgc                                       <  1:159881/68‑2 (MQ=255)
                                acGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAg                                    >  1:1962048/1‑71 (MQ=255)
                                 cGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTg                                                                  <  1:1888435/40‑1 (MQ=255)
                                     tCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCa                                     >  1:3464236/1‑65 (MQ=255)
                                     tCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCa                                     >  1:3565001/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:3104310/63‑1 (MQ=255)
                                                 ggTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTAt                           <  1:2732990/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:2369895/1‑70 (MQ=255)
                                                               tCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGt      >  1:1241453/1‑70 (MQ=255)
                                                                cACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAg               >  1:2554934/1‑60 (MQ=255)
                                                                  cGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAtt  <  1:269066/71‑1 (MQ=255)
                                                                  cGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAtt  <  1:202207/71‑1 (MQ=255)
                                                                  cGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAtt  <  1:2117088/71‑1 (MQ=255)
                                                                  |                                                                      
CAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATT  >  minE/834569‑834705

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: