Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,037 G→C 19.1% A1165A (GCG→GCC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,0370GC19.1% 262.2 / 49.5 115A1165A (GCG→GCCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (35/58);  new base C (7/15);  total (42/73)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

GGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCA  >  minE/834969‑835107
                                                                    |                                                                      
ggTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCAcgcg                                                                        >  1:4350576/1‑69 (MQ=255)
 gTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAg                                                                     <  1:1291937/71‑1 (MQ=255)
  tGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAgg                                                                    <  1:859508/71‑1 (MQ=255)
    cgctGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGaa                                                                  <  1:412202/69‑1 (MQ=255)
      ctGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACg                                                                >  1:447976/1‑71 (MQ=255)
       tGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTTCGGCAGGGATGACCATGATGTACCACGCGCAGGAACg                                                                >  1:4915988/1‑70 (MQ=255)
        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)
             cccGCGCGGTGGTCAACCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCAttat                                                          >  1:45183/1‑67 (MQ=255)
              ccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGt                                                               <  1:1621651/64‑1 (MQ=255)
              ccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGt                                                               <  1:3891855/64‑1 (MQ=255)
              ccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGt                                                               <  1:3707454/64‑1 (MQ=255)
                gTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTAt                                                             <  1:2908133/64‑1 (MQ=255)
                 tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                         <  1:3040612/67‑1 (MQ=255)
                 tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                         <  1:2860641/67‑1 (MQ=255)
                 tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                         <  1:679706/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)
                        gTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACg                                                                >  1:1033878/1‑53 (MQ=255)
                        gTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTgg                                              >  1:4585703/1‑71 (MQ=255)
                        gTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCa                                                                      <  1:2518405/47‑1 (MQ=16)
                         tCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCAt                                                             <  1:3224699/55‑1 (MQ=255)
                          cagccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGtt                                            >  1:4710332/1‑71 (MQ=255)
                          cagccaacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAgg                                                                    >  1:2456129/1‑47 (MQ=16)
                          cagccaacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACgc                                                               >  1:5210683/1‑51 (MQ=16)
                           agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc                                           >  1:2082186/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)
                           agccaacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTc                                           >  1:3183756/1‑71 (MQ=255)
                            gccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                                                        <  1:4065311/57‑1 (MQ=255)
                            gccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                                                        <  1:5384197/57‑1 (MQ=255)
                            gccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                                                        <  1:2997136/57‑1 (MQ=255)
                             ccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCgg                                         <  1:641250/71‑1 (MQ=255)
                             ccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCgg                                         <  1:4492449/71‑1 (MQ=255)
                             ccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGCTTCgg                                         <  1:2188890/71‑1 (MQ=255)
                               agcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa                                        <  1:5320507/70‑1 (MQ=255)
                               agcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa                                        >  1:1111952/1‑70 (MQ=255)
                                acGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGaa                                                                  >  1:4853511/2‑43 (MQ=12)
                                 cGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCgg                                         <  1:3043868/67‑1 (MQ=255)
                                  gtgtACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaagta                                    <  1:4552860/71‑4 (MQ=255)
                                  gtgtACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaagt                                     >  1:540723/1‑68 (MQ=255)
                                   tgtTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGt                                                          <  1:5269924/48‑1 (MQ=255)
                                   tgtTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGt                                                          <  1:4029514/48‑1 (MQ=255)
                                       ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGaaa                                      >  1:630580/1‑64 (MQ=255)
                                       ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGaaa                                      >  1:5110952/1‑64 (MQ=255)
                                       ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc                                >  1:1956105/1‑70 (MQ=255)
                                       ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc                                >  1:1409111/1‑70 (MQ=255)
                                       ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc                                >  1:5149989/1‑70 (MQ=255)
                                       ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc                                <  1:266157/70‑1 (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:4290125/71‑1 (MQ=255)
                                            aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                          <  1:2919289/71‑1 (MQ=255)
                                            aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                          >  1:2105595/1‑71 (MQ=255)
                                            aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                          <  1:3324322/71‑1 (MQ=255)
                                            aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                          >  1:5220430/1‑71 (MQ=255)
                                            aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc                          <  1:711814/71‑1 (MQ=255)
                                             gggATGACCATGATGTACCACGCGCAGTAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCCACAGCg                         <  1:740230/71‑1 (MQ=255)
                                             gggATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg                         <  1:1873344/71‑1 (MQ=255)
                                             gggATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg                         >  1:32498/1‑71 (MQ=255)
                                                aTGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCgtgg                      <  1:4157222/71‑1 (MQ=255)
                                                aTGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCgtgg                      <  1:1870127/71‑1 (MQ=255)
                                                 tGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaagt                                     <  1:5571842/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)
                                                  gACCATGATGTACCACGCGCAGGAACGTATCGTTAAcc                                                     <  1:1715527/38‑1 (MQ=255)
                                                   aCCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAc                                  >  1:3980459/1‑56 (MQ=255)
                                                    ccATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTggtggt                  <  1:3594538/71‑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:4719251/70‑1 (MQ=255)
                                                         atgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATcc              <  1:5501726/70‑1 (MQ=255)
                                                         atgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCa             <  1:2565994/71‑1 (MQ=255)
                                                          tgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCAt            <  1:1066530/71‑1 (MQ=255)
                                                          tgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCAt            <  1:755176/71‑1 (MQ=255)
                                                          tgTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCAt            <  1:2338144/71‑1 (MQ=255)
                                                            tACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATaa          >  1:4669266/1‑71 (MQ=255)
                                                              ccACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACt        <  1:1293248/71‑1 (MQ=255)
                                                               cACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGaaa                                      >  1:2583663/1‑40 (MQ=255)
                                                               cACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCa             <  1:4578529/65‑1 (MQ=255)
                                                                aCGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCg      >  1:316830/1‑71 (MQ=255)
                                                                 cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg                         <  1:1304411/51‑1 (MQ=255)
                                                                 cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg                         <  1:3891796/51‑1 (MQ=255)
                                                                 cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg                         <  1:347520/51‑1 (MQ=255)
                                                                 cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCgg     >  1:3817086/1‑71 (MQ=255)
                                                                 cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCgg     >  1:3692524/1‑71 (MQ=255)
                                                                 cgcgCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCgg     >  1:4904930/1‑71 (MQ=255)
                                                                    gcAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCa  >  1:179856/1‑71 (MQ=255)
                                                                    |                                                                      
GGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCA  >  minE/834969‑835107

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: