Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,738 C→T 7.4% L71L (CTG→TTG) ‡ 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,7380CT 7.4% 173.2 / 3.4 54L71L (CTG→TTG) ‡narInitrate reductase 1, gamma (cytochrome b(NR)) subunit
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (2/2);  total (28/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 = 7.87e-01

ATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGT  >  minE/837669‑837801
                                                                     |                                                               
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)
   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)
        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: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)
                                                   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)
                                                                     |                                                               
ATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGT  >  minE/837669‑837801

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: