Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,591,633 G→A 31.9% intergenic (+92/+97) 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,6330GA31.9% 80.6 / 75.7 101intergenic (+92/+97)yfeS/cysMconserved hypothetical protein/cysteine synthase B
Reads supporting (aligned to +/- strand):  ref base G (35/23);  new base A (21/10);  total (60/41)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCG  >  Exported/1591569‑1591699
                                                                |                                                                  
ttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTaa                                                              >  1:4081563‑M1/27‑71 (MQ=255)
ttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTaa                                                              >  1:2207788‑M1/27‑71 (MQ=255)
gTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTaa                                                              >  1:3292707/1‑71 (MQ=255)
 ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                 >  1:235094/1‑67 (MQ=39)
 ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTaa                                                              <  1:2770226/70‑1 (MQ=255)
 ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAAt                                                             >  1:4599416/1‑71 (MQ=255)
 caagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                                             >  1:209120‑M1/26‑70 (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)
  aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                                             >  1:2446982‑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:655331‑M1/23‑71 (MQ=255)
    gcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATCa                                                          >  1:4811518‑M1/23‑71 (MQ=255)
    aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc                                                           <  1:1177903/70‑1 (MQ=39)
    aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCa                                                          <  1:3609408/71‑1 (MQ=38)
    aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCa                                                          <  1:2427778/71‑1 (MQ=38)
     cgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                                         >  1:187267‑M1/22‑70 (MQ=255)
     cgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                                         >  1:5546817‑M1/22‑70 (MQ=255)
     cgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                                         >  1:3585021‑M1/22‑70 (MQ=255)
     cgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                                         <  1:519117‑M1/50‑2 (MQ=255)
      gtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAAtcagg                                                        <  1:1299606‑M1/51‑3 (MQ=255)
       tcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat                                                            >  1:3384337‑M1/20‑64 (MQ=255)
       tcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat                                                            >  1:3391076‑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)
          catcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATCa                                                          <  1:3876066‑M1/49‑1 (MQ=255)
           gcgcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATAtt                                                   >  1:1995596/1‑71 (MQ=39)
              cTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc                                                 <  1:4304830/70‑1 (MQ=40)
                gcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                                         >  1:3826488‑M1/11‑59 (MQ=255)
                  gTTTATGCCGGATTCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc                                                           <  1:3353321/56‑1 (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:4317859‑M1/50‑2 (MQ=255)
                     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:3033291/1‑71 (MQ=255)
                            gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGc                                  <  1:487477/71‑1 (MQ=255)
                             aTGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaa                                             <  1:4827526/59‑1 (MQ=11)
                             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:3676233/53‑2 (MQ=11)
                                     tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaat                                           <  1:5390949/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)
                                     tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAACCCTTGt                                      <  1:3206/58‑1 (MQ=255)
                                     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)
                                       aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                               >  1:2702139/1‑47 (MQ=11)
                                        aCGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTa                      <  1:2269605/71‑1 (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:1372159/1‑71 (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:1384668/1‑70 (MQ=255)
                                               aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                >  1:596644/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)
                                                            aaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCg  >  1:5153723/1‑71 (MQ=255)
                                                             aTCGTGCTAATTCAATATATTGCAGAAACCTTGTa                                     >  1:2087793/1‑35 (MQ=255)
                                                                |                                                                  
GTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCG  >  Exported/1591569‑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: