Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,215 G→A 5.5% V1225M (GTG→ATG) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,2150GA 5.5% 300.2 / 5.0 91V1225M (GTG→ATG) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (41/44);  new base A (3/2);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

TATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATG  >  minE/835146‑835285
                                                                     |                                                                      
tATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGt                                                                       >  1:1902457/1‑71 (MQ=255)
tATGCCCACCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGt                                                                       >  1:3935296/1‑71 (MQ=255)
 aTGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAAGCGCGATGAGTTTGTTGTAGTg                                                                      >  1:71732/1‑71 (MQ=255)
 aTGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTg                                                                      >  1:3152808/1‑71 (MQ=255)
     ccATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTa                                                                  <  1:230824/71‑1 (MQ=255)
       aTCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGt                                                                   >  1:1865410/1‑68 (MQ=255)
        tCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCg                                                                    >  1:3939148/1‑66 (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:3309420/70‑1 (MQ=255)
                 cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAAc                                                       <  1:2453992/70‑1 (MQ=255)
                 cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                      <  1:1816317/71‑1 (MQ=255)
                 cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                                      >  1:1743056/1‑71 (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:3964302/1‑71 (MQ=255)
                        aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                               >  1:1446072/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:3308921/1‑68 (MQ=255)
                             ctaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                             >  1:436350/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:3131654/1‑71 (MQ=255)
                                             tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                          >  1:2742322/1‑71 (MQ=255)
                                             tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                          >  1:2495360/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:716929/1‑70 (MQ=255)
                                                 aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                       <  1:2917359/70‑1 (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:2236259/71‑1 (MQ=255)
                                                  ccGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAgg                     <  1:1514824/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)
                                                                     |                                                                      
TATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATG  >  minE/835146‑835285

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: