Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,858,847 C→T 18.6% K245K (AAG→AAA cysI ← sulfite reductase, beta subunit, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,858,8470CT18.6% 130.0 / 21.4 59K245K (AAG→AAAcysIsulfite reductase, beta subunit, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (23/25);  new base T (8/3);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGG  >  Exported/1858781‑1858915
                                                                  |                                                                    
gTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGc                                                                  >  1:4558323/1‑71 (MQ=255)
      tttCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGtttt                                                            >  1:4877252/1‑71 (MQ=255)
             ttgCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGa                                                      <  1:2507445/70‑1 (MQ=255)
             ttgCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGa                                                      <  1:60340/70‑1 (MQ=255)
             ttgCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGa                                                      <  1:1665613/70‑1 (MQ=255)
                   tgtTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCAc                                               >  1:382797/1‑71 (MQ=255)
                     ttCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACg                                              <  1:174326/70‑1 (MQ=255)
                     ttCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGa                                             <  1:65929/71‑1 (MQ=255)
                     ttCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACAAGCTTGCCGtttt                                                            >  1:4569764/1‑56 (MQ=255)
                        gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAg                                           <  1:3141286/70‑1 (MQ=255)
                        gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGt                                          <  1:963666/71‑1 (MQ=255)
                        gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGt                                          <  1:4121802/71‑1 (MQ=255)
                        gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGt                                          <  1:2186528/71‑1 (MQ=255)
                        gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGt                                          <  1:2062524/71‑1 (MQ=255)
                        gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGt                                          <  1:1304954/71‑1 (MQ=255)
                         aTGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCg                                                       <  1:2324579/57‑1 (MQ=255)
                         aTGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCg                                                       <  1:5001757/57‑1 (MQ=255)
                         aTGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCg                                                       <  1:507387/57‑1 (MQ=255)
                         aTGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCg                                                       <  1:1492520/57‑1 (MQ=255)
                          tGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTc                                        <  1:1632057/71‑1 (MQ=255)
                                 cccaccgcccaccAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGt                                 <  1:1493628/71‑1 (MQ=255)
                                  ccaccgcccaccAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGtt                                <  1:3715132/71‑1 (MQ=255)
                                       gcccaccAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCG‑‑GATCGCCACGAAGTTCATGTCGTTGGCGTg                          >  1:4621927/1‑70 (MQ=255)
                                          caccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGtt                                         <  1:4064512/54‑1 (MQ=255)
                                          caccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGtt                                         <  1:2739174/54‑1 (MQ=255)
                                          caccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGtt                                         <  1:4186436/54‑1 (MQ=255)
                                          caccaACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTg                          >  1:2119353/1‑69 (MQ=255)
                                            ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGc                                                        >  1:4226924/1‑37 (MQ=255)
                                            ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa                                       >  1:5020576/1‑54 (MQ=255)
                                            ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa                                       >  1:3137108/1‑54 (MQ=255)
                                            ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa                                       >  1:2957937/1‑54 (MQ=255)
                                            ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa                                       >  1:2921073/1‑54 (MQ=255)
                                            ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa                                       >  1:4070660/1‑54 (MQ=255)
                                            ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa                                       >  1:4661521/1‑54 (MQ=255)
                                            ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa                                       >  1:1522611/1‑54 (MQ=255)
                                             caACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGtt                                         <  1:565259/51‑1 (MQ=255)
                                                  ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata                >  1:1275562/1‑71 (MQ=255)
                                                  ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata                >  1:1075576/1‑71 (MQ=255)
                                                  ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata                >  1:135824/1‑71 (MQ=255)
                                                  ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata                >  1:3186400/1‑71 (MQ=255)
                                                  ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata                >  1:605313/1‑71 (MQ=255)
                                                  ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata                >  1:5055894/1‑71 (MQ=255)
                                                  ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGa                  >  1:3245145/1‑69 (MQ=255)
                                                  ggTTAAAGCCCACCAGCTTGCCGTTTTCGGAGATCCCCACGAAGTTCATGTc                                   >  1:631277/1‑52 (MQ=39)
                                                   gTTAAAGCCCACCAGCTTGCCTTTTTCGGCGATCTCCACGAAGTTCATGTCGTTGGCGTGCAGa                      <  1:511506/64‑1 (MQ=255)
                                                   gTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGa                      <  1:3179119/64‑1 (MQ=255)
                                                   gTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGa                      <  1:3372639/64‑1 (MQ=255)
                                                   gTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGa                      <  1:2891826/64‑1 (MQ=255)
                                                     tAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATc                    >  1:3394163/1‑64 (MQ=255)
                                                      aaaGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata                >  1:1441128/1‑67 (MQ=255)
                                                      aaaGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGt            >  1:2663674/1‑71 (MQ=255)
                                                      aaaGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGt            >  1:2357715/1‑71 (MQ=255)
                                                       aaGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGt            <  1:232850/70‑1 (MQ=255)
                                                        aGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTc          <  1:1152356/71‑1 (MQ=255)
                                                             accaGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTgg     >  1:2798108/1‑71 (MQ=255)
                                                               caGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCg   >  1:1731380/1‑71 (MQ=255)
                                                               caGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCg   >  1:1355055/1‑71 (MQ=255)
                                                               caGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCg   >  1:95249/1‑71 (MQ=255)
                                                                aGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTATGTGGCgg  >  1:3220066/1‑71 (MQ=255)
                                                                  cTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGa                  <  1:3773389/53‑1 (MQ=255)
                                                                  cTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGa                  <  1:3059993/53‑1 (MQ=255)
                                                                  |                                                                    
GTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGG  >  Exported/1858781‑1858915

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