Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,303,907 T→A 13.2% intergenic (+80/‑104) 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,9070TA13.2% 132.6 / 10.8 53intergenic (+80/‑104)metJ/yiiXDNA‑binding transcriptional repressor, S‑adenosylmethionine‑binding/predicted peptidoglycan peptidase
Reads supporting (aligned to +/- strand):  ref base T (19/27);  new base A (3/4);  total (22/31)
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

CGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGC  >  minE/2303840‑2303977
                                                                     |                                                                      
cggcgCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAAttt                                                                      <  1:2723425/70‑1 (MQ=255)
cggcgCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTg                                                                     <  1:1430083/71‑1 (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: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:1192829/1‑71 (MQ=2)
                          gcTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAgg                                             >  1:2156774/1‑71 (MQ=2)
                             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: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)
                                             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:2077217/70‑1 (MQ=255)
                                                                ttGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGtt        <  1:2483091/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)
                                                                     |                                                                      
CGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGC  >  minE/2303840‑2303977

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: