Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,249 T→G 25.9% G569G (GGT→GGG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,2490TG25.9% 102.0 / 62.6 85G569G (GGT→GGGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (34/29);  new base G (16/6);  total (50/35)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

ACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGC  >  minE/833179‑833319
                                                                      |                                                                      
aCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                                                                        <  1:1917723/71‑1 (MQ=255)
aCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                                                                        <  1:575325/71‑1 (MQ=255)
 ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                                                                        >  1:3930466/1‑70 (MQ=255)
 ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTc                                                                       <  1:1400640/71‑1 (MQ=255)
 ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGGGCTGTGTCGGt                                                                        >  1:5052111/1‑69 (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)
            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:2477942/48‑1 (MQ=255)
                          tATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                                     <  1:2191825/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: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:3851377/1‑58 (MQ=255)
                                       tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                                              >  1:3857956/1‑58 (MQ=255)
                                       tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                                              >  1:49515/1‑58 (MQ=255)
                                       tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                                              >  1:2691954/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:188112/1‑58 (MQ=255)
                                            aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAg                                         >  1:1262426/1‑58 (MQ=255)
                                            aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAg                                         >  1:4690972/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)
                                                  atcttCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa                              >  1:3438767/4‑63 (MQ=255)
                                                  aTTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                                        <  1:5604210/37‑1 (MQ=255)
                                                  aTTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                                        <  1:3063288/37‑1 (MQ=255)
                                                   ttttCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                                        <  1:4352605/36‑1 (MQ=255)
                                                   ttttCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                                        <  1:2520814/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)
                                                               gtgtCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTgg         >  1:5529711/1‑71 (MQ=255)
                                                                tgtCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTgg         <  1:1003885/70‑1 (MQ=255)
                                                                tgtCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTggt        <  1:4121250/71‑2 (MQ=255)
                                                                tgtCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTgg         <  1:1119725/70‑1 (MQ=255)
                                                                  tCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                             >  1:1166919/1‑48 (MQ=255)
                                                                  tCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCaa                  <  1:2999885/59‑1 (MQ=255)
                                                                  tCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCaa                  <  1:1294849/59‑1 (MQ=255)
                                                                   cGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGc     >  1:5440272/1‑71 (MQ=255)
                                                                    ggTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                             >  1:3649610/1‑46 (MQ=255)
                                                                    ggTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                             >  1:3089178/1‑46 (MQ=255)
                                                                      tCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGcccc  <  1:3377181/71‑3 (MQ=255)
                                                                      |                                                                      
ACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGC  >  minE/833179‑833319

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: