Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,303,892:1 +C 11.4% intergenic (+65/‑119) metJ → / → yiiX DNA‑binding transcriptional repressor, S‑adenosylmethionine‑binding/predicted peptidoglycan peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,303,8921.C11.4% 159.2 / 19.8 44intergenic (+65/‑119)metJ/yiiXDNA‑binding transcriptional repressor, S‑adenosylmethionine‑binding/predicted peptidoglycan peptidase
Reads supporting (aligned to +/- strand):  ref base . (18/21);  new base C (3/2);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAAC  >  minE/2303824‑2303962
                                                                     |                                                                       
cTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑g                                                                      >  1:3421182/1‑71 (MQ=255)
   aaGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCaa                                                                   >  1:1859100/1‑71 (MQ=255)
        aaaaaGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTAt                                                              <  1:701422/71‑1 (MQ=255)
         aaaaGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTAtt                                                             >  1:1728957/1‑71 (MQ=255)
                cggcgCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAAttt                                                       <  1:2723425/70‑1 (MQ=255)
                cggcgCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTg                                                      <  1:1430083/71‑1 (MQ=255)
                          tGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTg                                            <  1:2176156/71‑1 (MQ=255)
                           ggAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAAtt                                                        <  1:331990/58‑1 (MQ=255)
                              aTCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAgg                                        <  1:1914871/71‑1 (MQ=255)
                               tCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGc                                       >  1:475944/1‑71 (MQ=255)
                                cGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTg                                                      >  1:3460736/1‑55 (MQ=255)
                                   cTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGa                                   <  1:684360/71‑1 (MQ=255)
                                    tGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACtt                                             <  1:1659941/60‑1 (MQ=255)
                                      aTGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGaa                                                >  1:661777/1‑55 (MQ=255)
                                      aTGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTa                                          <  1:2381663/61‑1 (MQ=255)
                                      aTGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATa                                 <  1:552111/70‑1 (MQ=255)
                                      aTGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATa                                 <  1:3063233/70‑1 (MQ=255)
                                       tGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCaa                                                                   >  1:956310/1‑35 (MQ=255)
                                       tGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAg                               >  1:486897/1‑71 (MQ=255)
                                         cgcTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTg                                                            <  1:1361191/40‑1 (MQ=255)
                                         cgcTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATaa                                <  1:551920/70‑1 (MQ=2)
                                         cgcTACGCTTATCAAGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGc                             <  1:2667106/71‑1 (MQ=255)
                                          gcTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATaa                                <  1:201637/69‑1 (MQ=2)
                                          gcTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAgg                              >  1:2156774/1‑71 (MQ=2)
                                          gcTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAgg                              >  1:1526918/1‑71 (MQ=2)
                                          gcTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAgg                              >  1:1192829/1‑71 (MQ=2)
                                             aCGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTAtt                                                             >  1:2397392/1‑35 (MQ=255)
                                             aCGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCa                                                    >  1:417656/1‑44 (MQ=255)
                                              cGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCa                        >  1:2152599/1‑71 (MQ=255)
                                               gCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCAc                       <  1:1729040/71‑1 (MQ=255)
                                                cTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCAcg                      <  1:2442162/71‑1 (MQ=255)
                                                  tATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATa                                 <  1:1503696/58‑1 (MQ=255)
                                                   aTCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGc                                       >  1:2954192/1‑51 (MQ=255)
                                                   aTCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCg                            <  1:1539513/62‑1 (MQ=255)
                                                    tCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTg                                                      >  1:1223942/1‑35 (MQ=255)
                                                       ggCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATc               <  1:3201526/71‑1 (MQ=255)
                                                         ccTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCg                            <  1:3737738/56‑1 (MQ=255)
                                                          cTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATaa                                >  1:3401733/1‑51 (MQ=255)
                                                          cTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATaa                                >  1:3733583/1‑51 (MQ=255)
                                                           tACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCgg                                    <  1:412488/46‑1 (MQ=255)
                                                           tACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg            >  1:1282061/1‑70 (MQ=255)
                                                           tACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg            >  1:1570548/1‑70 (MQ=255)
                                                             cGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGt                                           >  1:1149222/1‑37 (MQ=255)
                                                              gTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATa        <  1:1090771/71‑1 (MQ=255)
                                                                 atat‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCAcgcc                    <  1:2807952/56‑1 (MQ=255)
                                                                  tat‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaaca    <  1:1638124/71‑1 (MQ=255)
                                                                  tat‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaaca    <  1:83032/71‑1 (MQ=255)
                                                                  tat‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaac     >  1:1253341/1‑70 (MQ=255)
                                                                    t‑t‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaacaac  <  1:2315886/71‑1 (MQ=255)
                                                                    t‑t‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaacaac  >  1:2957717/1‑71 (MQ=255)
                                                                     |                                                                       
CTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAAC  >  minE/2303824‑2303962

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