Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,005 T→A 23.0% S241S (TCT→TCA narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE836,0050TA23.0% 145.7 / 43.7 74S241S (TCT→TCAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (18/39);  new base A (4/13);  total (22/52)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

GTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCG  >  minE/835935‑836075
                                                                      |                                                                      
gTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCt                                                                        <  1:239832/71‑1 (MQ=255)
 tggctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTg                                                                       <  1:2128510/71‑1 (MQ=255)
  ggctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCt                                                                        >  1:2222021/1‑69 (MQ=255)
   gctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgag                                                                     <  1:203394/71‑1 (MQ=255)
     tggcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAgaaa                                                                    >  1:1152407/1‑68 (MQ=255)
       gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaaat                                                                 >  1:1719179/1‑69 (MQ=255)
       gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa                                                                   <  1:3241472/69‑1 (MQ=255)
       gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa                                                                   <  1:3246154/69‑1 (MQ=255)
       gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa                                                                   <  1:1843920/69‑1 (MQ=255)
       gcGTATGTTCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAg                                                                  >  1:2915481/1‑70 (MQ=255)
       gcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGt                                                                 >  1:3747510/1‑71 (MQ=255)
        cGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgag                                                                     <  1:934189/66‑1 (MQ=255)
          tATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa                                                              <  1:1503231/71‑1 (MQ=255)
          tATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa                                                              <  1:1710470/71‑1 (MQ=255)
           aTGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt                                                             >  1:1890133/1‑71 (MQ=255)
           aTGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt                                                             >  1:309516/1‑71 (MQ=255)
            tgtgCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtc                                                            <  1:1561817/71‑1 (MQ=255)
            tgtgCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCAtc                                                            >  1:778851/1‑71 (MQ=255)
                cATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtcttc                                                         <  1:1723363/70‑1 (MQ=255)
                  tCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGc                                                      >  1:241643/1‑71 (MQ=255)
                    aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa                                                    <  1:2942732/71‑1 (MQ=255)
                    aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa                                                    <  1:501911/71‑1 (MQ=255)
                     cTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTAt                                                   >  1:3086932/1‑71 (MQ=255)
                       ggATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATcc                                                 >  1:2096601/1‑71 (MQ=255)
                         aTGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa                                                                   <  1:2546401/51‑1 (MQ=255)
                          tGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTga                                                                      <  1:136854/47‑1 (MQ=255)
                            tccGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCcacgaa                                            <  1:2509527/70‑3 (MQ=255)
                              cGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTAtt                                          <  1:1497057/71‑1 (MQ=255)
                               gTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTg                                         <  1:960377/71‑1 (MQ=255)
                                tACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgaga                                                                    <  1:1231986/43‑1 (MQ=255)
                                  cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt                                                             >  1:2328939/1‑48 (MQ=255)
                                  cAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGag                                       <  1:1155359/70‑2 (MQ=255)
                                    aaaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtcttc                                                         <  1:3307498/50‑1 (MQ=255)
                                    aaaaaaTCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaaat                                                                 <  1:2405111/42‑3 (MQ=255)
                                     aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTa                                            <  1:3062759/62‑1 (MQ=255)
                                     aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTg                                         <  1:1252982/65‑1 (MQ=255)
                                     aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCCTCTGCTATCCGCGTAtt                                          <  1:299188/64‑1 (MQ=255)
                                      aaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTATGCTATCCGCGTATTGAAGCggg                                  <  1:3819574/71‑1 (MQ=255)
                                         aTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCa                               <  1:3498115/71‑1 (MQ=255)
                                          tCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCggg                                  <  1:3554602/67‑1 (MQ=255)
                                          tCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCggg                                  <  1:2545390/67‑1 (MQ=255)
                                          tCTACTTCAACTGGAAGAGAGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                      >  1:3360921/1‑63 (MQ=255)
                                             aCTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtctt                                                          <  1:1433275/40‑1 (MQ=255)
                                              cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGcc                            <  1:2234600/69‑1 (MQ=255)
                                              cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGcc                            <  1:3729428/69‑1 (MQ=255)
                                              cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCg                           <  1:1050029/70‑1 (MQ=255)
                                              cTTCAACTCCAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAACCGGGTCAGcc                            <  1:224040/69‑1 (MQ=255)
                                                 cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                       >  1:1174327/1‑71 (MQ=255)
                                                 cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                       >  1:2590286/1‑71 (MQ=255)
                                                 cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                       >  1:39260/1‑71 (MQ=255)
                                                  aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt                      <  1:514625/71‑1 (MQ=255)
                                                  aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt                      <  1:1229407/71‑1 (MQ=255)
                                                   aCTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTac                                                   <  1:611707/41‑2 (MQ=12)
                                                     tGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGtgtg                   >  1:2997512/1‑71 (MQ=255)
                                                       gAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCt                 <  1:534380/71‑1 (MQ=255)
                                                          gagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                                        <  1:364190/45‑1 (MQ=255)
                                                          gagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                                        <  1:359474/45‑1 (MQ=255)
                                                          gagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAg              <  1:431094/71‑1 (MQ=255)
                                                           agCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAg              >  1:2151562/1‑70 (MQ=255)
                                                           agCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGa             <  1:78065/71‑1 (MQ=255)
                                                            gCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCgg                                   <  1:1834760/48‑1 (MQ=255)
                                                            gCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGa             <  1:2602887/70‑1 (MQ=255)
                                                             cGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgtgt                    <  1:2625560/62‑1 (MQ=255)
                                                                caaGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTc                <  1:1988407/62‑1 (MQ=255)
                                                                 aaGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCgg                                   >  1:1747970/1‑43 (MQ=255)
                                                                 aaGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCgg                                   >  1:2661141/1‑43 (MQ=255)
                                                                 aaGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCt        <  1:730330/70‑1 (MQ=255)
                                                                 aaGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtg       <  1:1467286/71‑1 (MQ=255)
                                                                  aGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtgt      >  1:1183589/1‑71 (MQ=255)
                                                                  aGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtgc      <  1:1451602/71‑2 (MQ=255)
                                                                  aGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtgc      <  1:57487/71‑2 (MQ=255)
                                                                  aGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtgc      <  1:1001021/71‑2 (MQ=255)
                                                                   gTCTGAGAAGTTCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtgtg     <  1:3809722/71‑1 (MQ=255)
                                                                    tCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtgtg     <  1:3616713/70‑1 (MQ=255)
                                                                     cTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtc   >  1:1855850/1‑71 (MQ=255)
                                                                     cTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtc   >  1:2418464/1‑71 (MQ=255)
                                                                      tGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcg  >  1:713415/1‑71 (MQ=255)
                                                                      |                                                                      
GTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCG  >  minE/835935‑836075

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: