Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,716 T→C 8.7% F63F (TTT→TTC 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,7160TC 8.7% 125.6 / 4.3 46F63F (TTT→TTCnarInitrate reductase 1, gamma (cytochrome b(NR)) subunit
Reads supporting (aligned to +/- strand):  ref base T (18/24);  new base C (2/2);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAA  >  minE/837646‑837785
                                                                      |                                                                     
gccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttttt                                                                       <  1:2995246/71‑1 (MQ=255)
gccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttttt                                                                       <  1:2504515/71‑1 (MQ=255)
gccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttttt                                                                       <  1:833053/71‑1 (MQ=255)
gccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttttt                                                                       <  1:2320832/71‑1 (MQ=255)
 ccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttttt                                                                       <  1:1803875/70‑1 (MQ=255)
  cagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGt                                                                     <  1:4119496/71‑1 (MQ=255)
    gATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttttt                                                                       <  1:2157281/67‑1 (MQ=255)
        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)
                        tGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGAttttt                                                                       <  1:2504560/47‑1 (MQ=255)
                        tGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGt                                                                     <  1:1234052/49‑1 (MQ=255)
                         gAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATg                                                  >  1:3353311/1‑67 (MQ=255)
                          aaCCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa                                              >  1:1468098/1‑70 (MQ=255)
                          aaCCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa                                              <  1:2689415/70‑1 (MQ=255)
                           aCCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcg                                            >  1:3558725/1‑71 (MQ=255)
                            ccTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa                                              >  1:473151/1‑68 (MQ=255)
                              tGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTa                                                    <  1:767103/60‑1 (MQ=255)
                               ggCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGt                                                                     >  1:1537379/1‑42 (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:3345214/1‑70 (MQ=255)
                                       aCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGAt                                 >  1:3365325/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:1051889/71‑1 (MQ=255)
                                                              gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGa         <  1:3092479/71‑1 (MQ=255)
                                                                  tttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAtt           >  1:135447/1‑65 (MQ=255)
                                                                     ttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa  >  1:217656/1‑71 (MQ=255)
                                                                      |                                                                     
GCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAA  >  minE/837646‑837785

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: