Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,303,925 C→A 14.2% intergenic (+98/‑86) 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,9250CA14.2% 136.5 / 8.3 50intergenic (+98/‑86)metJ/yiiXDNA‑binding transcriptional repressor, S‑adenosylmethionine‑binding/predicted peptidoglycan peptidase
Reads supporting (aligned to +/- strand):  ref base C (15/26);  new base A (3/4);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

TCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCTTTTTAT  >  minE/2303855‑2303994
                                                                        |                                                                     
tCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGc                                                                       >  1:475944/1‑71 (MQ=255)
    cTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGa                                                                   <  1:684360/71‑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:1282061/1‑70 (MQ=255)
                            tACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                            >  1:1570548/1‑70 (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:2483091/70‑1 (MQ=255)
                                                 ttGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGtt                         <  1:2077217/70‑1 (MQ=255)
                                                 ttGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGtt                         <  1:3316956/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)
                                                               aCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCg          >  1:116358/1‑71 (MQ=255)
                                                                 ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCt        >  1:3824046/1‑71 (MQ=255)
                                                                   gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCttt      >  1:3617957/1‑71 (MQ=255)
                                                                     aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCa             >  1:1760371/1‑62 (MQ=255)
                                                                       gCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCg          <  1:3133517/63‑1 (MQ=255)
                                                                        ccGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCTTTTtat  <  1:2631238/70‑1 (MQ=255)
                                                                        |                                                                     
TCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCTTTTTAT  >  minE/2303855‑2303994

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: