Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,858,220 A→G 14.1% F454F (TTT→TTC cysI ← sulfite reductase, beta subunit, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,858,2200AG14.1% 170.0 / 14.7 98F454F (TTT→TTCcysIsulfite reductase, beta subunit, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (68/16);  new base G (8/6);  total (76/22)
Fisher's exact test for biased strand distribution p-value = 7.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 101 total aligned reads displayed.

ACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCA  >  minE/1858155‑1858287
                                                                 |                                                                   
acacGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATaaaa                                                                  >  1:1022499/1‑69 (MQ=255)
acacGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATaaaa                                                                  >  1:2898307/1‑69 (MQ=255)
     cATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGAc                                                               >  1:969183/1‑67 (MQ=255)
     cATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGAc                                                               >  1:655431/1‑67 (MQ=255)
     cATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGAc                                                               >  1:436653/1‑67 (MQ=255)
     cATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGAc                                                               >  1:3106872/1‑67 (MQ=255)
     cATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGAc                                                               >  1:1013922/1‑67 (MQ=255)
     cATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGAc                                                               >  1:1529955/1‑67 (MQ=255)
     cATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGAc                                                               >  1:1267921/1‑67 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATa                                                                     >  1:2062211/1‑57 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCa                                                           >  1:2640986/1‑67 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:2358620/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:1557832/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:169640/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:1715641/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:1794595/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:1803687/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:187219/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:1993056/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:2220454/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:1164236/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:2235887/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:2588363/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:308484/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:3103606/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:1086609/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:3397217/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:382739/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:594556/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:1018447/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:1321720/1‑68 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAg                                                          >  1:1222489/1‑68 (MQ=255)
              atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATGAAAGACGGCAGGAAAc                                                     >  1:2384524/1‑68 (MQ=255)
              atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATGAAAGACGGCAGGAAAc                                                     >  1:489993/1‑68 (MQ=255)
              atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATGAAAGACGGCAGGAAACg                                                    >  1:2307240/1‑69 (MQ=255)
              atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATGAAAGACGGCAGGAAACg                                                    >  1:643944/1‑69 (MQ=255)
              atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATGAAAGACGGCAGGAAACg                                                    >  1:638083/1‑69 (MQ=255)
              atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATGAAAGACGGCAGGAAACg                                                    >  1:1312496/1‑69 (MQ=255)
              atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATGAAAGACGGCAGGAAACg                                                    >  1:622181/1‑69 (MQ=255)
              atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATGAAAGACGGCAGGAAACg                                                    >  1:1064851/1‑69 (MQ=255)
                 cTCATCGCTGCCACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGctc                                                 >  1:1731873/1‑69 (MQ=255)
                 cTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGctc                                                 >  1:1392550/1‑69 (MQ=255)
                 cTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGctc                                                 >  1:493213/1‑69 (MQ=255)
                 cTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGctc                                                 >  1:893251/1‑69 (MQ=255)
                 cTCATCGCTGACACCATGTTTCGCCATTAAATTACCGATGTTGTCGATAAAAGACGGCAGGAAACGctc                                                 >  1:901196/1‑69 (MQ=255)
                      cGCTGTCACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCt                                             <  1:1186873/68‑1 (MQ=255)
                      cGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCt                                             <  1:138217/68‑1 (MQ=255)
                      cGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCt                                             <  1:2406365/68‑1 (MQ=255)
                      cGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCt                                             <  1:1288389/68‑1 (MQ=255)
                       gCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACgg                                                             >  1:2970265/1‑51 (MQ=255)
                       gCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGc                                                            >  1:3116001/1‑52 (MQ=255)
                       gCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGc                                                            >  1:2226455/1‑52 (MQ=255)
                       gCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTc                                           >  1:954613/1‑69 (MQ=255)
                       gCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTc                                           >  1:2551486/1‑69 (MQ=255)
                       gCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTc                                           >  1:2923914/1‑69 (MQ=255)
                       gCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTc                                           >  1:2477825/1‑69 (MQ=255)
                       gCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTc                                           >  1:2807300/1‑69 (MQ=255)
                       gCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTc                                           >  1:3340825/1‑69 (MQ=255)
                       gCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTc                                           >  1:1908716/1‑69 (MQ=255)
                       gCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTc                                           >  1:1644307/1‑69 (MQ=255)
                                aTGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATc                                    >  1:3060024/1‑67 (MQ=255)
                                aTGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATc                                    >  1:467131/1‑67 (MQ=255)
                                aTGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATc                                    >  1:1940895/1‑67 (MQ=255)
                                aTGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATc                                    >  1:1483430/1‑67 (MQ=255)
                                aTGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATc                                    >  1:1523950/1‑67 (MQ=255)
                                          aTTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCaa                        <  1:1150830/69‑1 (MQ=255)
                                          aTTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCaa                        <  1:1547264/69‑1 (MQ=255)
                                          aTTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCaa                        <  1:2789868/69‑1 (MQ=255)
                                           ttAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCaa                        <  1:1201493/68‑1 (MQ=255)
                                           ttAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCaa                        <  1:1168807/68‑1 (MQ=255)
                                           ttAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCaa                        <  1:1613866/68‑1 (MQ=255)
                                           ttAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCaa                        <  1:2357077/68‑1 (MQ=255)
                                             aaaTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTc                     >  1:2699017/1‑69 (MQ=255)
                                             aaaTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTc                     >  1:53843/1‑69 (MQ=255)
                                             aaaTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTc                     >  1:1882097/1‑69 (MQ=255)
                                             aaaTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTc                     >  1:879571/1‑69 (MQ=255)
                                             aaaTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTc                     >  1:1209628/1‑69 (MQ=255)
                                             aaaTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTc                     >  1:1110587/1‑69 (MQ=255)
                                             aaaTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTc                     >  1:3025582/1‑69 (MQ=255)
                                                ttATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCggg                  <  1:1225488/69‑1 (MQ=255)
                                                ttATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCggg                  <  1:110167/69‑1 (MQ=255)
                                                ttATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCggg                  <  1:2160916/69‑1 (MQ=255)
                                                ttATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCggg                  <  1:1019604/69‑1 (MQ=255)
                                                  aTCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGa                 >  1:2052327/1‑68 (MQ=255)
                                                     gATGTTGTCGATGAAAGACGTCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGa             <  1:1340396/69‑1 (MQ=255)
                                                     gATGTTGTCGATGAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGa             <  1:309711/69‑1 (MQ=255)
                                                     gATGTTGTCGATGAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGa             <  1:155428/69‑1 (MQ=255)
                                                     gATGTTGTCGATGAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGa             <  1:1629765/69‑1 (MQ=255)
                                                     gATGTTGTCGATGAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGa             <  1:2034935/69‑1 (MQ=255)
                                                     gATGTTGTCGATGAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGa             <  1:1454831/69‑1 (MQ=255)
                                                         ttgtCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCggg                  <  1:2901769/60‑1 (MQ=255)
                                                                tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCTAGTCGGGAATGACACGCAAGcc   >  1:2490360/1‑68 (MQ=255)
                                                                tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCa       >  1:2900433/1‑64 (MQ=255)
                                                                tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa  >  1:1719343/1‑69 (MQ=255)
                                                                tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa  >  1:1714941/1‑69 (MQ=255)
                                                                tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa  >  1:1137864/1‑69 (MQ=255)
                                                                tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa  >  1:1733425/1‑69 (MQ=255)
                                                                tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa  >  1:477983/1‑69 (MQ=255)
                                                                tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa  >  1:2196522/1‑69 (MQ=255)
                                                                tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa  >  1:1372485/1‑69 (MQ=255)
                                                                 |                                                                   
ACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCA  >  minE/1858155‑1858287

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: