Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,045,825 A→G 7.8% intergenic (+587/+208) yqiK → / ← rfaE conserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,045,8250AG 7.8% 213.8 / 8.9 90intergenic (+587/+208)yqiK/rfaEconserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base A (35/48);  new base G (3/4);  total (38/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 108 total aligned reads displayed.

TTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTA  >  minE/2045756‑2045895
                                                                     |                                                                       
ttGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaa                                                                        <  1:4422296/71‑1 (MQ=255)
   cgcgGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaagaa                                                                     <  1:3636028/71‑1 (MQ=255)
   cgcgGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag                                                                       >  1:1446895/1‑69 (MQ=255)
         gTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGcc                                                              >  1:4650558/1‑71 (MQ=255)
            cGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                                      >  1:4429500/1‑61 (MQ=255)
                gcTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGta                                                         >  1:5442407/1‑69 (MQ=255)
                 cTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGtata                                                       <  1:1887124/70‑1 (MQ=255)
                    atcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGtat                                                        <  1:2333700/66‑1 (MQ=255)
                    atcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAgg                                                    >  1:3343442/1‑70 (MQ=255)
                    atcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAg                                                     <  1:5466060/69‑1 (MQ=255)
                     tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTg                                                           <  1:4827154/62‑1 (MQ=255)
                     tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTg                                                           <  1:1552405/62‑1 (MQ=255)
                     tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGaa                                                  <  1:2717963/71‑1 (MQ=255)
                       aacaaTCGGGGGGCAGCAAGGGGCTGAAACTGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGa                                                   <  1:892892/68‑1 (MQ=255)
                       aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag                                                                       >  1:822613/1‑49 (MQ=255)
                       aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAg                                                >  1:5640262/1‑71 (MQ=255)
                       aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAg                                                >  1:2466718/1‑71 (MQ=255)
                        ctaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGaaa                                                 <  1:4132819/67‑1 (MQ=255)
                        acaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTg                                                           <  1:711585/59‑1 (MQ=255)
                        acaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTg                                                           <  1:1271467/59‑1 (MQ=255)
                        acaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTg                                                           <  1:5090951/59‑1 (MQ=255)
                               gggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag                                                                       <  1:3182772/41‑1 (MQ=255)
                                  ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAgg                                                    <  1:2562837/56‑1 (MQ=255)
                                  ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAgg                                                    <  1:346786/56‑1 (MQ=255)
                                  ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAgg                                                    <  1:276982/56‑1 (MQ=255)
                                   gcagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAgg                                                    <  1:165441/55‑1 (MQ=255)
                                   gcagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGa                                                   <  1:733399/56‑1 (MQ=255)
                                    cagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAgg                                                                  >  1:2605999/1‑40 (MQ=40)
                                    cagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAg                                   >  1:830053/1‑71 (MQ=255)
                                     agcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTa                                           <  1:2531398/62‑1 (MQ=255)
                                     agcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTa                                           <  1:4329956/62‑1 (MQ=255)
                                      gcaAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                                      <  1:4963063/35‑1 (MQ=40)
                                      gcaAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                                      <  1:3957118/35‑1 (MQ=40)
                                      gcaAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                                      <  1:4260276/35‑1 (MQ=40)
                                      gcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAgg                                               <  1:3000704/57‑1 (MQ=255)
                                      gcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTatg                                         >  1:4328361/1‑63 (MQ=255)
                                        aaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGt                                             >  1:781119/1‑57 (MQ=255)
                                         aGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTa                                           >  1:3736880/1‑58 (MQ=255)
                                         aGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc                                <  1:657273/69‑1 (MQ=255)
                                           gggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGTCCTTGTATAAGGAAAGGGTTa                                           <  1:1977594/56‑1 (MQ=255)
                                           acgCTGAAACGGGAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGt                              >  1:303292/3‑70 (MQ=0)
                                            ggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGtat                                                        <  1:57129/42‑1 (MQ=255)
                                            ggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                           >  1:3187602/1‑71 (MQ=255)
                                            ggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                           >  1:5627124/1‑71 (MQ=255)
                                            cgCTGAAACGGGAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGtca                            <  1:1786500/70‑1 (MQ=1)
                                              cTGAAACGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGaaa                                                 >  1:3653992/1‑47 (MQ=14)
                                              cTGAAACGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCatg                                         <  1:144998/55‑1 (MQ=14)
                                              cTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATaa                                                      >  1:442617/1‑42 (MQ=255)
                                                gAAACGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGaa                                                  >  1:1492151/1‑44 (MQ=14)
                                                gAAACGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACtcat                           >  1:1889915/1‑67 (MQ=14)
                                                gAAACGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACtcat                           >  1:1896669/1‑67 (MQ=14)
                                                gAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGtt                                            <  1:1190765/50‑1 (MQ=255)
                                                 aaaCGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCatgat                                       <  1:4403377/54‑1 (MQ=14)
                                                 aaaCGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCatgat                                       <  1:3787364/54‑1 (MQ=14)
                                                 aaaCGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTcatcatcg                      <  1:565443/71‑3 (MQ=11)
                                                 aaaCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATAc                      >  1:4899675/1‑71 (MQ=255)
                                                  aaCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAg                                   >  1:3535738/1‑57 (MQ=255)
                                                   aCGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCAtgatga                                     >  1:3633700/1‑54 (MQ=14)
                                                   aCGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCAtgatga                                     >  1:3590871/1‑54 (MQ=14)
                                                   aCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAATGAAAGGGTTATGATGAAGCTCGt                              <  1:5193989/61‑1 (MQ=255)
                                                      ggAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCatgat                                       <  1:2622584/49‑1 (MQ=14)
                                                      ggAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtg                 >  1:1155474/1‑70 (MQ=2)
                                                      ggAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgtt                 >  1:3793091/1‑71 (MQ=255)
                                                       gAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGTCATCATACTGgt                  >  1:3680797/1‑70 (MQ=11)
                                                       gAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtgc                >  1:3332284/1‑69 (MQ=2)
                                                       gAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttg                <  1:2081641/71‑1 (MQ=255)
                                                        aaaGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtgc                >  1:4164927/1‑68 (MQ=2)
                                                        aaaGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtgc                >  1:3132334/1‑68 (MQ=2)
                                                          aGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGt                                             <  1:22867/40‑1 (MQ=0)
                                                             cccTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACt                     <  1:1253870/60‑1 (MQ=255)
                                                              tcTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGTCATCATACTGgtg                 <  1:3123410/63‑2 (MQ=11)
                                                              ccTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttgt          >  1:3907027/1‑70 (MQ=255)
                                                               cTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtc                             <  1:3373786/50‑1 (MQ=255)
                                                               cTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTa        <  1:4919004/71‑1 (MQ=255)
                                                                tCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCat                                          <  1:128161/36‑1 (MQ=12)
                                                                tCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCat                                          <  1:2965762/36‑1 (MQ=12)
                                                                tCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCat                                          <  1:1943990/36‑1 (MQ=12)
                                                                 cccGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCtac                               >  1:74665/1‑44 (MQ=9)
                                                                 cccGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGc                                  >  1:3127908/1‑43 (MQ=255)
                                                                 cccGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                       >  1:4898656/1‑54 (MQ=255)
                                                                  ccGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGTCATCATACTGgtga                >  1:3161910/1‑59 (MQ=2)
                                                                  ccGAGGAAGGGGCCATAAATAAGGAAAGGGTCatga                                        >  1:4550644/1‑36 (MQ=12)
                                                                  ccGAGGAAGGGGCCATAAATAAGGAAAGGGTCatga                                        >  1:1660542/1‑36 (MQ=12)
                                                                  ccGAGGAAGGGGCCATAAATAAGGAAAGGGTCatga                                        >  1:5472944/1‑36 (MQ=12)
                                                                  ccGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTcatcat                        <  1:3479904/52‑1 (MQ=9)
                                                                  ccGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTcatcat                        <  1:2184869/52‑1 (MQ=9)
                                                                  ccGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTcatcat                        <  1:213885/52‑1 (MQ=9)
                                                                  ccGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttgt          <  1:2337461/66‑1 (MQ=255)
                                                                   cGAGGAAGGGGCCATAAATAAGGAAAGGGTCatgat                                       <  1:2218638/36‑1 (MQ=9)
                                                                   cGAGGAAGGGGCCATAAATAAGGAAAGGGTCatgat                                       <  1:3598393/36‑1 (MQ=9)
                                                                   cGAGGAAGGGGCCATAAATAAGGAAAGGGTCatgat                                       <  1:3931266/36‑1 (MQ=9)
                                                                   cGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                       >  1:1020149/1‑52 (MQ=255)
                                                                    gaagaaGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                       >  1:525181/1‑51 (MQ=255)
                                                                    gaagaaGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt   <  1:4307788/71‑1 (MQ=255)
                                                                    gaagaaGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt   <  1:4698006/71‑1 (MQ=255)
                                                                    gaagaaGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGt    >  1:5507266/1‑70 (MQ=255)
                                                                     aagaagGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgt                  >  1:1201610/1‑55 (MQ=255)
                                                                     aagaagGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt   >  1:3440667/1‑70 (MQ=255)
                                                                     aagaagGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTa  <  1:1849759/71‑1 (MQ=255)
                                                                     aagaagGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTa  <  1:2702621/71‑1 (MQ=255)
                                                                     |                                                                       
TTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTA  >  minE/2045756‑2045895

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: