Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,844,947 A→G 8.5% R46R (CGT→CGC cysN ← sulfate adenylyltransferase, subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,844,9470AG 8.5% 371.0 / 18.8 129R46R (CGT→CGCcysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base A (52/66);  new base G (5/6);  total (57/72)
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.98e-01

Only 100 of 131 total aligned reads displayed.

ACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTAT  >  Exported/1844881‑1845014
                                                                  |                                                                   
aCGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAgacgac                                                                 >  1:3736415/1‑71 (MQ=255)
  gCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGAcc                                                                >  1:141990/1‑70 (MQ=255)
       cTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCa                                                           <  1:2810687/70‑1 (MQ=255)
       cTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCa                                                           <  1:3207204/70‑1 (MQ=255)
       cTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCa                                                           <  1:4651229/70‑1 (MQ=255)
       cTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCa                                                           <  1:4645114/70‑1 (MQ=255)
       cTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCa                                                           <  1:4333932/70‑1 (MQ=255)
       cTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCa                                                           <  1:4239191/70‑1 (MQ=255)
       cTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCa                                                           <  1:1782674/70‑1 (MQ=255)
       cTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCa                                                           <  1:1657927/70‑1 (MQ=255)
       cTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCa                                                           <  1:4214054/70‑1 (MQ=255)
          tCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCagag                                                        >  1:3131617/1‑70 (MQ=255)
          tCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCagag                                                        >  1:1544888/1‑70 (MQ=255)
               aTGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACtt                                                   <  1:929970/70‑1 (MQ=255)
               aTGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACtt                                                   <  1:2103348/70‑1 (MQ=255)
               aTGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACtt                                                   <  1:3422146/70‑1 (MQ=255)
               aTGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACtt                                                   <  1:4979159/70‑1 (MQ=255)
               aTGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATAGTGCAGCAGACGACCAATCAGAGTACtt                                                   <  1:3468942/70‑1 (MQ=255)
                tGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACtttt                                                 <  1:3735925/71‑1 (MQ=255)
                tGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACtttt                                                 <  1:4224143/71‑1 (MQ=255)
                tGCAGCGATGACAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACttt                                                  <  1:4849246/70‑1 (MQ=255)
                   aGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGcc                                              >  1:1387113/1‑71 (MQ=255)
                    gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg                                             <  1:3530537/71‑1 (MQ=255)
                    gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg                                             <  1:3127250/71‑1 (MQ=255)
                    gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg                                             <  1:900538/71‑1 (MQ=255)
                    gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg                                             >  1:2317683/1‑71 (MQ=255)
                    gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg                                             >  1:2120284/1‑71 (MQ=255)
                    gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg                                             >  1:2199484/1‑71 (MQ=255)
                     cGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCa                                                           <  1:1560608/56‑1 (MQ=255)
                     cGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCa                                                           <  1:1257070/56‑1 (MQ=255)
                     cGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg                                             <  1:1092429/70‑1 (MQ=255)
                      gATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCagag                                                        >  1:3073457/1‑58 (MQ=255)
                        tGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCAt                                         >  1:4957516/1‑71 (MQ=255)
                        tGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCAt                                         >  1:1250219/1‑71 (MQ=255)
                        tGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCAt                                         >  1:599302/1‑71 (MQ=255)
                        tGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCAt                                         >  1:2819545/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:3536551/71‑1 (MQ=255)
                             gCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCATCGACg                                    <  1:1841782/71‑1 (MQ=255)
                             gCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCATCGACg                                    <  1:2805181/71‑1 (MQ=255)
                              cTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGGGTACTTTTGCCGTCATCGACg                                    <  1:4455371/70‑1 (MQ=255)
                               tGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGc                                   <  1:2104815/70‑1 (MQ=255)
                                gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                        <  1:4663480/64‑1 (MQ=255)
                                gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                        <  1:204410/64‑1 (MQ=255)
                                gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                        <  1:1168310/64‑1 (MQ=255)
                                gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                        <  1:1477745/64‑1 (MQ=255)
                                gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                        <  1:3250941/64‑1 (MQ=255)
                                 aTCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                        <  1:4515000/63‑1 (MQ=255)
                                 aTCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                        <  1:4083309/63‑1 (MQ=255)
                                 aTCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCt                                  <  1:3101550/69‑1 (MQ=255)
                                 aTCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCt                                  <  1:2609820/69‑1 (MQ=255)
                                 aTCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCt                                  <  1:4120070/69‑1 (MQ=255)
                                   cTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                    >  1:4984116/1‑65 (MQ=255)
                                   cTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                    >  1:2640055/1‑65 (MQ=255)
                                   cTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                    >  1:5151043/1‑65 (MQ=255)
                                   cTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                    >  1:5156940/1‑65 (MQ=255)
                                   cTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                    >  1:3572549/1‑65 (MQ=255)
                                    ttCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                        <  1:2564054/60‑1 (MQ=255)
                                    ttCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                        <  1:4399680/60‑1 (MQ=255)
                                      cGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                      >  1:4754769/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:5154085/41‑1 (MQ=255)
                                       gTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCagag                                                        <  1:3778932/41‑1 (MQ=255)
                                       gTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCagag                                                        <  1:913811/41‑1 (MQ=255)
                                          gATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGt                         >  1:1639488/1‑69 (MQ=255)
                                          gATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGt                         >  1:3487178/1‑69 (MQ=255)
                                            tttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAg                      <  1:1784122/70‑1 (MQ=255)
                                            tttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAg                      <  1:1656710/70‑1 (MQ=255)
                                             ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                      >  1:825490/1‑53 (MQ=255)
                                             ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                      >  1:1307302/1‑53 (MQ=255)
                                             ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                      >  1:1329625/1‑53 (MQ=255)
                                             ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                      >  1:5194300/1‑53 (MQ=255)
                                             ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                      >  1:3774012/1‑53 (MQ=255)
                                             ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGa                     >  1:4343903/1‑70 (MQ=255)
                                             ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGa                     >  1:698746/1‑70 (MQ=255)
                                             ttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGa                     >  1:3495973/1‑70 (MQ=255)
                                              tGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGaa                    <  1:3812067/70‑1 (MQ=255)
                                              tGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGaa                    <  1:1829053/70‑1 (MQ=255)
                                                 cGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCATAAAACg                <  1:11128/71‑1 (MQ=255)
                                                 cGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACg                <  1:1103728/71‑1 (MQ=255)
                                                 cGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACg                <  1:1174046/71‑1 (MQ=255)
                                                  gggTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACg                >  1:506740/1‑70 (MQ=255)
                                                  gggTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACg                >  1:3429346/1‑70 (MQ=255)
                                                   ggTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                        <  1:118449/45‑1 (MQ=255)
                                                    gTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACgcag             <  1:807625/71‑1 (MQ=255)
                                                    gTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGACGTCATCGACGCTACCACAGGTCAGAAAACgcag             <  1:3782357/71‑1 (MQ=255)
                                                      aTCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCt                                  <  1:3195889/48‑1 (MQ=255)
                                                       tCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                <  1:3597245/49‑1 (MQ=255)
                                                       tCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                <  1:2362473/49‑1 (MQ=255)
                                                       tCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                <  1:1478242/49‑1 (MQ=255)
                                                         gTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCa                       >  1:3432887/1‑56 (MQ=255)
                                                         gTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCa                       >  1:2426314/1‑56 (MQ=255)
                                                         gTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCa                       >  1:4356694/1‑56 (MQ=255)
                                                         gTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCa                       >  1:2295246/1‑56 (MQ=255)
                                                           gcagcaGGCGACCAATCAGGGTACTTTTGCCGTCATCGACg                                    <  1:3627449/41‑1 (MQ=37)
                                                                 gacgacCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTAt  <  1:4600178/69‑1 (MQ=255)
                                                                  |                                                                   
ACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTAT  >  Exported/1844881‑1845014

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: