Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,045,826 A→G 8.1% intergenic (+588/+207) 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,8260AG 8.1% 163.7 / 6.5 62intergenic (+588/+207)yqiK/rfaEconserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base G (3/2);  total (31/31)
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.03e-01

GCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTT  >  minE/2045758‑2045894
                                                                    |                                                                     
gcgcGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag                                                                      <  1:913967/70‑1 (MQ=255)
          cGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTg                                                          <  1:1179650/71‑1 (MQ=255)
             cgcTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtat                                                       <  1:834630/71‑1 (MQ=255)
                   tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                                     <  1:3946819/52‑1 (MQ=255)
                   tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGaa                                                 <  1:1718862/71‑1 (MQ=255)
                   tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGaa                                                 <  1:1971785/71‑1 (MQ=255)
                    caacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                                     <  1:3901906/51‑1 (MQ=255)
                     aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAg                                               >  1:317614/1‑71 (MQ=255)
                     aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAg                                               >  1:2129772/1‑71 (MQ=255)
                         aTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGtt                                           <  1:4136953/71‑1 (MQ=255)
                         aTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGt                                            <  1:3753043/70‑1 (MQ=255)
                          tCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTa                                          >  1:1746001/1‑71 (MQ=255)
                          tCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTa                                          >  1:4079186/1‑71 (MQ=255)
                           cGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat                                         >  1:3696942/1‑71 (MQ=255)
                            ggggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat                                         <  1:722908/70‑1 (MQ=255)
                            ggggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat                                         >  1:1486467/1‑70 (MQ=255)
                            ggggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat                                         >  1:3568808/1‑70 (MQ=255)
                               gggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTatgatg                                     <  1:3509603/71‑1 (MQ=255)
                                ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTAtgatga                                    <  1:3685976/71‑1 (MQ=255)
                                   agcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtat                                                       <  1:2458235/49‑1 (MQ=255)
                                   agcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATaa                                                     >  1:1678597/1‑51 (MQ=255)
                                   agcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATaa                                                     >  1:50707/1‑51 (MQ=255)
                                    gcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtat                                                       <  1:2291861/48‑1 (MQ=255)
                                     caAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc                               <  1:975157/71‑1 (MQ=255)
                                      aaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCg                              >  1:1111609/1‑71 (MQ=255)
                                       aGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAgg                                              <  1:2603469/54‑1 (MQ=255)
                                       aGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGt                             >  1:2384862/1‑71 (MQ=255)
                                         gggCTGAATCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTAGTATAAGGAAAGGGt                                            <  1:3714754/54‑1 (MQ=18)
                                         gggCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtca                           >  1:3981406/1‑71 (MQ=255)
                                         acgCTGAAACGGGAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGtc                            <  1:1423385/69‑1 (MQ=0)
                                          ggCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                          >  1:1375826/1‑71 (MQ=255)
                                          ggCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                          >  1:3860975/1‑71 (MQ=255)
                                           gCTGAAACGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAgg                                                   <  1:574752/45‑1 (MQ=11)
                                           gCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTatga                                       >  1:4056918/1‑57 (MQ=255)
                                            cTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcatc                         >  1:2487077/1‑70 (MQ=255)
                                            cTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTCATGATGAAGCTCGtcatca                        <  1:3411888/71‑1 (MQ=255)
                                              gAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGta                                                        >  1:307769/1‑37 (MQ=255)
                                              gAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGt                             <  1:1798331/64‑1 (MQ=255)
                                              gAAACGGGAAAGCCCCTACCGAAGAAG‑GGGCCTTGTATAAgta                                                  >  1:77338/1‑41 (MQ=18)
                                               aaaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGa                                                  <  1:4084898/42‑1 (MQ=255)
                                               aaaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTatg                                        <  1:1446976/52‑1 (MQ=255)
                                               aaaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTatg                                        <  1:4055808/52‑1 (MQ=255)
                                                aaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACt                    <  1:83768/71‑1 (MQ=255)
                                                 aCGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGaaa                                                >  1:2312814/1‑42 (MQ=14)
                                                    ggAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtg                >  1:911323/1‑70 (MQ=2)
                                                    ggAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgtt                >  1:552725/1‑71 (MQ=255)
                                                      aaaGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGTCATCATACTGgtga               <  1:3016519/71‑3 (MQ=2)
                                                       aaGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg             >  1:183808/1‑71 (MQ=255)
                                                       aaGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg             >  1:3871001/1‑71 (MQ=255)
                                                       aaGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgt              >  1:3866165/1‑70 (MQ=255)
                                                        aGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGt                                            >  1:1868971/1‑39 (MQ=14)
                                                         gCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc                               >  1:3552264/1‑51 (MQ=255)
                                                         gCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttg               >  1:2698134/1‑67 (MQ=255)
                                                           cccTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttgt         >  1:3495394/1‑71 (MQ=255)
                                                           cccTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttgt         <  1:124610/71‑1 (MQ=255)
                                                            ccTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTgttgtt        <  1:355828/71‑1 (MQ=255)
                                                             cTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGtc                                           <  1:2982174/35‑2 (MQ=19)
                                                             cTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcatc                         <  1:957625/53‑1 (MQ=255)
                                                              tCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat                                         <  1:2025853/36‑1 (MQ=255)
                                                               cccGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg             >  1:2567221/1‑63 (MQ=255)
                                                               cccGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACt     >  1:3510575/1‑71 (MQ=255)
                                                                ccGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGaa                                   >  1:2553070/1‑40 (MQ=255)
                                                                ccGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                          <  1:711453/49‑1 (MQ=255)
                                                                ccGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACt     <  1:238083/70‑1 (MQ=255)
                                                                 cGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcatca                        <  1:2044391/50‑1 (MQ=255)
                                                                  gaagaaG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt  <  1:2295558/71‑1 (MQ=255)
                                                                   aagaag‑GGGCCTTGTATAAGGGAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGt   >  1:3893925/1‑69 (MQ=255)
                                                                   aagaag‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt  >  1:3216719/1‑70 (MQ=255)
                                                                    agaag‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                      >  1:2576634/1‑49 (MQ=255)
                                                                    agaag‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt  <  1:1974021/69‑1 (MQ=255)
                                                                    |                                                                     
GCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTT  >  minE/2045758‑2045894

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: ‑