Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,046 T→C 12.4% R1168R (CGT→CGC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,0460TC12.4% 284.1 / 25.2 113R1168R (CGT→CGCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (42/57);  new base C (4/10);  total (46/67)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

GACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTAT  >  minE/834977‑835114
                                                                     |                                                                    
gACTGCCCGTGCGGTTGTCAGGCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTa                                                                     >  1:4388736/1‑71 (MQ=255)
gACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTa                                                                     >  1:2540680/1‑71 (MQ=255)
     cccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCg                                                                  <  1:661469/69‑1 (MQ=255)
     cccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                                <  1:4116101/71‑1 (MQ=255)
     cccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGt                                                                 >  1:2246713/1‑70 (MQ=255)
     cccGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCAttat                                                                 <  1:3019946/70‑4 (MQ=255)
     cccGCGCGGTGGTCAACCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCAttat                                                                 >  1:45183/1‑67 (MQ=255)
      ccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGt                                                                      <  1:3707454/64‑1 (MQ=255)
      ccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGt                                                                      <  1:1621651/64‑1 (MQ=255)
        gTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTAt                                                                    <  1:2908133/64‑1 (MQ=255)
         tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                                <  1:2860641/67‑1 (MQ=255)
         tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                                <  1:679706/67‑1 (MQ=255)
         tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                                <  1:3040612/67‑1 (MQ=255)
         tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcc                                                            >  1:2168996/1‑71 (MQ=255)
          gCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcct                                                           >  1:1216563/1‑71 (MQ=255)
           cGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctg                                                          >  1:4223642/1‑71 (MQ=255)
           cGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctg                                                          >  1:3045363/1‑71 (MQ=255)
              ttGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctgcct                                                       <  1:3365965/71‑1 (MQ=255)
              ttGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctgcct                                                       <  1:1044911/71‑1 (MQ=255)
                gTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTgg                                                     >  1:4585703/1‑71 (MQ=255)
                 tCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                                <  1:745313/59‑1 (MQ=255)
                 tCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCAt                                                                    <  1:3224699/55‑1 (MQ=255)
                  cagccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGtt                                                   >  1:4710332/1‑71 (MQ=255)
                   agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc                                                  >  1:839341/1‑71 (MQ=255)
                   agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc                                                  >  1:832669/1‑71 (MQ=255)
                   agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc                                                  >  1:168575/1‑71 (MQ=255)
                   agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc                                                  >  1:2082186/1‑71 (MQ=255)
                   agccaacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTc                                                  >  1:3183756/1‑71 (MQ=255)
                    gccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                                                               <  1:4065311/57‑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)
                       agcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa                                               <  1:5320507/70‑1 (MQ=255)
                       aacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCg                                                 <  1:1771948/66‑1 (MQ=255)
                       aacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaa                                              <  1:4544682/69‑1 (MQ=255)
                         cGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCgg                                                <  1:3043868/67‑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)
                           tgtTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGt                                                                 <  1:4901642/48‑1 (MQ=255)
                           tgtTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGt                                                                 <  1:5269924/48‑1 (MQ=255)
                           tgtTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGt                                                                 <  1:4029514/48‑1 (MQ=255)
                             ttCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa                                               <  1:3109862/64‑1 (MQ=255)
                             ttCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAcc                                        >  1:3714975/1‑71 (MQ=255)
                               ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGaaa                                             >  1:5110952/1‑64 (MQ=255)
                               ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGaaa                                             >  1:630580/1‑64 (MQ=255)
                               ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc                                       >  1:5149989/1‑70 (MQ=255)
                               ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc                                       >  1:1956105/1‑70 (MQ=255)
                               ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc                                       <  1:266157/70‑1 (MQ=255)
                               ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc                                       >  1:1409111/1‑70 (MQ=255)
                                cGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCaa                                     <  1:100192/71‑1 (MQ=255)
                                   cAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAg                                  <  1:5439251/71‑1 (MQ=255)
                                    aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                                 <  1:4290125/71‑1 (MQ=255)
                                    aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                                 >  1:5220430/1‑71 (MQ=255)
                                    aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                                 <  1:3324322/71‑1 (MQ=255)
                                    aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                                 >  1:2105595/1‑71 (MQ=255)
                                    aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                                 <  1:2919289/71‑1 (MQ=255)
                                    aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                                 <  1:711814/71‑1 (MQ=255)
                                     gggATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                                                               <  1:4034070/40‑1 (MQ=255)
                                     gggATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg                                <  1:1873344/71‑1 (MQ=255)
                                     gggATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg                                >  1:32498/1‑71 (MQ=255)
                                        aTGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCgtgg                             <  1:1870127/71‑1 (MQ=255)
                                         tGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaagt                                            <  1:2520382/55‑3 (MQ=255)
                                         tGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaagt                                            <  1:5560455/55‑3 (MQ=255)
                                          gACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCgg                                                <  1:3391191/50‑1 (MQ=255)
                                          gACCATGATGTACCACGCGCAGGAACGTATCGTTAAcc                                                            <  1:4424211/38‑1 (MQ=255)
                                          gACCATGATGTACCACGCGCAGGAACGTATCGTTAAcc                                                            <  1:4516214/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)
                                                 atgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATcc                     <  1:4719251/70‑1 (MQ=255)
                                                 atgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATcc                     <  1:5501726/70‑1 (MQ=255)
                                                 atgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCa                    <  1:2565994/71‑1 (MQ=255)
                                                  tgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa                                               <  1:1175358/43‑1 (MQ=255)
                                                  tgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCAt                   <  1:1066530/71‑1 (MQ=255)
                                                  tgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCAt                   <  1:755176/71‑1 (MQ=255)
                                                  tgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCAt                   <  1:2338144/71‑1 (MQ=255)
                                                   gTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCAt                   >  1:1364143/1‑70 (MQ=255)
                                                    tACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATaa                 >  1:4669266/1‑71 (MQ=255)
                                                       cACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGaaa                                             >  1:2583663/1‑40 (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:1304411/51‑1 (MQ=255)
                                                         cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg                                <  1:347520/51‑1 (MQ=255)
                                                         cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCgg            >  1:3692524/1‑71 (MQ=255)
                                                         cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCgg            >  1:3817086/1‑71 (MQ=255)
                                                         cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCgg            >  1:4904930/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:2103356/1‑71 (MQ=255)
                                                                gAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCg     >  1:4177929/1‑71 (MQ=255)
                                                                    gTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTAt  >  1:3750120/1‑70 (MQ=255)
                                                                    gTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTAt  >  1:3971130/1‑70 (MQ=255)
                                                                     tATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCgtgg                             <  1:2803009/42‑1 (MQ=255)
                                                                     |                                                                    
GACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTAT  >  minE/834977‑835114

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: