Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,011 G→A 17.2% K243K (AAG→AAA 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,0110GA17.2% 150.9 / 20.9 64K243K (AAG→AAAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (18/35);  new base A (2/9);  total (20/44)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGT  >  minE/835942‑836077
                                                                     |                                                                  
gcGTATGTTCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAg                                                                    >  1:2915481/1‑70 (MQ=255)
gcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGt                                                                   >  1:3747510/1‑71 (MQ=255)
   tATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa                                                                <  1:1503231/71‑1 (MQ=255)
   tATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa                                                                <  1:1710470/71‑1 (MQ=255)
    aTGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt                                                               >  1:309516/1‑71 (MQ=255)
    aTGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt                                                               >  1:1890133/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)
                     tccGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCcacgaa                                              <  1:2509527/70‑3 (MQ=255)
                       cGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTAtt                                            <  1:1497057/71‑1 (MQ=255)
                        gTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTg                                           <  1:960377/71‑1 (MQ=255)
                           cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt                                                               >  1:2328939/1‑48 (MQ=255)
                           cAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGag                                         <  1:1155359/70‑2 (MQ=255)
                             aaaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtcttc                                                           <  1:3307498/50‑1 (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:3729428/69‑1 (MQ=255)
                                       cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGcc                              <  1:2234600/69‑1 (MQ=255)
                                       cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCg                             <  1:1050029/70‑1 (MQ=255)
                                       cTTCAACTCCAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAACCGGGTCAGcc                              <  1:224040/69‑1 (MQ=255)
                                          cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                         >  1:39260/1‑71 (MQ=255)
                                          cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                         >  1:2590286/1‑71 (MQ=255)
                                          cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                         >  1:1174327/1‑71 (MQ=255)
                                           aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt                        <  1:1229407/71‑1 (MQ=255)
                                           aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt                        <  1:514625/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)
                                                                 agaAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtctgt  <  1:3299085/71‑1 (MQ=255)
                                                                   aaGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGAcc                          <  1:3748508/45‑1 (MQ=255)
                                                                     |                                                                  
GCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGT  >  minE/835942‑836077

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: