Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,966 C→T 48.7% H808H (CAC→CAT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,9660CT48.7% 5.3 / 105.2 80H808H (CAC→CATnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (12/29);  new base T (16/23);  total (28/52)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01

CATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCA  >  minE/833897‑834030
                                                                     |                                                                
catTTTGCCGACGGCGACCTGGTACGAAAAAGAAGACATGAATACTTCGGATATGCATCCGTTTATTCAcc                                                                 >  1:2473935/1‑71 (MQ=255)
caatttGCCGACGGCGACCTGGTACGAAAAAGACTACATGAATACTTCGGATATGCATCCGTTTATTCAcc                                                                 >  1:1915154/4‑71 (MQ=255)
 atTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAccc                                                                >  1:1279436/1‑71 (MQ=255)
      gCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTg                                                            <  1:3934286/70‑1 (MQ=255)
           cGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcg                                                      >  1:182574/1‑71 (MQ=255)
             gCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcggc                                                    <  1:4778080/71‑1 (MQ=255)
             gCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcggc                                                    <  1:2071345/71‑1 (MQ=255)
              cGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTg                                                            <  1:1610348/62‑1 (MQ=255)
                  cTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTc                                                >  1:3764463/1‑70 (MQ=255)
                  cTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCg                                               <  1:495935/71‑1 (MQ=255)
                     gTACGAAAAAGACTATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGa                                              <  1:2753154/69‑1 (MQ=255)
                     gTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCttt                                                           >  1:2928721/1‑54 (MQ=255)
                     gTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGa                                              <  1:774894/69‑1 (MQ=255)
                     gTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATc                                            >  1:5188932/1‑71 (MQ=255)
                     gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCTGTCGAt                                             <  1:5336978/70‑1 (MQ=255)
                     gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGAt                                             <  1:4590972/70‑1 (MQ=255)
                     gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGAt                                             <  1:2241915/70‑1 (MQ=255)
                     gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGAt                                             <  1:3132333/70‑1 (MQ=255)
                      tACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCt                                                             <  1:1476946/53‑1 (MQ=255)
                       aCGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCg                                          >  1:93619/1‑71 (MQ=255)
                         gAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAcc                                                                 <  1:533759/46‑1 (MQ=255)
                          aaaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGcc                                       >  1:159557/1‑71 (MQ=255)
                           aaaaGACGATATGAACATCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGcgt                                      <  1:5196116/71‑3 (MQ=255)
                           aaaaGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGcgt                                      <  1:1165135/71‑3 (MQ=255)
                           aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGcc                                       <  1:2971644/70‑1 (MQ=255)
                            aaaGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTg                                     <  1:527492/71‑1 (MQ=255)
                                acgacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcggc                                                    <  1:4214709/52‑1 (MQ=255)
                                 cgatATGAACACCTCGGATATGCATCCGTTTATTCATCCGCttt                                                           >  1:796807/1‑42 (MQ=25)
                                 cgatATGAACACCTCGGATATGCATCCGTTTATTCATCCGCttt                                                           >  1:1265174/1‑42 (MQ=25)
                                 cgatATGAACACCTCGGATATGCATCCGTTTATTCATCCGCttt                                                           >  1:2625660/1‑42 (MQ=25)
                                  gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                <  1:2757281/70‑1 (MQ=255)
                                  gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                <  1:673802/70‑1 (MQ=255)
                                  gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                <  1:5605469/70‑1 (MQ=255)
                                   ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcgg                                                     <  1:4542792/46‑1 (MQ=255)
                                   ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcgg                                                     <  1:3226024/46‑1 (MQ=255)
                                   ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcgg                                                     <  1:4017036/46‑1 (MQ=255)
                                   ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcgg                                                     <  1:3953829/46‑1 (MQ=255)
                                   acATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa                                              <  1:295201/55‑1 (MQ=255)
                                    cATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTg                                                            >  1:4451928/1‑40 (MQ=255)
                                     aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcggcg                                                   >  1:1377363/1‑48 (MQ=255)
                                     aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcggcg                                                   >  1:761349/1‑48 (MQ=255)
                                     aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcagcg                                                   >  1:1760693/1‑48 (MQ=255)
                                     aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCg                                               <  1:52347/52‑1 (MQ=37)
                                     aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGa                                              <  1:3625764/53‑1 (MQ=37)
                                      tGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATc                                            <  1:2116227/54‑1 (MQ=255)
                                      tGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGa                             <  1:2070046/69‑1 (MQ=255)
                                       gAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGaa                            >  1:4454508/1‑69 (MQ=255)
                                        aaTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                >  1:1726496/1‑64 (MQ=255)
                                         aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc                        <  1:3035158/71‑1 (MQ=255)
                                            cTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa                                              <  1:2661546/46‑1 (MQ=255)
                                            cTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa                                              <  1:317352/46‑1 (MQ=255)
                                            cTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa                                              <  1:642867/46‑1 (MQ=255)
                                             ctCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGaa                                 <  1:2503919/57‑1 (MQ=255)
                                             ctCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCAGCGTGGGaa                                 <  1:2231905/57‑1 (MQ=255)
                                              tCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATc                                            >  1:3815887/1‑46 (MQ=255)
                                               cGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                  >  1:5187583/1‑71 (MQ=255)
                                                ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTggg                 <  1:1846389/71‑1 (MQ=255)
                                                ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                  >  1:3689816/1‑70 (MQ=255)
                                                   tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGTGaaa              <  1:1481207/71‑1 (MQ=255)
                                                   tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGaaa              <  1:2915590/71‑1 (MQ=255)
                                                   tatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGa                <  1:3322576/69‑1 (MQ=255)
                                                     tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCgg                                         <  1:3744128/42‑1 (MQ=38)
                                                     tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTggg                                   <  1:2517544/48‑1 (MQ=37)
                                                     tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTggg                                   <  1:4323523/48‑1 (MQ=37)
                                                     tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTggg                                   <  1:1551834/48‑1 (MQ=37)
                                                     tGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGc                               <  1:2124322/52‑1 (MQ=255)
                                                      gCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGa                      >  1:2439283/1‑60 (MQ=255)
                                                       cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCt           <  1:3401193/70‑1 (MQ=255)
                                                       cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa          <  1:4377260/71‑1 (MQ=255)
                                                       cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa          <  1:771430/71‑1 (MQ=255)
                                                         tCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc                        <  1:78578/55‑1 (MQ=255)
                                                           cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaaa      <  1:1469837/71‑1 (MQ=255)
                                                           cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACa        >  1:4287315/1‑69 (MQ=255)
                                                           cGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCt                                      >  1:2225820/1‑39 (MQ=255)
                                                           cGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACa        >  1:3522696/1‑69 (MQ=255)
                                                            gTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAAt             <  1:4112380/63‑1 (MQ=255)
                                                             tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGa                             <  1:2210751/46‑1 (MQ=255)
                                                             tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc                        <  1:5515961/51‑1 (MQ=255)
                                                              ttATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAgg    >  1:4407565/1‑69 (MQ=255)
                                                               tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                  >  1:5665279/1‑55 (MQ=255)
                                                               tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                  >  1:5369519/1‑55 (MQ=255)
                                                               tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggta  >  1:4662936/1‑68 (MQ=255)
                                                               tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCa  >  1:3288003/1‑71 (MQ=255)
                                                               tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCa  >  1:5064259/1‑71 (MQ=255)
                                                                     |                                                                
CATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCA  >  minE/833897‑834030

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: