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

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

TATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAA  >  minE/835146‑835283
                                                                    |                                                                     
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)
          tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                                       >  1:2168586/1‑59 (MQ=255)
          tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                                       >  1:1560120/1‑59 (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: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)
                        aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                                       >  1:584510/1‑45 (MQ=255)
                        aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                             >  1:1446072/1‑71 (MQ=255)
                        aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                             >  1:3964302/1‑71 (MQ=255)
                        aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                             >  1:1524225/1‑71 (MQ=255)
                         actactATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTgg                                            >  1:3977686/1‑71 (MQ=255)
                            actaTGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGt                                           >  1:2601759/1‑69 (MQ=255)
                            actaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcg                                                    >  1:2101860/1‑59 (MQ=255)
                            actaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcg                                                    >  1:761130/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:2495360/1‑71 (MQ=255)
                                             tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                        >  1:2742322/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:1604771/64‑1 (MQ=255)
                                               tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGc                             <  1:1360076/64‑1 (MQ=255)
                                               tAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATg                         <  1:3862420/68‑1 (MQ=255)
                                                aaCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                       >  1:1926929/1‑69 (MQ=255)
                                                aaCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                       >  1:1319084/1‑69 (MQ=255)
                                                aaCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                       >  1:3994764/1‑69 (MQ=255)
                                                 aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                     <  1:2785046/70‑1 (MQ=255)
                                                 aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                     >  1:716929/1‑70 (MQ=255)
                                                 aCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                     <  1:2917359/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:3551702/1‑67 (MQ=255)
                                                    gcgATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCa                     >  1:3286444/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:4134512/70‑1 (MQ=255)
                                                                  gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTaa    <  1:1597182/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)
                                                                    |                                                                     
TATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAA  >  minE/835146‑835283

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: