Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,146 T→G 16.0% S202A (TCG→GCG) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,1460TG16.0% 196.2 / 40.4 100S202A (TCG→GCG) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (54/30);  new base G (6/10);  total (60/40)
Fisher's exact test for biased strand distribution p-value = 5.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

Only 100 of 110 total aligned reads displayed.

AAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCC  >  minE/832077‑832216
                                                                     |                                                                      
aaaaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCAt                                                                        <  1:2540267/70‑1 (MQ=255)
           cccGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACg                                                              >  1:3754798/1‑69 (MQ=255)
            ccGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCt                                                            >  1:1641594/1‑70 (MQ=255)
                 ccGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtctct                                                      >  1:3760901/1‑71 (MQ=255)
                    tgtTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCg                                                    >  1:4606854/1‑70 (MQ=255)
                    tgtTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCg                                                    >  1:5561017/1‑70 (MQ=255)
                      ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCAt                                                                        <  1:4333045/48‑1 (MQ=255)
                      ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCAt                                                                        <  1:3673273/48‑1 (MQ=255)
                      ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                  >  1:2834931/1‑70 (MQ=255)
                      ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTg                                                 <  1:4234690/71‑1 (MQ=255)
                        gCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGc                                                   >  1:4435639/1‑67 (MQ=255)
                          cggTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTg                                             <  1:4056302/70‑1 (MQ=255)
                          cggTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTg                                             <  1:5681275/70‑1 (MQ=255)
                           ggTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTgg                                            >  1:3516300/1‑71 (MQ=255)
                              ttCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCgg                                         >  1:309090/1‑71 (MQ=255)
                              ttCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCgg                                         >  1:797299/1‑71 (MQ=255)
                               tctcGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCgg                                         >  1:5516512/1‑70 (MQ=255)
                                   gCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTgct                                               >  1:907774/1‑58 (MQ=255)
                                   gCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTgct                                               >  1:5337072/1‑58 (MQ=255)
                                   gCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTgct                                               >  1:3490314/1‑58 (MQ=255)
                                   gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGc                                                             <  1:4074580/46‑1 (MQ=255)
                                   gCCAATTCCGGCAACGTCGATGGTTTCTTACGCATCGGGTGCACGc                                                             <  1:5222526/46‑1 (MQ=255)
                                     cAATTCCGGCAATGTCGATGGTTTCTTACGCATCggg                                                                    >  1:4835259/1‑37 (MQ=255)
                                     cAATTCCGGCAATGTCGATGGTTTCTTACGCATCg                                                                      >  1:832534/1‑35 (MQ=255)
                                     cAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGc                                                             >  1:2409515/1‑44 (MQ=255)
                                      aaTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTg                                             <  1:5503313/59‑1 (MQ=255)
                                      aaTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTg                                             <  1:110474/59‑1 (MQ=255)
                                      aaTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTg                                             <  1:5495926/59‑1 (MQ=255)
                                       atcccGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGttt                                <  1:921405/68‑1 (MQ=255)
                                       atcccGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGttt                                <  1:3742798/68‑1 (MQ=255)
                                       aTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACt                                     <  1:927527/66‑1 (MQ=255)
                                       aTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCtt                                <  1:4752251/71‑1 (MQ=255)
                                       aTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCtt                                <  1:3642011/71‑1 (MQ=255)
                                       aTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCt                                 >  1:4678871/1‑70 (MQ=255)
                                       aTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCt                                 >  1:4341106/1‑70 (MQ=255)
                                        ttCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAt                                               >  1:898989/1‑55 (MQ=255)
                                        ttCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAt                                               >  1:3654107/1‑55 (MQ=255)
                                          ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                >  1:4061991/1‑52 (MQ=255)
                                          ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAg                             >  1:2197797/1‑71 (MQ=255)
                                          ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAg                             >  1:592971/1‑71 (MQ=255)
                                           cGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTa                               <  1:1411874/68‑1 (MQ=255)
                                           cGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTa                               <  1:1759196/68‑1 (MQ=255)
                                            ggCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCt                           <  1:5232913/71‑1 (MQ=255)
                                             gCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCtt                          >  1:1360984/1‑71 (MQ=255)
                                             gCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCtt                          >  1:1552504/1‑71 (MQ=255)
                                             gCAATGTCGATGGTTTCTTACGCATCGGGTGAACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCt                           >  1:3361782/1‑70 (MQ=255)
                                              cAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGTCGGTACTTGCTTAAGCTTc                         <  1:5181192/71‑1 (MQ=255)
                                              cAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTc                         <  1:1379888/71‑1 (MQ=255)
                                              cAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTc                         <  1:2256318/71‑1 (MQ=255)
                                                aTGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGt                                        <  1:3089293/54‑1 (MQ=255)
                                                aTGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACt                                     >  1:2705791/1‑57 (MQ=255)
                                                aTGTCGATGGTTTCTTACGCATCGGGTGAACGCTATCTCTCGCTGATTGGCGGTACt                                     >  1:1341708/1‑57 (MQ=255)
                                                 tGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                >  1:974481/1‑45 (MQ=255)
                                                 tGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTAc                      >  1:1451138/1‑71 (MQ=255)
                                                 tGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTAc                      <  1:4317805/71‑1 (MQ=255)
                                                   tCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTa                               <  1:1159936/60‑1 (MQ=255)
                                                   tCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTa                               <  1:145463/60‑1 (MQ=255)
                                                   tCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                  >  1:3910083/1‑41 (MQ=255)
                                                   tCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                >  1:2618927/1‑43 (MQ=255)
                                                   tCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                >  1:5602990/1‑43 (MQ=255)
                                                   tCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGa                    >  1:1965010/1‑71 (MQ=255)
                                                   tCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGa                    >  1:4296815/1‑71 (MQ=255)
                                                    cGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                >  1:4349373/1‑42 (MQ=255)
                                                    cGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                >  1:2749841/1‑42 (MQ=255)
                                                     gATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACt                  >  1:2149607/1‑71 (MQ=255)
                                                     gATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACt                  >  1:2666366/1‑71 (MQ=255)
                                                     gATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACt                  >  1:2233995/1‑71 (MQ=255)
                                                       tGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                >  1:1453050/1‑39 (MQ=255)
                                                       tGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTgg                <  1:2098656/71‑1 (MQ=255)
                                                        ggTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTgg                >  1:4249428/1‑70 (MQ=255)
                                                         gTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGt               <  1:5306882/70‑1 (MQ=255)
                                                         gTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTa              >  1:1920521/1‑71 (MQ=255)
                                                         gTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTa              >  1:1816686/1‑71 (MQ=255)
                                                          tttCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTAc             >  1:586300/1‑71 (MQ=255)
                                                           ttCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACt                  <  1:1504544/65‑1 (MQ=255)
                                                           ttCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACt                  <  1:966277/65‑1 (MQ=255)
                                                           ttCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCtt                                <  1:636632/51‑1 (MQ=255)
                                                           ttCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCtt                                <  1:828072/51‑1 (MQ=255)
                                                            tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTAc                                      >  1:2949811/1‑44 (MQ=255)
                                                            tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACt                  >  1:5542322/1‑64 (MQ=255)
                                                            tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACt            <  1:428408/70‑1 (MQ=255)
                                                            tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACt            <  1:4581840/70‑1 (MQ=255)
                                                            tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTg           >  1:3482011/1‑71 (MQ=255)
                                                              ttACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTc                         >  1:291284/1‑55 (MQ=255)
                                                                aCGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCtt                                >  1:4341842/1‑46 (MQ=255)
                                                                aCGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCtt                                >  1:4700812/1‑46 (MQ=255)
                                                                aCGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGc          <  1:4645435/68‑1 (MQ=255)
                                                                 cGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTAc                      <  1:5182964/55‑1 (MQ=255)
                                                                 cGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACt      <  1:4281585/71‑1 (MQ=255)
                                                                  gcttCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGc          <  1:798671/63‑1 (MQ=255)
                                                                  gCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt     <  1:3784066/71‑1 (MQ=255)
                                                                    aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGAc       >  1:3373802/1‑67 (MQ=255)
                                                                    aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGAc       >  1:2123407/1‑67 (MQ=255)
                                                                    aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGc   >  1:3964838/1‑71 (MQ=255)
                                                                     tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCtt                                <  1:2032844/41‑1 (MQ=255)
                                                                     tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCtt                                <  1:4734857/41‑1 (MQ=255)
                                                                     tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCtt                                <  1:38016/41‑1 (MQ=255)
                                                                     tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc  <  1:4910086/71‑1 (MQ=255)
                                                                     tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc  <  1:2835482/71‑1 (MQ=255)
                                                                     tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc  <  1:121244/71‑1 (MQ=255)
                                                                     |                                                                      
AAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCC  >  minE/832077‑832216

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: