Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,591,630 A→T 10.8% intergenic (+89/+100) yfeS → / ← cysM conserved hypothetical protein/cysteine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,591,6300AT10.8% 293.8 / 24.2 111intergenic (+89/+100)yfeS/cysMconserved hypothetical protein/cysteine synthase B
Reads supporting (aligned to +/- strand):  ref base A (58/41);  new base T (4/8);  total (62/49)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 120 total aligned reads displayed.

GCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCG  >  Exported/1591560‑1591699
                                                                      |                                                                     
taagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACaaaa                                                                       >  1:2930468‑M1/36‑71 (MQ=255)
taagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACaaaa                                                                       <  1:673634‑M1/36‑1 (MQ=255)
 aagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACaaaa                                                                       <  1:3316933‑M1/36‑1 (MQ=255)
 aagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAAt                                                                      <  1:3698561‑M1/37‑1 (MQ=255)
 aagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAAt                                                                      <  1:5355782‑M1/37‑1 (MQ=255)
 aagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAAt                                                                      <  1:701828‑M1/37‑1 (MQ=255)
 aagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAAt                                                                      <  1:4347264‑M1/37‑1 (MQ=255)
  agacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATc                                                                     >  1:3931994‑M1/34‑71 (MQ=255)
  aacaacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATc                                                                     >  1:4471248/1‑71 (MQ=255)
  aacaacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATc                                                                     >  1:1993540/1‑71 (MQ=255)
   gacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAAtca                                                                    >  1:1472609‑M1/33‑70 (MQ=255)
   gacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATc                                                                     >  1:4360767‑M1/33‑70 (MQ=255)
   gacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAAtca                                                                    <  1:2232205‑M1/39‑2 (MQ=255)
         ttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACaaaa                                                                       >  1:670278‑M1/27‑62 (MQ=255)
         ttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTaa                                                              >  1:2207788‑M1/27‑71 (MQ=255)
         gTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTaa                                                              >  1:3292707/1‑71 (MQ=255)
          ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAt                                                                      <  1:2143460/62‑1 (MQ=39)
          ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                 >  1:235094/1‑67 (MQ=39)
          ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTaa                                                              <  1:2770226/70‑1 (MQ=255)
          ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAAt                                                             >  1:4599416/1‑71 (MQ=255)
          caagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                                             <  1:4873076‑M1/46‑2 (MQ=255)
           tAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAAtt                                                            >  1:2591310/1‑71 (MQ=40)
           tAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAAtt                                                            >  1:2474244/1‑71 (MQ=40)
           aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                                             >  1:3968611‑M1/25‑69 (MQ=255)
           aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                                             >  1:2920001‑M1/25‑69 (MQ=255)
            agcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                                             <  1:2228188‑M1/46‑2 (MQ=255)
            agcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACgctaaatc                                                           >  1:4238166‑M1/24‑68 (MQ=255)
            agcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaat                                                            <  1:2089850‑M1/47‑3 (MQ=255)
            agcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaat                                                            <  1:2082957‑M1/47‑3 (MQ=255)
             gcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATCa                                                          >  1:4811518‑M1/23‑71 (MQ=255)
             aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc                                                           <  1:1177903/70‑1 (MQ=39)
             aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCa                                                          <  1:2427778/71‑1 (MQ=38)
             aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCa                                                          <  1:3609408/71‑1 (MQ=38)
              cgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                                         >  1:5546817‑M1/22‑70 (MQ=255)
              cgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                                         >  1:187267‑M1/22‑70 (MQ=255)
              cgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                                         <  1:519117‑M1/50‑2 (MQ=255)
               gtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAAtcagg                                                        <  1:1299606‑M1/51‑3 (MQ=255)
                tcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat                                                            >  1:3391076‑M1/20‑64 (MQ=255)
                tcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat                                                            >  1:3384337‑M1/20‑64 (MQ=255)
                tcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAatcaggc                                                       >  1:3504489‑M1/20‑68 (MQ=255)
                tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACGTACAAAATCACGCTAAatcaggc                                                       <  1:3454996‑M1/52‑4 (MQ=255)
                tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                                             >  1:462860‑M1/20‑64 (MQ=255)
                tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAatcaggc                                                       >  1:5118810‑M1/20‑68 (MQ=255)
                tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAatcaggc                                                       >  1:4696563‑M1/20‑68 (MQ=255)
                tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAatcagga                                                       >  1:1080715‑M1/20‑68 (MQ=255)
                 cgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaat                                                            >  1:754138‑M1/19‑63 (MQ=255)
                  gTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAatata                                                     <  1:3988206/71‑1 (MQ=39)
                   tGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAatatat                                                    <  1:2293409/71‑1 (MQ=39)
                       cTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc                                                 <  1:4304830/70‑1 (MQ=40)
                         gcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                                         >  1:3826488‑M1/11‑59 (MQ=255)
                           gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaat                                            >  1:2029427/1‑70 (MQ=11)
                           gTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc                                                           <  1:2608126/56‑1 (MQ=35)
                              tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAc                                           >  1:2781268/1‑69 (MQ=255)
                              tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCt                                         >  1:2820729/1‑71 (MQ=255)
                              tATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaatca                                         >  1:305977/1‑68 (MQ=11)
                              tATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaatc                                          <  1:4377490/70‑3 (MQ=11)
                              gatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                                         <  1:228643‑M1/50‑2 (MQ=255)
                                     gATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGc                                  >  1:2859712/1‑71 (MQ=11)
                                     gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCt                                                                <  1:321210/41‑1 (MQ=32)
                                     gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAgg                                   >  1:161687/1‑70 (MQ=255)
                                     gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGc                                  <  1:487477/71‑1 (MQ=255)
                                     gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGc                                  >  1:3033291/1‑71 (MQ=255)
                                      aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc                                                 >  1:1538705/1‑55 (MQ=255)
                                       tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCt                                <  1:4258487/71‑1 (MQ=255)
                                        gcggcgTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTa                                     >  1:4840718/1‑65 (MQ=11)
                                         cggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAc                                           >  1:5689638/1‑58 (MQ=255)
                                            cttAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATaa                           <  1:1974287/69‑1 (MQ=255)
                                             gTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg                          >  1:1172367/1‑71 (MQ=255)
                                             gTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg                          <  1:291041/71‑1 (MQ=255)
                                              tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaat                                           <  1:5390949/53‑2 (MQ=11)
                                              tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaat                                           <  1:3676233/53‑2 (MQ=11)
                                              tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAg                                    <  1:1589502/60‑1 (MQ=11)
                                              tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCt                                <  1:3975788/64‑1 (MQ=11)
                                              tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCa                                                          <  1:750435/38‑1 (MQ=40)
                                              tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg                          <  1:579092/70‑1 (MQ=255)
                                                aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCaa                                                         >  1:3982589/1‑37 (MQ=40)
                                                aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCaa                                                         >  1:2654557/1‑37 (MQ=40)
                                                aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTg                                                  >  1:2682872/1‑44 (MQ=38)
                                                aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg                        >  1:3221703/1‑70 (MQ=255)
                                                aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg                        >  1:3680747/1‑70 (MQ=255)
                                                aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                               >  1:2013041/1‑47 (MQ=11)
                                                   gCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTAgc                    >  1:2944531/1‑71 (MQ=11)
                                                   gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaatt                                           <  1:4338094/48‑3 (MQ=11)
                                                     cTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgcg                   >  1:234498/1‑70 (MQ=255)
                                                     attATCCGGCCTACAAAATCGTGCTAATTCAATATATTGGAGAAACCTTg                                       <  1:4279856/49‑1 (MQ=255)
                                                       tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTa                                     >  1:3572312/1‑50 (MQ=255)
                                                       tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                >  1:1725128/1‑71 (MQ=255)
                                                       tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                >  1:1676644/1‑71 (MQ=255)
                                                        ttCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                             <  1:329634/40‑1 (MQ=11)
                                                        aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCtt                                        <  1:136632/46‑1 (MQ=255)
                                                        aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                >  1:596644/1‑70 (MQ=255)
                                                        aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                >  1:1384668/1‑70 (MQ=255)
                                                         tCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTAGCGCATCa              <  1:4161649/71‑1 (MQ=11)
                                                                tACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTg                               <  1:559899/47‑1 (MQ=255)
                                                                tACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgcgc                  <  1:5312983/60‑1 (MQ=255)
                                                                  cAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATaa                           >  1:5266805/1‑49 (MQ=255)
                                                                  cAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCa                 >  1:2144260/1‑59 (MQ=255)
                                                                    aaaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtt       >  1:4591587/1‑67 (MQ=255)
                                                                     aaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg                        >  1:5019737/1‑49 (MQ=255)
                                                                      aTCGTGCTAATTCAATATATTGCAGAAACCTTGTa                                     >  1:2087793/1‑35 (MQ=255)
                                                                      |                                                                     
GCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCG  >  Exported/1591560‑1591699

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: