Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,591,637 T→A 18.6% intergenic (+96/+93) 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,6370TA18.6% 227.5 / 36.7 91intergenic (+96/+93)yfeS/cysMconserved hypothetical protein/cysteine synthase B
Reads supporting (aligned to +/- strand):  ref base T (43/31);  new base A (7/10);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

GTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCA  >  Exported/1591569‑1591706
                                                                    |                                                                     
ttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTaa                                                                     >  1:2207788‑M1/27‑71 (MQ=255)
ttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTaa                                                                     >  1:4081563‑M1/27‑71 (MQ=255)
gTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTaa                                                                     >  1:3292707/1‑71 (MQ=255)
 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:2920001‑M1/25‑69 (MQ=255)
  aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                                                    >  1:2446982‑M1/25‑69 (MQ=255)
  aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                                                    >  1:3968611‑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: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)
                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:228643‑M1/50‑2 (MQ=255)
                     gatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                                                <  1:4317859‑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)
                             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)
                                    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)
                                     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:2702139/1‑47 (MQ=11)
                                       aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                      >  1:2013041/1‑47 (MQ=11)
                                        aCGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTa                             <  1:2269605/71‑1 (MQ=11)
                                          gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaatt                                                  <  1:4338094/48‑3 (MQ=11)
                                            attATCCGGCCTACAAAATCGTGCTAATTCAATATATTGGAGAAACCTTg                                              <  1:4279856/49‑1 (MQ=255)
                                              tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTa                                            >  1:3572312/1‑50 (MQ=255)
                                              tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                       >  1:1372159/1‑71 (MQ=255)
                                              tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                       >  1:1676644/1‑71 (MQ=255)
                                              tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                       >  1:1725128/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)
                                                                 tGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGt    <  1:5566025/71‑1 (MQ=255)
                                                                 tGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGt    <  1:61967/71‑1 (MQ=255)
                                                                   cTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtca  <  1:2307079/71‑1 (MQ=255)
                                                                   cTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtca  <  1:1215325/71‑1 (MQ=255)
                                                                   cTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtc   <  1:1925143/70‑1 (MQ=255)
                                                                    tAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGttt             >  1:3941392/1‑59 (MQ=255)
                                                                    tAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGttt             >  1:3445288/1‑59 (MQ=255)
                                                                    |                                                                     
GTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCA  >  Exported/1591569‑1591706

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: