Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,685,303 T→C 18.2% Y42H (TAT→CAT)  sseA → 3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,685,3030TC18.2% 97.7 / 16.5 44Y42H (TAT→CAT) sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (11/25);  new base C (5/3);  total (16/28)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCC  >  Exported/1685233‑1685370
                                                                      |                                                                   
cccGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGt                                                                     <  1:1886231/71‑1 (MQ=255)
 ccGGAAATTCAGATTATCGATGCCGGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTa                                                                    <  1:3131137/71‑1 (MQ=255)
  cGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGCAt                                                                   <  1:1962992/71‑1 (MQ=255)
   ggAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATc                                                                  >  1:2150502/1‑71 (MQ=255)
          cAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATg                                                            <  1:261106/70‑1 (MQ=255)
           aGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATgg                                                           <  1:529065/70‑1 (MQ=255)
           aGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATgg                                                           <  1:2809994/70‑1 (MQ=255)
                 tCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACatat                                                     <  1:2548652/70‑1 (MQ=255)
                       cccGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGgc                                              >  1:1338797/1‑71 (MQ=255)
                        ccGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTAt                                                                   <  1:1158583/49‑1 (MQ=255)
                        ccGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTAt                                                                   <  1:459778/49‑1 (MQ=255)
                        ccGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTAt                                                                   <  1:247374/49‑1 (MQ=255)
                        ccGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTAt                                                                   <  1:2402950/49‑1 (MQ=255)
                        ccGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTAt                                                                   <  1:2083697/49‑1 (MQ=255)
                           cATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCt                                                                 <  1:3140603/48‑1 (MQ=255)
                                cGTCGCCTGGACAGGAGGAGCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGttt                                     <  1:1708135/71‑1 (MQ=255)
                                  tCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGttttt                                   <  1:913630/71‑1 (MQ=255)
                                    gCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGtttttt                                  <  1:2201177/70‑1 (MQ=255)
                                     ccTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACata                                                      >  1:2257445/1‑49 (MQ=255)
                                     ccTGGACAGGAGGATCGTAACGTTGCTCAGGAGCATCTGAATGGACATATTcc                                                  <  1:1874599/53‑1 (MQ=255)
                                        ggACAGGAGGATCGTAACGTTGCTCAGGAGCATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGatat                             >  1:263942/1‑71 (MQ=255)
                                          aCAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATc                            >  1:1242493/1‑70 (MQ=255)
                                          aCAGGAGGATCGTAACGTTGCTCAGGAGCATCTGAATGGACat                                                       <  1:3139981/43‑1 (MQ=255)
                                            aggaggATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCg                           >  1:2570978/1‑69 (MQ=255)
                                            aggaggATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCg                           >  1:588821/1‑69 (MQ=255)
                                            aggaggATCGTAACGTTGCTCAGGAGCATCTGAATGGACatat                                                     >  1:1127061/1‑43 (MQ=255)
                                            aggaggATCGTAACGTTGCTCAGGAGCATCTGAATGGACatat                                                     >  1:2501366/1‑43 (MQ=255)
                                              gaggaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAgc                       <  1:1803145/71‑1 (MQ=255)
                                              gaggaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAgc                       <  1:1798885/71‑1 (MQ=255)
                                              gaggaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAgc                       >  1:2529086/1‑71 (MQ=255)
                                                 gaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTc                                                   <  1:965393/40‑1 (MQ=255)
                                                 gaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTc                                                   <  1:2389655/40‑1 (MQ=255)
                                                 gaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTc                                                   <  1:3705291/40‑1 (MQ=255)
                                                 gaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCt                    <  1:3237934/71‑1 (MQ=255)
                                                 gaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCt                    <  1:775032/71‑1 (MQ=255)
                                                     gTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCt                >  1:2997785/1‑71 (MQ=255)
                                                        aCGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGAt             <  1:2518296/71‑1 (MQ=255)
                                                        aCGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGAt             >  1:577740/1‑71 (MQ=255)
                                                           ttGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCg                           <  1:2829671/54‑1 (MQ=255)
                                                              cTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCacac        <  1:3389062/70‑1 (MQ=255)
                                                              cTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCacac        <  1:3497530/70‑1 (MQ=255)
                                                                 aGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAg                        >  1:527517/1‑51 (MQ=255)
                                                                 aGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTcc    >  1:1205152/1‑71 (MQ=255)
                                                                   gAGCATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTcccc  >  1:3585313/1‑71 (MQ=255)
                                                                   gAGCATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTcccc  >  1:2827037/1‑71 (MQ=255)
                                                                      tATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATc                            <  1:3158695/42‑1 (MQ=255)
                                                                      tATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCtt               <  1:567408/55‑2 (MQ=255)
                                                                      |                                                                   
CCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCC  >  Exported/1685233‑1685370

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: