Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,231 T→C 24.2% I563I (ATT→ATC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,2310TC24.2% 154.2 / 50.9 91I563I (ATT→ATCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (35/34);  new base C (15/7);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

TCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCG  >  minE/833161‑833300
                                                                      |                                                                     
tcgtcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGAtt                                                                       >  1:5324934/1‑71 (MQ=255)
        gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGAtttt                                                                     >  1:2960422/1‑65 (MQ=255)
        gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGAtttt                                                                     >  1:4236686/1‑65 (MQ=255)
        gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCgg                                                               >  1:5051978/1‑71 (MQ=255)
        gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCg                                                                >  1:4316618/1‑70 (MQ=255)
        gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCg                                                                >  1:5672695/1‑70 (MQ=255)
           gggCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTg                                                                  <  1:1363383/65‑1 (MQ=255)
           gggCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTg                                                                  <  1:5569653/65‑1 (MQ=255)
            ggCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCtgt                                                           <  1:1035629/71‑1 (MQ=255)
             gCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCtgtg                                                          >  1:3725412/1‑71 (MQ=255)
               ttaaCCACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCTGATCTTCTGCGGTTgtgt                                                         <  1:3337644/68‑1 (MQ=255)
                gAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCg                                                       >  1:3202465/1‑71 (MQ=255)
                 aaCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                                                    <  1:4707401/57‑1 (MQ=255)
                 aaCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                                                    <  1:4140865/57‑1 (MQ=255)
                 aaCCACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCTGATCTTCTGCGGTTGTGTCgg                                                      >  1:2572011/1‑71 (MQ=255)
                  aCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                                                     <  1:1917723/71‑1 (MQ=255)
                  aCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                                                     <  1:575325/71‑1 (MQ=255)
                  aCCACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCTGATCTTCTGCGGTTGTGTCggg                                                     >  1:3094736/1‑70 (MQ=255)
                   ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTg                                                                  >  1:4676367/1‑57 (MQ=255)
                   ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTg                                                                  >  1:2864839/1‑57 (MQ=255)
                   ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                                                     >  1:3930466/1‑70 (MQ=255)
                   ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTc                                                    <  1:1400640/71‑1 (MQ=255)
                   ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGGGCTGTGTCGGt                                                     >  1:5052111/1‑69 (MQ=255)
                   ccACTGGTATCACCTCGATATGAACTATCGGGGTCTGATCAACATGCTGATTTTCTg                                                                  >  1:5203024/1‑57 (MQ=255)
                   ccACTGGTATCACCTCGAAATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTg                                                                  >  1:3727123/1‑57 (MQ=255)
                   ccACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCTGATCTTCTGCGGTTGTGTCgggc                                                    <  1:4169593/71‑3 (MQ=255)
                   ccACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCTGATCTTCTGCGGTTGTGTCgggc                                                    <  1:1205448/71‑3 (MQ=255)
                    cACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCa                                                   <  1:1693458/71‑1 (MQ=255)
                     aCTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                  <  1:3824254/71‑1 (MQ=255)
                     aCTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCTGATCTTCTGCGGTTGTGTCGGGCa                                                   >  1:5614519/1‑70 (MQ=255)
                         gTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCagag                                                <  1:763299/69‑1 (MQ=255)
                            tCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                                                    <  1:4202694/46‑1 (MQ=255)
                            tCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                                                    <  1:5452341/46‑1 (MQ=255)
                            tCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGc                                               <  1:3514023/67‑1 (MQ=255)
                            tCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCCGCTGTGTCGGTCAGAGc                                               <  1:577276/67‑1 (MQ=255)
                            tCACCTCGATATGAACTATCGTGATCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGc                                               <  1:3016432/67‑1 (MQ=255)
                             cACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGc                                                              <  1:1675534/51‑1 (MQ=255)
                             cACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGc                                                              <  1:5107298/51‑1 (MQ=255)
                              aCCTCGATATGAACTATCGTGGTCTGATCAACATGCTGAtttt                                                                     >  1:2173763/1‑43 (MQ=255)
                              aCCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTg                                                                  >  1:5588664/1‑46 (MQ=255)
                              aCCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCgg                                             >  1:1604126/1‑67 (MQ=255)
                                  cGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                     >  1:832257/1‑71 (MQ=255)
                                   gATATGAACTATCGTGGTCTGATCAACATGCTGAtttt                                                                     >  1:3695424/1‑38 (MQ=255)
                                   gATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                                                    >  1:5221555/1‑39 (MQ=255)
                                   gATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                                                    >  1:884748/1‑39 (MQ=255)
                                   gATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTg                                    >  1:1200263/1‑71 (MQ=255)
                                    atatGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCa                                                   <  1:267822/55‑1 (MQ=255)
                                    atatGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCgg                                             >  1:3879120/1‑61 (MQ=255)
                                    atatGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGAgg                                             >  1:3948894/1‑61 (MQ=255)
                                     tatGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTggg                                  <  1:4653097/71‑1 (MQ=255)
                                       tGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                  <  1:2083719/53‑1 (MQ=255)
                                        gAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCggg                                            >  1:47930/1‑58 (MQ=255)
                                          aCTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCAc                             >  1:3275974/1‑71 (MQ=255)
                                           cTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGc                                               >  1:3281777/1‑52 (MQ=255)
                                            tATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                  <  1:2191825/48‑1 (MQ=255)
                                            tATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                  <  1:2477942/48‑1 (MQ=255)
                                            tATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                  <  1:2566668/48‑1 (MQ=255)
                                             aTCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTAt                          >  1:1231917/1‑71 (MQ=255)
                                              tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTa                           <  1:2326530/69‑1 (MQ=255)
                                              tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                         >  1:407809/1‑71 (MQ=255)
                                              tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                         >  1:2480508/1‑71 (MQ=255)
                                              tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGAGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                         <  1:472933/71‑1 (MQ=255)
                                                   gTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                    <  1:4566378/71‑1 (MQ=255)
                                                   gTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                    <  1:1916905/71‑1 (MQ=255)
                                                    tCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                    >  1:806225/1‑70 (MQ=255)
                                                     cTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCa                                                   >  1:3725610/1‑38 (MQ=255)
                                                     cTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCa                                                   >  1:1198842/1‑38 (MQ=255)
                                                     cTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCa                                                   >  1:4474075/1‑38 (MQ=255)
                                                      tGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                  <  1:1527186/38‑1 (MQ=255)
                                                        aTCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAg                 >  1:2821760/1‑69 (MQ=255)
                                                        aTCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGa               <  1:2944028/71‑1 (MQ=255)
                                                        aTCAACATGCTGATTTTCTGCGGCTGTGTCGGTCagag                                                <  1:2509602/38‑1 (MQ=255)
                                                         tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                           >  1:3857956/1‑58 (MQ=255)
                                                         tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                           >  1:49515/1‑58 (MQ=255)
                                                         tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                           >  1:5152176/1‑58 (MQ=255)
                                                         tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                           >  1:2691954/1‑58 (MQ=255)
                                                         tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                           >  1:3851377/1‑58 (MQ=255)
                                                          cAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGgc                                 >  1:5581015/1‑51 (MQ=255)
                                                           aataTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaa            <  1:4184778/68‑1 (MQ=255)
                                                           aataTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaa            <  1:3016967/68‑1 (MQ=255)
                                                           aataTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaa             >  1:2981445/4‑70 (MQ=255)
                                                           aaCATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaa             <  1:3493599/70‑1 (MQ=255)
                                                           aaCATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaa             <  1:2167987/70‑1 (MQ=255)
                                                              aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAg                      >  1:4690972/1‑58 (MQ=255)
                                                              aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAg                      >  1:188112/1‑58 (MQ=255)
                                                              aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAg                      >  1:1262426/1‑58 (MQ=255)
                                                              aTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa           >  1:3406557/1‑69 (MQ=255)
                                                              aTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa           >  1:4464413/1‑69 (MQ=255)
                                                              aTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa           >  1:4304371/1‑69 (MQ=255)
                                                               tGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCggggg                                          >  1:1931444/1‑37 (MQ=255)
                                                               tGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCggggg                                          >  1:1402031/1‑37 (MQ=255)
                                                               tGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACt                            <  1:5102118/51‑1 (MQ=255)
                                                               tGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa           >  1:4945892/1‑68 (MQ=255)
                                                               tGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTg        >  1:3034003/1‑71 (MQ=255)
                                                                    aTTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                     <  1:5604210/37‑1 (MQ=255)
                                                                    aTTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                     <  1:3063288/37‑1 (MQ=255)
                                                                     ttttCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                     <  1:2520814/36‑1 (MQ=255)
                                                                     ttttCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                     <  1:4352605/36‑1 (MQ=255)
                                                                      tttCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTgg                                   >  1:2816984/1‑37 (MQ=255)
                                                                      tttCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCg  >  1:1400618/1‑70 (MQ=255)
                                                                      |                                                                     
TCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCG  >  minE/833161‑833300

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: