Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,591,605 G→A 26.4% intergenic (+64/+125) yfeS → / ← cysM conserved hypothetical protein/cysteine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,591,6050GA26.4% 135.4 / 55.5 87intergenic (+64/+125)yfeS/cysMconserved hypothetical protein/cysteine synthase B
Reads supporting (aligned to +/- strand):  ref base G (37/27);  new base A (11/12);  total (48/39)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

ATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGAT  >  minE/1591537‑1591672
                                                                    |                                                                   
aTTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                                                   <  1:1853068/71‑1 (MQ=255)
aTTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                                                   <  1:3758380/71‑1 (MQ=255)
 ttAAATACTGCATTTGTCGTCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaa                                                                  <  1:3191428/71‑1 (MQ=255)
 ttAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaa                                                                  <  1:941153/71‑1 (MQ=255)
  tAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                 >  1:386799/1‑71 (MQ=255)
      tACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGcc                                                             <  1:2763832/71‑1 (MQ=255)
      tACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGcc                                                             <  1:3947953/71‑1 (MQ=255)
           gcgtaggcctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTa                                                          >  1:2336942‑M1/48‑69 (MQ=255)
           gcgtaggcctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATc                                                        <  1:3770279‑M1/24‑1 (MQ=255)
           gcgtaggcctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATc                                                        <  1:3496256‑M1/24‑1 (MQ=255)
           gcgtaggcctgataagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCtt                                                           >  1:2146674‑M1/48‑68 (MQ=255)
           gcgtaggcctgataagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTa                                                          >  1:2597252‑M1/48‑69 (MQ=255)
           cATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGcc                                                             >  1:1538075/1‑66 (MQ=255)
            gataggcctgataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTAt                                                         >  1:2039726‑M1/47‑69 (MQ=255)
            cgtaggcctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATcc                                                       <  1:2167983‑M1/25‑1 (MQ=255)
            cgtaggcctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATcc                                                       >  1:1834397‑M1/47‑71 (MQ=255)
             gtaggcctgataagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCg                                                      >  1:2838742‑M1/46‑71 (MQ=255)
              ttGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCgg                                                     >  1:3899039/1‑71 (MQ=255)
              taggcctgataagacgctgcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCcga                                                     >  1:1896200‑M1/45‑70 (MQ=255)
              taggcctgataagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCcga                                                     >  1:1208360‑M1/45‑70 (MQ=255)
                 tCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTTTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                  <  1:1372906/71‑1 (MQ=255)
                 gcctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCt                                                  >  1:1841386‑M1/42‑71 (MQ=255)
                  cctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTa                                                 <  1:2784971‑M1/31‑1 (MQ=255)
                   ctgataagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTa                                                 <  1:2708781‑M1/31‑1 (MQ=255)
                    tgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTAc                                                <  1:1235923‑M1/32‑1 (MQ=255)
                    tgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTAc                                                <  1:1576870‑M1/32‑1 (MQ=255)
                    tgataagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTAc                                                <  1:2287383‑M1/32‑1 (MQ=255)
                    tgataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTaaa                                                                 <  1:617188‑M1/15‑1 (MQ=255)
                    gcagcaACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                    >  1:2064764/1‑66 (MQ=255)
                     gataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACa                                               <  1:3765823‑M1/33‑1 (MQ=255)
                      ataagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCg                                                      >  1:287855‑M1/37‑62 (MQ=255)
                      ataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACaa                                              <  1:2019613‑M1/34‑1 (MQ=255)
                       gcaACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                            <  1:3877167/71‑1 (MQ=255)
                        caacaaCTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                                                   <  1:2133516/47‑1 (MQ=255)
                         agacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAAt                                           >  1:1342290‑M1/34‑70 (MQ=255)
                          gacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATc                                          >  1:2433217‑M1/33‑70 (MQ=255)
                          gacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATc                                          >  1:1099010‑M1/33‑70 (MQ=255)
                          acaacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCt                                                            <  1:815401/52‑1 (MQ=255)
                           caacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGt                                        <  1:3351350/71‑1 (MQ=255)
                           acgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAatcac                                        >  1:1139882‑M1/32‑69 (MQ=255)
                            cgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATaaa                                                                 <  1:243830‑M1/15‑1 (MQ=255)
                            cgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTaaa                                                                 <  1:1168714‑M1/15‑1 (MQ=255)
                            aacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                       >  1:3868984/1‑71 (MQ=255)
                            aacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                       <  1:3597778/71‑1 (MQ=255)
                            aacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                       >  1:2877013/1‑71 (MQ=255)
                             gcgttaagcgtcgcatccggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAaatcacg                                       <  1:3529163‑M1/41‑4 (MQ=255)
                              cTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                      >  1:1170133/1‑70 (MQ=255)
                                gTTAAAAAAGTGCGCTTTGTTTATGCCGTATgcggcg                                                                     <  1:2186433/37‑1 (MQ=255)
                                 ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATc                                                        >  1:2143882/1‑49 (MQ=255)
                                 ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTaa                                   <  1:2719/70‑1 (MQ=255)
                                 taagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                  <  1:873801‑M1/46‑2 (MQ=255)
                                  aagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTa                                                          >  1:1947796‑M1/25‑46 (MQ=255)
                                  aagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat                                 >  1:2248164‑M1/25‑69 (MQ=255)
                                  aagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat                                 >  1:2856372‑M1/25‑69 (MQ=255)
                                  aagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat                                 <  1:3521988‑M1/47‑3 (MQ=255)
                                  aagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat                                 >  1:1298180‑M1/25‑69 (MQ=255)
                                  aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaat                                 <  1:3225617‑M1/47‑3 (MQ=255)
                                   agcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACgctaaatc                                >  1:351329‑M1/24‑68 (MQ=255)
                                   agcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaat                                 >  1:541617‑M1/24‑68 (MQ=255)
                                   agcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACgctaaatc                                >  1:800747‑M1/24‑68 (MQ=255)
                                    gcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATcc                                                       >  1:3308408‑M1/23‑47 (MQ=255)
                                     cgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACaa                                              >  1:969151‑M1/22‑55 (MQ=255)
                                     aaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCaa                              >  1:312931/1‑71 (MQ=35)
                                      gtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTAt                                                         <  1:1628412‑M1/23‑1 (MQ=255)
                                      gtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAAtcagg                             >  1:279019‑M1/21‑69 (MQ=255)
                                       tcgcatcaggcgttgatgtCGGATGCGGCATTAACGCCTTATCCGACCTACAAAATCACGCTAAatcaggc                            >  1:1712698‑M1/20‑68 (MQ=255)
                                          tGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc                                <  1:120748/64‑1 (MQ=35)
                                          catcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATCa                               >  1:2714985‑M1/17‑65 (MQ=255)
                                           atcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                  >  1:189275‑M1/16‑60 (MQ=255)
                                             gcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc                      <  1:2956338/71‑1 (MQ=39)
                                             caggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                              <  1:1515700‑M1/50‑2 (MQ=255)
                                             caggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                              <  1:573430‑M1/50‑2 (MQ=255)
                                             caggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                              <  1:3023827‑M1/50‑2 (MQ=255)
                                                ttGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCaa                              <  1:1793520/60‑1 (MQ=35)
                                                 tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa                  >  1:178762/1‑71 (MQ=11)
                                                  gttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACa                                               >  1:2163208‑M1/9‑41 (MQ=255)
                                                  gTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGaaa                 >  1:3629449/1‑71 (MQ=255)
                                                   ttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat                                 >  1:421014‑M1/8‑52 (MQ=255)
                                                    ttATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAcc               >  1:2415884/1‑71 (MQ=255)
                                                     tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAcc               <  1:2674011/70‑1 (MQ=255)
                                                      atgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATCa                               >  1:2038921‑M1/5‑53 (MQ=255)
                                                       tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTg            <  1:3094598/71‑1 (MQ=255)
                                                        gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTg            >  1:368045/1‑70 (MQ=255)
                                                            gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaatt                >  1:7191/1‑60 (MQ=11)
                                                            gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATAGTGCAAATTCAATACATTGCAGaatt                >  1:521196/1‑60 (MQ=11)
                                                            gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAgg        >  1:914177/1‑70 (MQ=255)
                                                              tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa                  >  1:2794054/1‑58 (MQ=11)
                                                              tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTa          >  1:3191892/1‑66 (MQ=255)
                                                                cggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGa   >  1:1079344/1‑71 (MQ=255)
                                                                 ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGAt  <  1:3787615/71‑1 (MQ=255)
                                                                 ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGAt  <  1:1043999/71‑1 (MQ=255)
                                                                 ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGATAATTCAATATATTGCAGAAACCTTGTAGGCCTGa   >  1:1335392/1‑70 (MQ=255)
                                                                   cgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCtt             <  1:902035/58‑1 (MQ=255)
                                                                    |                                                                   
ATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGAT  >  minE/1591537‑1591672

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: