Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,303,888 C→G 10.4% intergenic (+61/‑123) 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,8880CG10.4% 144.7 / 5.8 48intergenic (+61/‑123)metJ/yiiXDNA‑binding transcriptional repressor, S‑adenosylmethionine‑binding/predicted peptidoglycan peptidase
Reads supporting (aligned to +/- strand):  ref base C (21/22);  new base G (3/2);  total (24/24)
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.25e-01

GAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATA  >  minE/2303819‑2303956
                                                                     |                                                                      
gAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCa                                                                       >  1:2230404/1‑71 (MQ=255)
 aaTACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCa                                                                       >  1:1263848/1‑70 (MQ=255)
  aTACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCata                                                                     <  1:2136359/71‑1 (MQ=255)
     cTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTg                                                                >  1:3421182/1‑71 (MQ=255)
        aaGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCaa                                                             >  1:1859100/1‑71 (MQ=255)
             aaaaaGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTAt                                                        <  1:701422/71‑1 (MQ=255)
              aaaaGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTAtt                                                       >  1:1728957/1‑71 (MQ=255)
                     cggcgCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAAttt                                                 <  1:2723425/70‑1 (MQ=255)
                     cggcgCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTg                                                <  1:1430083/71‑1 (MQ=255)
                          cgGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCa                                                                       >  1:3918579/1‑45 (MQ=255)
                               tGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTg                                      <  1:2176156/71‑1 (MQ=255)
                                ggAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAAtt                                                  <  1:331990/58‑1 (MQ=255)
                                   aTCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAgg                                  <  1:1914871/71‑1 (MQ=255)
                                    tCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGc                                 >  1:475944/1‑71 (MQ=255)
                                     cGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTg                                                >  1:3460736/1‑55 (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:552111/70‑1 (MQ=255)
                                           aTGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATa                           <  1:3063233/70‑1 (MQ=255)
                                            tGCGCTACGCTTATCAGGCCTACGTCATATTGCaa                                                             >  1:956310/1‑35 (MQ=255)
                                            tGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAg                         >  1:486897/1‑71 (MQ=255)
                                              cgcTACGCTTATCAGGCCTACGTCATATTGCAATTTATTg                                                      <  1:1361191/40‑1 (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:2156774/1‑71 (MQ=2)
                                               gcTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAgg                        >  1:1192829/1‑71 (MQ=2)
                                               gcTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAgg                        >  1:1526918/1‑71 (MQ=2)
                                                  aCGCTTATCAGGCCTACGTCATATTGCAATTTAtt                                                       >  1:2397392/1‑35 (MQ=255)
                                                  aCGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCa                                              >  1:417656/1‑44 (MQ=255)
                                                   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)
                                                         tCAGGCCTACGTCATATTGCAATTTATTGAATTTg                                                >  1:1223942/1‑35 (MQ=255)
                                                            ggCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATc         <  1:3201526/71‑1 (MQ=255)
                                                              ccTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCg                      <  1:3737738/56‑1 (MQ=255)
                                                               cTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATaa                          >  1:3401733/1‑51 (MQ=255)
                                                               cTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATaa                          >  1:3733583/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)
                                                                     |                                                                      
GAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATA  >  minE/2303819‑2303956

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: