Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,209 G→A 5.7% V1223I (GTT→ATT) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,2090GA 5.7% 291.7 / 5.0 88V1223I (GTT→ATT) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (44/39);  new base A (3/2);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 110 total aligned reads displayed.

CGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG  >  minE/835142‑835278
                                                                   |                                                                     
cggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGTTTCTAACCGCGATGAGTttgttg                                                                    <  1:245151/71‑1 (MQ=255)
cggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                                    <  1:2732991/71‑1 (MQ=255)
cggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGGttgttg                                                                    <  1:6775/71‑1 (MQ=255)
 ggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                                    >  1:3279475/1‑70 (MQ=255)
 ggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgt                                                                      >  1:1233875/1‑68 (MQ=255)
 ggcTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                                                   <  1:3688752/71‑1 (MQ=255)
    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)
            tCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttg                                                                       <  1:2799776/56‑1 (MQ=255)
            tCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttg                                                                       <  1:1606182/56‑1 (MQ=255)
            tCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCg                                                             >  1:3939148/1‑66 (MQ=255)
             cTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGa                                                        >  1:1912363/1‑70 (MQ=255)
              tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                                    <  1:1387017/57‑1 (MQ=255)
              tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                                                   <  1:2129653/58‑1 (MQ=255)
              tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                                  >  1:1560120/1‑59 (MQ=255)
              tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                                  >  1:2168586/1‑59 (MQ=255)
              tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTaa                                                          <  1:1844840/67‑1 (MQ=255)
              tGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGAt                                                       >  1:192605/1‑70 (MQ=255)
               ggCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                                    <  1:3085191/56‑1 (MQ=255)
                gCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGAt                                                       <  1:2851865/68‑1 (MQ=255)
                gCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaa                                                    <  1:604381/71‑1 (MQ=255)
                 cATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaag                                                   <  1:4072356/71‑1 (MQ=255)
                   tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                                    <  1:1714691/52‑1 (MQ=255)
                   tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                                    <  1:2334095/52‑1 (MQ=255)
                   tACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgttg                                                                    <  1:1260439/52‑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:2788969/1‑71 (MQ=255)
                     cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                               <  1:1816317/71‑1 (MQ=255)
                     cGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACa                                               >  1:1743056/1‑71 (MQ=255)
                      ggCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATgaagaa                                                 <  1:232103/68‑1 (MQ=255)
                       gCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttg                                                                       >  1:1984950/1‑45 (MQ=255)
                       gCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACAtt                                             <  1:674326/71‑1 (MQ=255)
                       gCTTTAACTACTACGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAACGtt                                             <  1:1072415/71‑1 (MQ=255)
                        cTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTg                                            <  1:1492233/71‑1 (MQ=255)
                         tttAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                                                   <  1:2800791/47‑1 (MQ=255)
                            aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgt                                                                      >  1:157987/1‑41 (MQ=255)
                            aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTttgt                                                                      >  1:3011530/1‑41 (MQ=255)
                            aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTa                                                                  >  1:584510/1‑45 (MQ=255)
                            aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                        >  1:1524225/1‑71 (MQ=255)
                            aaCTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTg                                        >  1:3964302/1‑71 (MQ=255)
                              ctactaTGGCACCGTAGGTTCTAACCGCGATGAGTTtgttgt                                                                   <  1:1107993/42‑1 (MQ=255)
                                actaCGGCACCGTCGGATCGAACCGCGATGAGTTCATCATGATCCGCAAGATGAAGAAcg                                               >  1:761130/1‑59 (MQ=255)
                                 ctaTGGCACCGTAGGTTCTAACCGCGATGAGTttgt                                                                      >  1:3940512/1‑36 (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:3131654/1‑71 (MQ=255)
                                                 tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                   >  1:2495360/1‑71 (MQ=255)
                                                 tCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                   >  1:2742322/1‑71 (MQ=255)
                                                  cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGa                   <  1:4099664/70‑1 (MQ=255)
                                                  cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                  >  1:1286771/1‑71 (MQ=255)
                                                  cTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                  >  1:3929115/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:3994764/1‑69 (MQ=255)
                                                    aaCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                  >  1:1926929/1‑69 (MQ=255)
                                                    aaCCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGAc                  >  1:1319084/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)
                                                      ccGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAgg              <  1:2236259/71‑1 (MQ=255)
                                                      ccGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAgg              <  1:756013/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)
                                                         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)
                                                                   |                                                                     
CGGCTATGCCCATCTGGCATACGGCTTTAACTACTATGGCACCGTAGGTTCTAACCGCGATGAGTTTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG  >  minE/835142‑835278

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: