Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,844,959 A→G 7.8% T42T (ACT→ACC cysN ← sulfate adenylyltransferase, subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,844,9590AG 7.8% 423.7 / 18.7 141T42T (ACT→ACCcysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base A (48/82);  new base G (5/6);  total (53/88)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 152 total aligned reads displayed.

TCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGT  >  Exported/1844891‑1845019
                                                                    |                                                            
     aTGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACtt                                                        <  1:4979159/70‑1 (MQ=255)
     aTGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACtt                                                        <  1:2103348/70‑1 (MQ=255)
     aTGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACtt                                                        <  1:3422146/70‑1 (MQ=255)
      tGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACtttt                                                      <  1:4224143/71‑1 (MQ=255)
      tGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACtttt                                                      <  1:3735925/71‑1 (MQ=255)
         aGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGcc                                                   >  1:1387113/1‑71 (MQ=255)
          gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg                                                  >  1:2317683/1‑71 (MQ=255)
          gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg                                                  <  1:3127250/71‑1 (MQ=255)
          gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg                                                  >  1:3452815/1‑71 (MQ=255)
          gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg                                                  <  1:900538/71‑1 (MQ=255)
          gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg                                                  <  1:3530537/71‑1 (MQ=255)
           cGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg                                                  <  1:1092429/70‑1 (MQ=255)
            gATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCagag                                                             >  1:5180144/1‑58 (MQ=255)
            gATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCagag                                                             >  1:3073457/1‑58 (MQ=255)
              tGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCAt                                              >  1:2819545/1‑71 (MQ=255)
              tGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCAt                                              >  1:599302/1‑71 (MQ=255)
              tGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCAt                                              >  1:4957516/1‑71 (MQ=255)
              tGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCAt                                              >  1:1250219/1‑71 (MQ=255)
               gagaGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACttt                                                       >  1:1749413/1‑61 (MQ=255)
                agagCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                             <  1:1014456/70‑1 (MQ=255)
                 gagCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACGAATCAGAGTACTTTTGCCGTCATCGa                                           <  1:4759899/71‑1 (MQ=255)
                   gCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCATCGACg                                         <  1:2805181/71‑1 (MQ=255)
                   gCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCATCGACg                                         <  1:1841782/71‑1 (MQ=255)
                    cTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCATCGACg                                         <  1:4455371/70‑1 (MQ=255)
                      gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                             <  1:204410/64‑1 (MQ=255)
                      gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                             <  1:3250941/64‑1 (MQ=255)
                      gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                             <  1:3786551/64‑1 (MQ=255)
                      gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                             <  1:1325347/64‑1 (MQ=255)
                      gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                             <  1:4663480/64‑1 (MQ=255)
                      gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                             <  1:1598814/64‑1 (MQ=255)
                      gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                             <  1:1477745/64‑1 (MQ=255)
                       aTCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                             <  1:4083309/63‑1 (MQ=255)
                       aTCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                             <  1:4515000/63‑1 (MQ=255)
                       aTCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCt                                       <  1:3101550/69‑1 (MQ=255)
                       aTCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCt                                       <  1:2609820/69‑1 (MQ=255)
                       aTCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                     >  1:4406642/1‑71 (MQ=255)
                         cTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                         >  1:2652145/1‑65 (MQ=255)
                         cTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                         >  1:5151043/1‑65 (MQ=255)
                         cTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                         >  1:5156940/1‑65 (MQ=255)
                         cTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                         >  1:1795352/1‑65 (MQ=255)
                         cTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                         >  1:3572549/1‑65 (MQ=255)
                          ttCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                             <  1:4399680/60‑1 (MQ=255)
                            cGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                           >  1:4125950/1‑60 (MQ=255)
                            cGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                           >  1:795834/1‑60 (MQ=255)
                            cGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                           >  1:2190321/1‑60 (MQ=255)
                            cGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                           >  1:182027/1‑60 (MQ=255)
                             gTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCagag                                                             <  1:913811/41‑1 (MQ=255)
                             gTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCagag                                                             <  1:3778932/41‑1 (MQ=255)
                                gATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGt                              >  1:3487178/1‑69 (MQ=255)
                                gATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGt                              >  1:1639488/1‑69 (MQ=255)
                                gATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGt                              >  1:2157243/1‑69 (MQ=255)
                                   ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                           >  1:1329625/1‑53 (MQ=255)
                                   ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                           >  1:1307302/1‑53 (MQ=255)
                                   ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                           >  1:3774012/1‑53 (MQ=255)
                                   ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                           >  1:4718777/1‑53 (MQ=255)
                                   ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                           >  1:825490/1‑53 (MQ=255)
                                   ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                           >  1:2305924/1‑53 (MQ=255)
                                   ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGa                          >  1:1160425/1‑70 (MQ=255)
                                   ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGa                          >  1:4343903/1‑70 (MQ=255)
                                   ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGa                          >  1:3495973/1‑70 (MQ=255)
                                   ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGa                          >  1:698746/1‑70 (MQ=255)
                                    tGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGaa                         <  1:3812067/70‑1 (MQ=255)
                                       cGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCATAAAACg                     <  1:11128/71‑1 (MQ=255)
                                       cGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACg                     <  1:1174046/71‑1 (MQ=255)
                                       cGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACg                     <  1:3177612/71‑1 (MQ=255)
                                        gggTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACg                     >  1:506740/1‑70 (MQ=255)
                                         ggTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                             <  1:118449/45‑1 (MQ=255)
                                            aTCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCt                                       <  1:3195889/48‑1 (MQ=255)
                                             tCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                     <  1:83677/49‑1 (MQ=255)
                                             tCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                     <  1:94083/49‑1 (MQ=255)
                                             tCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                     <  1:2362473/49‑1 (MQ=255)
                                             tCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                     <  1:3597245/49‑1 (MQ=255)
                                             tCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                     <  1:1478242/49‑1 (MQ=255)
                                               gTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCa                            >  1:4280640/1‑56 (MQ=255)
                                               gTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCa                            >  1:2426314/1‑56 (MQ=255)
                                               gTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCa                            >  1:2295246/1‑56 (MQ=255)
                                               gTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCa                            >  1:3432887/1‑56 (MQ=255)
                                                 gcagcaGGCGACCAATCAGGGTACTTTTGCCGTCATCGACg                                         <  1:3627449/41‑1 (MQ=37)
                                                       gacgacCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTAt       <  1:4600178/69‑1 (MQ=255)
                                                         cgacCAATCTGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTAtgttg   <  1:2105752/71‑1 (MQ=255)
                                                         cgacCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTAtgttg   <  1:3054453/71‑1 (MQ=255)
                                                         cgacCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTAtgttg   <  1:754378/71‑1 (MQ=255)
                                                         cgacCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTAtgttg   <  1:234487/71‑1 (MQ=255)
                                                          gacCAATCAGAGTACTTTTGCCGTCATCGACGCTACcac                                  <  1:4632720/39‑1 (MQ=255)
                                                          gacCAATCAGAGTACTTTTGCCGTCATCGACGCTACcac                                  <  1:3279101/39‑1 (MQ=255)
                                                          gacCAATCAGAGTACTTTTGCCGTCATCGACGCTACcac                                  <  1:500828/39‑1 (MQ=255)
                                                          gacCAATCAGAGTACTTTTGCCGTCATCGACGCTACcac                                  <  1:2519273/39‑1 (MQ=255)
                                                          gacCAATCAGAGTACTTTTGCCGTCATCGACGCTACcac                                  <  1:2485824/39‑1 (MQ=255)
                                                          gacCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGcagcag               <  1:3520098/58‑1 (MQ=255)
                                                          gacCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGcagcag               <  1:4854976/58‑1 (MQ=255)
                                                          gacCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGcagcag               <  1:2620601/58‑1 (MQ=255)
                                                          gacCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTAtgttgt  >  1:4011021/1‑71 (MQ=255)
                                                          gacCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTAtgttgt  >  1:5126619/1‑71 (MQ=255)
                                                          gacCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTAtgttgt  >  1:2866789/1‑71 (MQ=255)
                                                           acCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAgg                               <  1:485576/41‑1 (MQ=255)
                                                           acCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAgg                               <  1:201902/41‑1 (MQ=255)
                                                           acCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAgg                               <  1:1763885/41‑1 (MQ=255)
                                                           acCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAgg                               <  1:4739633/41‑1 (MQ=255)
                                                             cAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCttt          <  1:4788164/60‑1 (MQ=255)
                                                             cAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCttt          <  1:1247457/60‑1 (MQ=255)
                                                                    |                                                            
TCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGT  >  Exported/1844891‑1845019

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: