Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,303,913 C→T 13.8% intergenic (+86/‑98) 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,9130CT13.8% 137.8 / 10.2 51intergenic (+86/‑98)metJ/yiiXDNA‑binding transcriptional repressor, S‑adenosylmethionine‑binding/predicted peptidoglycan peptidase
Reads supporting (aligned to +/- strand):  ref base C (14/28);  new base T (3/4);  total (19/32)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACC  >  minE/2303850‑2303981
                                                                 |                                                                    
tGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTg                                                               <  1:2176156/71‑1 (MQ=255)
    aTCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAgg                                                           <  1:1914871/71‑1 (MQ=255)
     tCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGc                                                          >  1:475944/1‑71 (MQ=255)
         cTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGa                                                      <  1:684360/71‑1 (MQ=255)
          tGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACtt                                                                <  1:1659941/60‑1 (MQ=255)
            aTGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGaa                                                                   >  1:661777/1‑55 (MQ=255)
            aTGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTa                                                             <  1:2381663/61‑1 (MQ=255)
            aTGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATa                                                    <  1:3063233/70‑1 (MQ=255)
            aTGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATa                                                    <  1:552111/70‑1 (MQ=255)
             tGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAg                                                  >  1:486897/1‑71 (MQ=255)
               cgcTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATaa                                                   <  1:551920/70‑1 (MQ=2)
               cgcTACGCTTATCAAGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGc                                                <  1:2667106/71‑1 (MQ=255)
                gcTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATaa                                                   <  1:201637/69‑1 (MQ=2)
                gcTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAgg                                                 >  1:1526918/1‑71 (MQ=2)
                gcTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAgg                                                 >  1:2156774/1‑71 (MQ=2)
                gcTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAgg                                                 >  1:1192829/1‑71 (MQ=2)
                    cGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCa                                           >  1:2152599/1‑71 (MQ=255)
                     gCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCAc                                          <  1:1729040/71‑1 (MQ=255)
                      cTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCAcg                                         <  1:2442162/71‑1 (MQ=255)
                        tATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATa                                                    <  1:1503696/58‑1 (MQ=255)
                         aTCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGc                                                          >  1:2954192/1‑51 (MQ=255)
                         aTCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCg                                               <  1:1539513/62‑1 (MQ=255)
                             ggCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                  <  1:3201526/71‑1 (MQ=255)
                               ccTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCg                                               <  1:3737738/56‑1 (MQ=255)
                                cTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATaa                                                   >  1:3733583/1‑51 (MQ=255)
                                cTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATaa                                                   >  1:3401733/1‑51 (MQ=255)
                                 tACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCgg                                                       <  1:412488/46‑1 (MQ=255)
                                 tACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                               >  1:1570548/1‑70 (MQ=255)
                                 tACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                               >  1:1282061/1‑70 (MQ=255)
                                   cGTCATATTGCAATTTATTGAATTTGCACGAACTTGt                                                              >  1:1149222/1‑37 (MQ=255)
                                    gTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATa                           <  1:1090771/71‑1 (MQ=255)
                                       atatTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCAcgcc                                       <  1:2807952/56‑1 (MQ=255)
                                        tatTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaaca                       <  1:83032/71‑1 (MQ=255)
                                        tatTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaaca                       <  1:1638124/71‑1 (MQ=255)
                                        tatTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaac                        >  1:1253341/1‑70 (MQ=255)
                                          ttGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaacaac                     >  1:2957717/1‑71 (MQ=255)
                                          ttGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaacaac                     <  1:2315886/71‑1 (MQ=255)
                                              gCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAaaaaagc                      <  1:2064167/68‑4 (MQ=14)
                                               cAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAg                  >  1:1310132/1‑71 (MQ=14)
                                                 acttattaAAATTGCATGAACTTGTAGGACGGATAAGGCGTTCAc                                          <  1:3327258/43‑1 (MQ=14)
                                                  tttATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGc                 >  1:1598515/1‑69 (MQ=11)
                                                      ttGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGtt            <  1:3316956/70‑1 (MQ=255)
                                                      ttGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGtt            <  1:2483091/70‑1 (MQ=255)
                                                      ttGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGtt            <  1:2077217/70‑1 (MQ=255)
                                                       tGAATTTGCACGAACTTGTAGGCCGGATAAGGCGtt                                             >  1:3563336/1‑36 (MQ=40)
                                                       tGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCAc                                          <  1:3057011/39‑1 (MQ=40)
                                                       tGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGt          >  1:3129690/1‑71 (MQ=255)
                                                        gAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                  <  1:3421444/46‑1 (MQ=38)
                                                         aaTTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTc         <  1:3287368/70‑1 (MQ=255)
                                                           tttGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGc      <  1:403187/71‑1 (MQ=255)
                                                               cACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGAcc  <  1:659233/71‑1 (MQ=255)
                                                                 cGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTg           <  1:2362268/60‑1 (MQ=255)
                                                                 cGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGc      >  1:1357632/1‑65 (MQ=255)
                                                                 |                                                                    
TGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACC  >  minE/2303850‑2303981

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: