Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,032 T→A 23.2% R250R (CGT→CGA 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,0320TA23.2% 132.0 / 40.1 69R250R (CGT→CGAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (18/35);  new base A (7/9);  total (25/44)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCT  >  minE/835965‑836100
                                                                   |                                                                    
cGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAAtt                                                                   <  1:2113350/71‑1 (MQ=255)
  tACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTAt                                                                    >  1:5108599/1‑68 (MQ=255)
  tACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTAt                                                                    >  1:1815495/1‑68 (MQ=255)
  tACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                                                                 <  1:2635482/71‑1 (MQ=255)
    cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGaa                                                                <  1:2470981/70‑1 (MQ=255)
    cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                                               >  1:975136/1‑71 (MQ=255)
        aaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCggg                                                           >  1:4372422/1‑71 (MQ=255)
         aaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGt                                                          <  1:1879659/71‑1 (MQ=255)
              tACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAAt                                                                    >  1:892203/1‑56 (MQ=255)
                cTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGagtc                                                              >  1:640696/1‑57 (MQ=255)
                cTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAAcc                                                     <  1:3325358/69‑1 (MQ=255)
                 ttCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTAtt                                                                   <  1:4046841/54‑1 (MQ=255)
                   cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAg                                                       >  1:1748580/1‑64 (MQ=255)
                   cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAg                                                       >  1:688396/1‑64 (MQ=255)
                   cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGAGGGTCAg                                                       >  1:5459720/1‑64 (MQ=255)
                    aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                                                <  1:1274678/70‑1 (MQ=255)
                       tGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCTGCGTATTGa                                                                 <  1:4650891/50‑1 (MQ=255)
                       tGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                                                                 <  1:566801/50‑1 (MQ=255)
                       tGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                                                                 <  1:3360763/50‑1 (MQ=255)
                       tGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                                                                 <  1:1284842/50‑1 (MQ=255)
                       tGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                                                                 <  1:42145/50‑1 (MQ=255)
                       tGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAAtt                                                                   <  1:4161880/48‑1 (MQ=255)
                        ggAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGc                                           >  1:1181434/1‑71 (MQ=255)
                         gAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                                               >  1:3881858/1‑50 (MQ=255)
                         gAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                                               >  1:4733433/1‑50 (MQ=255)
                         gAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                                               >  1:1831097/1‑50 (MQ=255)
                         gAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGc                                           <  1:2373275/70‑1 (MQ=255)
                              gCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGaa                                     <  1:3124529/71‑1 (MQ=255)
                              gCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGaa                                     <  1:2603101/71‑1 (MQ=255)
                                 gTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                                                <  1:3512480/57‑1 (MQ=255)
                                  tAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                                                                 <  1:4125431/39‑1 (MQ=255)
                                  tAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGTGTCAGCCGACCGTGTGCTCAGAAACCt                                 <  1:3052700/71‑1 (MQ=255)
                                    aGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCg                                                             >  1:807190/1‑41 (MQ=255)
                                    aGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtgt                               >  1:1209816/1‑71 (MQ=255)
                                    aGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtgc                               >  1:1711302/1‑70 (MQ=255)
                                      tCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGc                                           <  1:3714398/57‑1 (MQ=255)
                                      tCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGc                                           <  1:3871235/57‑1 (MQ=255)
                                      tCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGc                                           <  1:4415156/57‑1 (MQ=255)
                                        tGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCgg                                                            >  1:4145982/1‑38 (MQ=255)
                                        tGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCgg                                                            >  1:3008409/1‑38 (MQ=255)
                                          agaAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcggt                         <  1:4389993/71‑1 (MQ=255)
                                          agaAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcgg                          <  1:508103/70‑1 (MQ=255)
                                            aaatGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCgtgggtc                        >  1:777092/4‑70 (MQ=255)
                                             aatGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGaa                                     <  1:1557606/54‑1 (MQ=255)
                                             aatGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGaa                                     <  1:4917612/54‑1 (MQ=255)
                                             aatGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCgtgggtcgc                      <  1:1656730/69‑2 (MQ=255)
                                             aGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGGCGGTCGt                      <  1:523969/71‑1 (MQ=255)
                                              atGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCgtgggtcgca                     >  1:3327909/2‑69 (MQ=255)
                                               tGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTAt                    <  1:401398/71‑1 (MQ=255)
                                                gCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAAcc                                  <  1:1254123/56‑1 (MQ=255)
                                                  aTCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATcct                 <  1:5139350/71‑2 (MQ=255)
                                                  aTCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCg                 <  1:1523005/71‑1 (MQ=255)
                                                  aTCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCg                 <  1:4443711/71‑1 (MQ=255)
                                                  aTCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAAACTGTGTCGGTCGTATCCg                 <  1:566046/71‑1 (MQ=255)
                                                   tcttctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGtgtg                                            <  1:2499695/43‑1 (MQ=255)
                                                    cttctGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGa                                      <  1:4473373/48‑1 (MQ=255)
                                                     ttctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTa              >  1:5082527/1‑71 (MQ=255)
                                                      tctGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTAt             >  1:3191624/1‑71 (MQ=255)
                                                      tctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTAt             >  1:58117/1‑71 (MQ=255)
                                                       ctGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCgtgggtcgca                     <  1:2475869/62‑3 (MQ=255)
                                                       ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGtt               <  1:5263314/68‑1 (MQ=255)
                                                       ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGtt               <  1:675981/68‑1 (MQ=255)
                                                       ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTa              <  1:5161340/69‑1 (MQ=255)
                                                          cTATCCGCGTATTTAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCt           <  1:4448824/69‑1 (MQ=255)
                                                             tCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCtt          <  1:2269130/67‑1 (MQ=255)
                                                             tCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCtt          <  1:2319512/67‑1 (MQ=255)
                                                               cgcgTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGt     <  1:1879847/70‑1 (MQ=255)
                                                                gcgTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGc   >  1:2062906/1‑71 (MQ=255)
                                                                 cgTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCt  >  1:208666/1‑71 (MQ=255)
                                                                   |                                                                    
CGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCT  >  minE/835965‑836100

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: