Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,052 T→G 9.8% V1170V (GTT→GTG) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,0520TG 9.8% 241.9 / 19.2 101V1170V (GTT→GTG) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (39/52);  new base G (2/8);  total (41/60)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

CCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACGCCG  >  minE/834982‑835121
                                                                      |                                                                     
cccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                                       <  1:4116101/71‑1 (MQ=255)
    tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                                       <  1:679706/67‑1 (MQ=255)
    tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                                       <  1:3040612/67‑1 (MQ=255)
    tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                                       <  1:2860641/67‑1 (MQ=255)
    tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcc                                                                   >  1:2168996/1‑71 (MQ=255)
     gCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcct                                                                  >  1:1216563/1‑71 (MQ=255)
      cGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctg                                                                 >  1:3045363/1‑71 (MQ=255)
         ttGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctgcct                                                              <  1:1044911/71‑1 (MQ=255)
         ttGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctgcct                                                              <  1:3365965/71‑1 (MQ=255)
         ttGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctgcct                                                              <  1:2602551/71‑1 (MQ=255)
           gTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTgg                                                            >  1:4585703/1‑71 (MQ=255)
            tCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                                       <  1:745313/59‑1 (MQ=255)
             cagccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGtt                                                          >  1:4710332/1‑71 (MQ=255)
              agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc                                                         >  1:832669/1‑71 (MQ=255)
              agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc                                                         >  1:168575/1‑71 (MQ=255)
              agccaacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTc                                                         >  1:3183756/1‑71 (MQ=255)
               gccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                                                                      <  1:2997136/57‑1 (MQ=255)
               gccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                                                                      <  1:5384197/57‑1 (MQ=255)
               gccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                                                                      <  1:4065311/57‑1 (MQ=255)
                ccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCgg                                                       <  1:4492449/71‑1 (MQ=255)
                ccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCgg                                                       <  1:641250/71‑1 (MQ=255)
                ccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGCTTCgg                                                       <  1:2188890/71‑1 (MQ=255)
                ccaacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCgg                                                       <  1:3500068/71‑1 (MQ=255)
                  agcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa                                                      >  1:1111952/1‑70 (MQ=255)
                  aacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCg                                                        <  1:1771948/66‑1 (MQ=255)
                  aacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaa                                                     <  1:4544682/69‑1 (MQ=255)
                    cGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCgg                                                       <  1:5510869/67‑1 (MQ=255)
                    cGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCgg                                                       <  1:2805172/67‑1 (MQ=255)
                     gtgtACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaagta                                                  <  1:4552860/71‑4 (MQ=255)
                     gtgtACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaagt                                                   >  1:540723/1‑68 (MQ=255)
                        ttCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa                                                      <  1:3109862/64‑1 (MQ=255)
                        ttCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAcc                                               >  1:3714975/1‑71 (MQ=255)
                          ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGaaa                                                    >  1:5110952/1‑64 (MQ=255)
                          ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc                                              <  1:266157/70‑1 (MQ=255)
                          ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc                                              >  1:1956105/1‑70 (MQ=255)
                          ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc                                              >  1:1409111/1‑70 (MQ=255)
                          ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc                                              >  1:5149989/1‑70 (MQ=255)
                           cGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCaa                                            <  1:100192/71‑1 (MQ=255)
                           cGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCa                                             <  1:2256105/70‑1 (MQ=255)
                              cAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAg                                         <  1:5439251/71‑1 (MQ=255)
                               aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCaa                                            >  1:3655850/1‑67 (MQ=255)
                               aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                                        >  1:5220430/1‑71 (MQ=255)
                               aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                                        <  1:2919289/71‑1 (MQ=255)
                               aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                                        <  1:711814/71‑1 (MQ=255)
                               aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                                        >  1:2105595/1‑71 (MQ=255)
                               aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                                        <  1:4290125/71‑1 (MQ=255)
                                gggATGACCATGATGTACCACGCGCAGTAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCCACAGCg                                       <  1:740230/71‑1 (MQ=255)
                                gggATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                                                                      <  1:4034070/40‑1 (MQ=255)
                                gggATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAGATTACCCAACAGCg                                       <  1:1989126/71‑1 (MQ=255)
                                gggATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg                                       >  1:32498/1‑71 (MQ=255)
                                gggATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg                                       <  1:1873344/71‑1 (MQ=255)
                                   aTGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCgtgg                                    <  1:4157222/71‑1 (MQ=255)
                                   aTGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCgtgg                                    <  1:1870127/71‑1 (MQ=255)
                                    tGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaagt                                                   <  1:5560455/55‑3 (MQ=255)
                                    tGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaagt                                                   <  1:5571842/55‑3 (MQ=255)
                                    tGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaagt                                                   <  1:2520382/55‑3 (MQ=255)
                                     gACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCgg                                                       <  1:3391191/50‑1 (MQ=255)
                                     gACCATGATGTACCACGCGCAGGAACGTATCGTTAAcc                                                                   <  1:1715527/38‑1 (MQ=255)
                                     gACCATGATGTACCACGCGCAGGAACGTATCGTTAAcc                                                                   <  1:4516214/38‑1 (MQ=255)
                                     gACCATGATGTACCACGCGCAGGAACGTATCGTTAAcc                                                                   <  1:4424211/38‑1 (MQ=255)
                                      aCCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAc                                                >  1:3980459/1‑56 (MQ=255)
                                       ccATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTggtggt                                <  1:3594538/71‑1 (MQ=255)
                                          tgatgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa                                                      <  1:5377771/46‑1 (MQ=255)
                                          tgatgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGtggtgg                                 >  1:3078849/1‑67 (MQ=255)
                                          tgatgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATc                             >  1:4629983/1‑71 (MQ=255)
                                           gatgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTggtggt                                <  1:5001685/67‑1 (MQ=255)
                                            atgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATcc                            <  1:5501726/70‑1 (MQ=255)
                                            atgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATcc                            <  1:4719251/70‑1 (MQ=255)
                                            atgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCa                           <  1:2565994/71‑1 (MQ=255)
                                             tgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa                                                      <  1:1175358/43‑1 (MQ=255)
                                             tgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCAt                          <  1:2338144/71‑1 (MQ=255)
                                             tgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCAt                          <  1:755176/71‑1 (MQ=255)
                                             tgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCAt                          <  1:1066530/71‑1 (MQ=255)
                                              gTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCAt                          >  1:1364143/1‑70 (MQ=255)
                                                 ccACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACt                      <  1:1293248/71‑1 (MQ=255)
                                                  cACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCa                           <  1:4578529/65‑1 (MQ=255)
                                                   aCGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCg                    >  1:316830/1‑71 (MQ=255)
                                                    cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCgt                                      >  1:656371/1‑52 (MQ=255)
                                                    cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg                                       <  1:3891796/51‑1 (MQ=255)
                                                    cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg                                       <  1:347520/51‑1 (MQ=255)
                                                    cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg                                       <  1:1304411/51‑1 (MQ=255)
                                                    cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCgg                   >  1:4904930/1‑71 (MQ=255)
                                                    cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCgg                   >  1:3817086/1‑71 (MQ=255)
                                                    cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCgg                   >  1:3692524/1‑71 (MQ=255)
                                                       gcAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCa                >  1:179856/1‑71 (MQ=255)
                                                          ggAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCAccc             >  1:1694899/1‑71 (MQ=255)
                                                          ggAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTACCCg            >  1:307274/1‑71 (MQ=255)
                                                           gAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTAt                              <  1:793558/53‑1 (MQ=255)
                                                           gAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCg            >  1:4177929/1‑71 (MQ=255)
                                                           gAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCg            >  1:2103356/1‑71 (MQ=255)
                                                               gTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTAt         >  1:3971130/1‑70 (MQ=255)
                                                                tATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCgtgg                                    <  1:3339139/42‑1 (MQ=255)
                                                                tATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCgtgg                                    <  1:2803009/42‑1 (MQ=255)
                                                                  tCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACg     <  1:2958272/71‑1 (MQ=255)
                                                                   cGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACg     <  1:3288410/70‑1 (MQ=255)
                                                                     ttAACCTGCCTGGTTCGGAAATTATCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTAt         <  1:5608505/64‑1 (MQ=255)
                                                                     ttAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACGCCg  <  1:4973534/71‑1 (MQ=255)
                                                                      tAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTa                               >  1:4839872/1‑41 (MQ=255)
                                                                      tAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACGCCg  >  1:2134356/1‑70 (MQ=255)
                                                                      tAACCTGCCTGGTTCAGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACGCCg  >  1:3192553/1‑70 (MQ=255)
                                                                      |                                                                     
CCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACGCCG  >  minE/834982‑835121

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: