Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,220 T→C 5.5% R1226R (CGT→CGC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,2200TC 5.5% 259.5 / 5.7 90R1226R (CGT→CGCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (38/47);  new base C (3/2);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

CCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAAT  >  minE/835151‑835290
                                                                     |                                                                      
ccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTa                                                                       <  1:230824/71‑1 (MQ=255)
  aTCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGt                                                                        >  1:1865410/1‑68 (MQ=255)
    cTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGa                                                                    >  1:1912363/1‑70 (MQ=255)
     tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTaa                                                                      <  1:1844840/67‑1 (MQ=255)
     tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGAt                                                                   >  1:192605/1‑70 (MQ=255)
       gCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGAt                                                                   <  1:2851865/68‑1 (MQ=255)
       gCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaa                                                                <  1:604381/71‑1 (MQ=255)
        cATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaag                                                               <  1:4072356/71‑1 (MQ=255)
         aTACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaaga                                                              <  1:182280/71‑1 (MQ=255)
            gggCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAc                                                            <  1:2679682/69‑1 (MQ=255)
            cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                                             >  1:601140/1‑69 (MQ=255)
            cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                                             >  1:3651660/1‑69 (MQ=255)
            cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                                                            <  1:2453992/70‑1 (MQ=255)
            cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                                                            <  1:3309420/70‑1 (MQ=255)
            cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                           >  1:1743056/1‑71 (MQ=255)
            cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                           <  1:1816317/71‑1 (MQ=255)
            cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                           >  1:2788969/1‑71 (MQ=255)
             ggCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                                             <  1:232103/68‑1 (MQ=255)
              gCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                                                         <  1:674326/71‑1 (MQ=255)
              gCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGtt                                                         <  1:1072415/71‑1 (MQ=255)
               cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                           >  1:2875070/1‑68 (MQ=255)
               cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTg                                                        <  1:1492233/71‑1 (MQ=255)
                   aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                                    >  1:1524225/1‑71 (MQ=255)
                   aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                                    >  1:1446072/1‑71 (MQ=255)
                   aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                                    >  1:3964302/1‑71 (MQ=255)
                    actactATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTgg                                                   >  1:3977686/1‑71 (MQ=255)
                       actaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                                  >  1:2601759/1‑69 (MQ=255)
                       actaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcg                                                           >  1:761130/1‑59 (MQ=255)
                       actaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcg                                                           >  1:2101860/1‑59 (MQ=255)
                        ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                                  >  1:436350/1‑68 (MQ=255)
                        ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                                  >  1:3308921/1‑68 (MQ=255)
                        ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                                <  1:3464078/70‑1 (MQ=255)
                        ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                                <  1:1400548/70‑1 (MQ=255)
                           tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATga                                                                 <  1:3062234/50‑1 (MQ=255)
                           tGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                                                            >  1:3293592/1‑55 (MQ=255)
                           cggCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATgaagaa                                                             >  1:1825414/2‑54 (MQ=255)
                             gCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                                             >  1:3418545/1‑52 (MQ=255)
                              cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg                                         <  1:834316/71‑1 (MQ=255)
                              cACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg                                         <  1:3299478/71‑1 (MQ=255)
                                ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                           <  1:70176/51‑1 (MQ=255)
                                ccGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATg                                            <  1:2883082/66‑1 (MQ=255)
                                ccGTAGGTTCTAACCGCGATGAGTTGTTGTAGTGCGTAAGATGAAGAACa                                                           <  1:3281758/50‑1 (MQ=255)
                                   tAGGTTCTAACCGCGATGAGTTTGTTGTATTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAgg                                     >  1:4022729/1‑70 (MQ=255)
                                   tAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACt                                                     >  1:2933549/1‑54 (MQ=255)
                                   tAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAgg                                     >  1:233002/1‑70 (MQ=255)
                                      gTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGa                                        <  1:2801162/64‑1 (MQ=255)
                                        tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                               >  1:2742322/1‑71 (MQ=255)
                                        tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                               >  1:2495360/1‑71 (MQ=255)
                                        tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                               >  1:3131654/1‑71 (MQ=255)
                                         cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGa                                              >  1:3203561/1‑55 (MQ=255)
                                         cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                               <  1:4099664/70‑1 (MQ=255)
                                         cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                              >  1:3929115/1‑71 (MQ=255)
                                         cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                              >  1:1286771/1‑71 (MQ=255)
                                          tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                                                         <  1:2065101/43‑1 (MQ=255)
                                          tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACt                                                     <  1:2998402/47‑1 (MQ=255)
                                          tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGc                                    <  1:1360076/64‑1 (MQ=255)
                                          tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGc                                    <  1:1604771/64‑1 (MQ=255)
                                          tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATg                                <  1:3862420/68‑1 (MQ=255)
                                           aaCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                              >  1:1319084/1‑69 (MQ=255)
                                           aaCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                              >  1:1926929/1‑69 (MQ=255)
                                           aaCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                              >  1:3994764/1‑69 (MQ=255)
                                            aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                            <  1:2917359/70‑1 (MQ=255)
                                            aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                            >  1:716929/1‑70 (MQ=255)
                                            aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                            <  1:2785046/70‑1 (MQ=255)
                                            aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAg                           >  1:992980/1‑71 (MQ=255)
                                             ccGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAgg                          <  1:756013/71‑1 (MQ=255)
                                             ccGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAgg                          <  1:1514824/71‑1 (MQ=255)
                                             ccGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAgg                          <  1:2236259/71‑1 (MQ=255)
                                              cgcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCg                                         <  1:3237832/55‑1 (MQ=255)
                                               gcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGAt                                             <  1:1387209/50‑1 (MQ=255)
                                               gcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGaa                                       >  1:2899431/1‑56 (MQ=255)
                                               gcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                            >  1:3286444/1‑67 (MQ=255)
                                               gcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                            >  1:3551702/1‑67 (MQ=255)
                                                cgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTa                        >  1:3874902/1‑70 (MQ=255)
                                                  aTGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAgg                          <  1:792997/66‑1 (MQ=255)
                                                    gAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAAt                                 >  1:4036455/1‑57 (MQ=255)
                                                     aGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGga                   <  1:265130/70‑1 (MQ=255)
                                                       tttGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTa                                                <  1:3866701/39‑1 (MQ=255)
                                                       tttGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAg                                               <  1:529778/40‑1 (MQ=255)
                                                       tttGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGagag                <  1:1979568/71‑1 (MQ=255)
                                                         tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCg              >  1:888426/1‑71 (MQ=255)
                                                             gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa           <  1:1597182/70‑1 (MQ=255)
                                                             gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa           <  1:4134512/70‑1 (MQ=255)
                                                             gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa           <  1:4080190/70‑1 (MQ=255)
                                                              tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaaaa         <  1:804020/71‑1 (MQ=255)
                                                              tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa           <  1:3671209/69‑1 (MQ=255)
                                                                gTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCa                                   <  1:371431/43‑1 (MQ=255)
                                                                gTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAAt        >  1:981673/1‑70 (MQ=255)
                                                                gTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATg       <  1:323748/71‑1 (MQ=255)
                                                                  gCGTAAGATGAAGAACATTGACTGGTTAGATGGCg                                         <  1:3822226/35‑1 (MQ=255)
                                                                   cGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAgg                          >  1:830165/1‑49 (MQ=255)
                                                                     tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAAt  <  1:1762831/71‑1 (MQ=255)
                                                                     tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAAt  <  1:2971032/71‑1 (MQ=255)
                                                                     |                                                                      
CCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAAT  >  minE/835151‑835290

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: