Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,722 T→G 8.3% G65G (GGT→GGG narI → nitrate reductase 1, gamma (cytochrome b(NR)) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE837,7220TG 8.3% 116.1 / 7.0 48G65G (GGT→GGGnarInitrate reductase 1, gamma (cytochrome b(NR)) subunit
Reads supporting (aligned to +/- strand):  ref base T (22/22);  new base G (2/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

CTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAA  >  minE/837654‑837792
                                                                    |                                                                      
cTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTc                                                                       <  1:1231715/70‑1 (MQ=255)
      cGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCActtctt                                                                >  1:1498950/1‑71 (MQ=255)
          aaGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGt                                                            <  1:3770314/71‑1 (MQ=255)
               aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                                       <  1:1074279/71‑1 (MQ=255)
               aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                                       >  1:1116994/1‑71 (MQ=255)
               aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                                       <  1:703249/71‑1 (MQ=255)
               aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                                       <  1:709286/71‑1 (MQ=255)
                 gAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATg                                                         >  1:3353311/1‑67 (MQ=255)
                  aaCCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa                                                     <  1:2689415/70‑1 (MQ=255)
                  aaCCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa                                                     >  1:1468098/1‑70 (MQ=255)
                   aCCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcg                                                   >  1:3558725/1‑71 (MQ=255)
                    ccTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa                                                     >  1:473151/1‑68 (MQ=255)
                      tGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTa                                                           <  1:767103/60‑1 (MQ=255)
                       ggCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg                                                >  1:2619216/1‑70 (MQ=255)
                        tggTCGAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAAcgccgc                                               <  1:3267859/68‑1 (MQ=255)
                           tCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGt                                                                        <  1:912764/42‑1 (MQ=255)
                           tCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                           <  1:4073534/71‑1 (MQ=255)
                             gAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGa                                         >  1:1150103/1‑71 (MQ=255)
                              aaCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccgc                                               <  1:712000/64‑1 (MQ=255)
                               aCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGAt                                        >  1:3365325/1‑70 (MQ=255)
                               aCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGAt                                        >  1:3345214/1‑70 (MQ=255)
                               aCCTGTTCCATATCGGGATTCTGGGGATTTTTGTAGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATg                                       >  1:3738402/1‑71 (MQ=255)
                                  tGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTAt                                    >  1:2807950/1‑71 (MQ=255)
                                   gTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTg                                                      >  1:4077518/1‑52 (MQ=255)
                                    ttCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCgg                                                             <  1:1268234/44‑1 (MQ=255)
                                    ttCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGATGt                                      <  1:2914066/67‑1 (MQ=255)
                                        atatCGGGATTCTGGGGATTTTTGTCGGTCACttcttc                                                               <  1:446137/38‑1 (MQ=255)
                                        atatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAg                                <  1:2977024/69‑1 (MQ=255)
                                        atatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGcc                              <  1:2957944/71‑1 (MQ=255)
                                         tatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccgc                                               <  1:744199/53‑1 (MQ=255)
                                           tCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGc                               <  1:2446692/67‑1 (MQ=255)
                                           tCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTg                            <  1:402644/70‑1 (MQ=255)
                                            cGGAATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGc                          >  1:618234/1‑71 (MQ=255)
                                             gggATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAc                                                    >  1:1741810/1‑44 (MQ=255)
                                             gggATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCt                         >  1:2940257/1‑71 (MQ=255)
                                              ggATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTg                            <  1:198432/67‑1 (MQ=255)
                                                   ttGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGa                                         >  1:252109/2‑49 (MQ=16)
                                                   ttGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGa                                         >  1:1597727/2‑49 (MQ=16)
                                                   cTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa                                                     >  1:3290693/1‑37 (MQ=255)
                                                    tGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGt                                      <  1:2581900/51‑1 (MQ=255)
                                                      gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTg                 <  1:352018/70‑1 (MQ=255)
                                                      gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGa                <  1:3092479/71‑1 (MQ=255)
                                                      gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGa                <  1:1051889/71‑1 (MQ=255)
                                                          tttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAtt                  >  1:135447/1‑65 (MQ=255)
                                                             ttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa         >  1:217656/1‑71 (MQ=255)
                                                                  ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTg            >  1:3533683/1‑63 (MQ=255)
                                                                  ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGa           >  1:1408935/1‑64 (MQ=255)
                                                                  ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGa     >  1:935723/1‑70 (MQ=255)
                                                                   gTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaa   <  1:1068878/71‑1 (MQ=255)
                                                                    tCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaaa  >  1:1883278/1‑71 (MQ=255)
                                                                    tCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaaa  <  1:966850/71‑1 (MQ=255)
                                                                    |                                                                      
CTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAA  >  minE/837654‑837792

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