Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,740 G→A 8.1% L71L (CTG→CTA) ‡ 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,7400GA 8.1% 153.7 / 3.6 49L71L (CTG→CTA) ‡narInitrate reductase 1, gamma (cytochrome b(NR)) subunit
Reads supporting (aligned to +/- strand):  ref base G (24/21);  new base A (2/2);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

AACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGT  >  minE/837672‑837801
                                                                    |                                                             
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)
     ggCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg                                                         >  1:2619216/1‑70 (MQ=255)
      tggTCGAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAAcgccgc                                                        <  1:3267859/68‑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)
                  ttCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGATGt                                               <  1:2914066/67‑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)
                                                    aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGa                             >  1:3648140/1‑51 (MQ=255)
                                                    aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAAt          <  1:3360085/70‑1 (MQ=255)
                                                    aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAAt          <  1:2728200/70‑1 (MQ=255)
                                                    aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATg         <  1:133515/71‑1 (MQ=255)
                                                        cttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGaa                        <  1:752778/52‑1 (MQ=255)
                                                        cttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaaa           >  1:3100994/1‑65 (MQ=255)
                                                        cttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCaa     >  1:291927/1‑71 (MQ=255)
                                                          tcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATg   >  1:15578/1‑71 (MQ=255)
                                                           cGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGt  >  1:738442/1‑71 (MQ=255)
                                                              tATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAtt                           <  1:3389569/43‑1 (MQ=255)
                                                                tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaa             <  1:2909760/55‑1 (MQ=255)
                                                                    |                                                             
AACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGT  >  minE/837672‑837801

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