Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,303,886 G→A 10.4% intergenic (+59/‑125) 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,8860GA10.4% 148.3 / 6.0 48intergenic (+59/‑125)metJ/yiiXDNA‑binding transcriptional repressor, S‑adenosylmethionine‑binding/predicted peptidoglycan peptidase
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base A (3/2);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATA  >  minE/2303819‑2303956
                                                                   |                                                                        
gAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCa                                                                       >  1:2230404/1‑71 (MQ=255)
 aaTACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCa                                                                       >  1:1263848/1‑70 (MQ=255)
  aTACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCata                                                                     <  1:2136359/71‑1 (MQ=255)
     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)
                          cgGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCa                                                                       >  1:3918579/1‑45 (MQ=255)
                            gAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGtg                                                                        >  1:3016309/1‑41 (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:1570548/1‑70 (MQ=255)
                                                                tACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg      >  1:1282061/1‑70 (MQ=255)
                                                                  cGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGt                                     >  1:1149222/1‑37 (MQ=255)
                                                                   gTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATa  <  1:1090771/71‑1 (MQ=255)
                                                                   |                                                                        
GAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATAT‑T‑GCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATA  >  minE/2303819‑2303956

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