Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,240 C→T 27.6% G566G (GGC→GGT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,2400CT27.6% 132.3 / 55.2 87G566G (GGC→GGTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base T (17/7);  total (51/36)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

GGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCA  >  minE/833173‑833302
                                                                   |                                                              
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)
     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)
       ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                                                       >  1:3930466/1‑70 (MQ=255)
       ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTc                                                      <  1:1400640/71‑1 (MQ=255)
       ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGGGCTGTGTCGGt                                                       >  1:5052111/1‑69 (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)
                tCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGc                                                 <  1:3514023/67‑1 (MQ=255)
                tCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCCGCTGTGTCGGTCAGAGc                                                 <  1:577276/67‑1 (MQ=255)
                tCACCTCGATATGAACTATCGTGATCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGc                                                 <  1:3016432/67‑1 (MQ=255)
                 cACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGc                                                                <  1:5107298/51‑1 (MQ=255)
                 cACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGc                                                                <  1:1675534/51‑1 (MQ=255)
                  aCCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCgg                                               >  1:1604126/1‑67 (MQ=255)
                      cGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                       >  1:832257/1‑71 (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:2566668/48‑1 (MQ=255)
                                tATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                    <  1:2477942/48‑1 (MQ=255)
                                 aTCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTAt                            >  1:1231917/1‑71 (MQ=255)
                                  tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTa                             <  1:2326530/69‑1 (MQ=255)
                                  tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                           >  1:2480508/1‑71 (MQ=255)
                                  tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                           >  1:407809/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:4474075/1‑38 (MQ=255)
                                         cTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCa                                                     >  1:1198842/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:5152176/1‑58 (MQ=255)
                                             tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                             >  1:49515/1‑58 (MQ=255)
                                             tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                             >  1:2691954/1‑58 (MQ=255)
                                             tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                             >  1:3851377/1‑58 (MQ=255)
                                             tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                             >  1:3857956/1‑58 (MQ=255)
                                              cAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGgc                                   >  1:5581015/1‑51 (MQ=255)
                                               aataTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaa              <  1:3016967/68‑1 (MQ=255)
                                               aataTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaa              <  1:4184778/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:1262426/1‑58 (MQ=255)
                                                  aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAg                        >  1:4690972/1‑58 (MQ=255)
                                                  aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAg                        >  1:188112/1‑58 (MQ=255)
                                                  aTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa             >  1:4304371/1‑69 (MQ=255)
                                                  aTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa             >  1:4464413/1‑69 (MQ=255)
                                                  aTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa             >  1:3406557/1‑69 (MQ=255)
                                                   tGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCggggg                                            >  1:1402031/1‑37 (MQ=255)
                                                   tGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCggggg                                            >  1:1931444/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)
                                                        atcttCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa             >  1:3438767/4‑63 (MQ=255)
                                                        aTTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                       <  1:3063288/37‑1 (MQ=255)
                                                        aTTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                       <  1:5604210/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)
                                                          cttCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGccca    >  1:949659/2‑69 (MQ=255)
                                                            tCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCa  <  1:5222416/70‑1 (MQ=255)
                                                              tGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCa                                >  1:1644431/1‑38 (MQ=255)
                                                              tGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc            >  1:2936667/1‑58 (MQ=255)
                                                              tGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc            >  1:4963301/1‑58 (MQ=255)
                                                               gCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaa               >  1:3734738/1‑54 (MQ=255)
                                                                 ggCTGTTTCGGTCAGAGCGGGGGCGGCTGGGCGCTCTATGTa                         >  1:2750740/1‑42 (MQ=38)
                                                                 ggCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                      >  1:3400123/1‑45 (MQ=255)
                                                                 ggCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                      >  1:3760021/1‑45 (MQ=255)
                                                                   |                                                              
GGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCA  >  minE/833173‑833302

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: