Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,903,259 T→A 46.3% R262R (CGA→CGT nrdD ← anaerobic ribonucleoside‑triphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,903,2590TA46.3% 19.9 / 141.6 121R262R (CGA→CGTnrdDanaerobic ribonucleoside‑triphosphate reductase
Reads supporting (aligned to +/- strand):  ref base T (29/36);  new base A (41/15);  total (70/51)
Fisher's exact test for biased strand distribution p-value = 1.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 129 total aligned reads displayed.

GTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGATGTTTGTTGTA  >  minE/2903206‑2903327
                                                     |                                                                    
gtgtTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAg                                                       >  1:2988975/1‑69 (MQ=255)
 tgtTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGc                                                      >  1:1711875/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGTTACGAGAGTTAGCGTAGcc                                                     >  1:826740/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:3383766/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:2428882/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:2402322/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:266952/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:2855595/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:2887426/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:16131/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:3013902/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:3287095/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:3305356/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:422011/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:2230939/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:575154/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:701712/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:1219124/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:791284/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:851720/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:852117/1‑69 (MQ=255)
  gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:924106/1‑69 (MQ=255)
  gtTTACTTCGTACTACAGCGACTGGAAGGCATCGTAGCACTCTTTGATCGTACGAGAGTTAGCGTAGcc                                                     >  1:723793/1‑69 (MQ=255)
  gtTTACTTCGAACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                                                     >  1:3147553/1‑69 (MQ=255)
     tACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTc                                                                     >  1:2921489/1‑50 (MQ=255)
             aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc                                           <  1:1937609/68‑1 (MQ=255)
             aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc                                           <  1:1038472/68‑1 (MQ=255)
             aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc                                           <  1:2576157/68‑1 (MQ=255)
             aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc                                           <  1:483424/68‑1 (MQ=255)
             aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc                                           <  1:953349/68‑1 (MQ=255)
             aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc                                           <  1:1161713/68‑1 (MQ=255)
             aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc                                           <  1:12670/68‑1 (MQ=255)
             aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc                                           <  1:126753/68‑1 (MQ=255)
             aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc                                           <  1:1349412/68‑1 (MQ=255)
             aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc                                           <  1:1862862/68‑1 (MQ=255)
             aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc                                           <  1:1659195/68‑1 (MQ=255)
             aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc                                           <  1:1708827/68‑1 (MQ=255)
              cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg                                         >  1:1376271/1‑69 (MQ=255)
              cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg                                         >  1:1028205/1‑69 (MQ=255)
              cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg                                         >  1:863879/1‑69 (MQ=255)
              cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg                                         >  1:641670/1‑69 (MQ=255)
              cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg                                         >  1:1434078/1‑69 (MQ=255)
              cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg                                         >  1:1683010/1‑69 (MQ=255)
              cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg                                         >  1:2734640/1‑69 (MQ=255)
              cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg                                         >  1:2545112/1‑69 (MQ=255)
              cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg                                         >  1:1873364/1‑69 (MQ=255)
                ccAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGa                                       <  1:494129/69‑1 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGc                                                      <  1:3160495/49‑1 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGc                                                      <  1:339627/49‑1 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGc                                                      <  1:2015213/49‑1 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGc                                                      <  1:2311643/49‑1 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGc                                                      <  1:1172598/49‑1 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGc                                                      <  1:847428/49‑1 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGc                                                      <  1:2585359/49‑1 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGtt                                   >  1:1855437/1‑68 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTc                                  >  1:821371/1‑69 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTc                                  >  1:555539/1‑69 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTc                                  >  1:492272/1‑69 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTc                                  >  1:469807/1‑69 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTc                                  >  1:2092518/1‑69 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTc                                  >  1:1126952/1‑69 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTc                                  >  1:2555338/1‑69 (MQ=255)
                     gACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTCC‑‑GATGTTCCa                                >  1:1355188/1‑69 (MQ=255)
                      aCTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTcc                                 >  1:2268286/1‑69 (MQ=255)
                            aGGCATCGTAGCACTCTTTTATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCActc                             <  1:442067/67‑1 (MQ=255)
                            aGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCActc                             <  1:2427500/67‑1 (MQ=255)
                            aGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCActc                             <  1:1754781/67‑1 (MQ=255)
                            aGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCActc                             <  1:640643/67‑1 (MQ=255)
                            aGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCActc                             <  1:34782/67‑1 (MQ=255)
                             ggCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGc                                                      <  1:2679703/41‑1 (MQ=255)
                                   gTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGt                    >  1:1439006/1‑69 (MQ=255)
                                   gTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGt                    >  1:3422291/1‑69 (MQ=255)
                                   gTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGt                    >  1:3312357/1‑69 (MQ=255)
                                   gTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGt                    >  1:600151/1‑69 (MQ=255)
                                      gCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGTGGTCtt                 <  1:5337/69‑1 (MQ=255)
                                      gCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCtt                 <  1:3050542/69‑1 (MQ=255)
                                      gCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCtt                 <  1:2881448/69‑1 (MQ=255)
                                      gCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCtt                 <  1:3220564/69‑1 (MQ=255)
                                      gCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCtt                 <  1:1925689/69‑1 (MQ=255)
                                      gCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCtt                 <  1:1096237/69‑1 (MQ=255)
                                      gCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCtt                 <  1:133245/69‑1 (MQ=255)
                                      gCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCtt                 <  1:2975840/69‑1 (MQ=255)
                                      gCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCtt                 <  1:3342761/69‑1 (MQ=255)
                                       cACTCTTTTATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc                                           <  1:1429397/42‑1 (MQ=255)
                                        aCTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGc               <  1:1573006/69‑1 (MQ=255)
                                        aCTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCAGGGATGTTCCACTCTTCTGCGGTCTTGc               <  1:2068620/69‑1 (MQ=255)
                                           ctTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGAt            >  1:340467/1‑69 (MQ=255)
                                           ctTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGAt            >  1:1258251/1‑69 (MQ=255)
                                           ctTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGAt            >  1:2975322/1‑69 (MQ=255)
                                            tttGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCTGTCTTGCGAtg           >  1:697266/1‑69 (MQ=255)
                                            tttGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGAtg           >  1:1968272/1‑69 (MQ=255)
                                            tttGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGAtg           >  1:804218/1‑69 (MQ=255)
                                                 tGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGAtgtttg       <  1:60121/68‑1 (MQ=255)
                                                     tCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGATGTTTGTTGTa  <  1:1860895/69‑1 (MQ=255)
                                                     tCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGATGTTTGTTGTa  <  1:176522/69‑1 (MQ=255)
                                                     tCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGATGTTTGTTGTa  <  1:1774887/69‑1 (MQ=255)
                                                     tCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGATGTTTGTTGTa  <  1:1182091/69‑1 (MQ=255)
                                                     tCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGATGTTTGTTGTa  <  1:2293286/69‑1 (MQ=255)
                                                     tCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGATGTTTGTTGTa  <  1:494339/69‑1 (MQ=255)
                                                     tCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGATGTTTGTTGTa  <  1:1629195/69‑1 (MQ=255)
                                                     |                                                                    
GTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGATGTTCCACTCTTCTGCGGTCTTGCGATGTTTGTTGTA  >  minE/2903206‑2903327

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: