Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,045,839 T→A 9.3% intergenic (+601/+194) 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,8390TA 9.3% 241.7 / 13.6 86intergenic (+601/+194)yqiK/rfaEconserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base T (47/31);  new base A (3/5);  total (50/36)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTA  >  minE/2045771‑2045909
                                                                     |                                                                      
cgcTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtat                                                                      <  1:834630/71‑1 (MQ=255)
      tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGaa                                                                <  1:1971785/71‑1 (MQ=255)
      tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGaa                                                                <  1:1718862/71‑1 (MQ=255)
        aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAg                                                              >  1:2129772/1‑71 (MQ=255)
        aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAg                                                              >  1:317614/1‑71 (MQ=255)
            aTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGtt                                                          <  1:4136953/71‑1 (MQ=255)
            aTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGt                                                           <  1:3753043/70‑1 (MQ=255)
             tCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTa                                                         >  1:4079186/1‑71 (MQ=255)
             tCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTa                                                         >  1:1746001/1‑71 (MQ=255)
              cGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat                                                        >  1:3696942/1‑71 (MQ=255)
               ggggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat                                                        >  1:3568808/1‑70 (MQ=255)
               ggggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat                                                        >  1:1486467/1‑70 (MQ=255)
               ggggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat                                                        <  1:722908/70‑1 (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)
                                                        gaag‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc                                              <  1:37304/39‑1 (MQ=255)
                                                        gaag‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAg              >  1:4011946/1‑71 (MQ=255)
                                                        gaag‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAg              >  1:2986336/1‑71 (MQ=255)
                                                        gaag‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAg              >  1:1174285/1‑71 (MQ=255)
                                                        gaag‑GGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTcatcat                                      <  1:4103083/47‑1 (MQ=255)
                                                           g‑gggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATagt                                   >  1:123436/1‑45 (MQ=255)
                                                           g‑gggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACt                                   >  1:2752281/1‑47 (MQ=255)
                                                           g‑gggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACt                                   >  1:491167/1‑47 (MQ=255)
                                                           g‑gggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGttt           <  1:3006968/71‑1 (MQ=255)
                                                           g‑gggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGttt           >  1:2320496/1‑71 (MQ=255)
                                                           g‑gggCCTTGTATAAGGAAAGGGTTAGGATGAAGCTCGTCATCATACt                                   >  1:1259975/1‑47 (MQ=255)
                                                             gggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAg              <  1:2864633/67‑1 (MQ=255)
                                                             gggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTc          <  1:1175351/71‑1 (MQ=255)
                                                                ccTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTccc        >  1:428173/1‑70 (MQ=255)
                                                                  ttGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                                         <  1:3881834/35‑1 (MQ=255)
                                                                  ttGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg                            >  1:856458/1‑48 (MQ=255)
                                                                  ttGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg                            <  1:3264500/48‑1 (MQ=255)
                                                                  ttGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg                            >  1:2584406/1‑48 (MQ=255)
                                                                  ttGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGAc     >  1:2699984/1‑71 (MQ=255)
                                                                  ttGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGAc     >  1:1245287/1‑71 (MQ=255)
                                                                   tGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGAc     >  1:3612270/1‑70 (MQ=255)
                                                                   tGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACt    >  1:4139899/1‑71 (MQ=255)
                                                                    gTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGAc     >  1:1268584/1‑69 (MQ=255)
                                                                    gTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACt    >  1:917346/1‑70 (MQ=255)
                                                                     tataAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgtt                          <  1:3670603/47‑1 (MQ=255)
                                                                     tataAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgtt                          <  1:822598/47‑1 (MQ=255)
                                                                     tataAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgtt                          <  1:4076667/47‑1 (MQ=255)
                                                                     tataAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTa  <  1:2703701/71‑1 (MQ=255)
                                                                     tataAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTa  <  1:2365116/71‑1 (MQ=255)
                                                                     |                                                                      
CGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTA  >  minE/2045771‑2045909

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