Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,303,856 C→T 17.0% intergenic (+29/‑155) 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,8560CT17.0% 101.9 / 13.1 47intergenic (+29/‑155)metJ/yiiXDNA‑binding transcriptional repressor, S‑adenosylmethionine‑binding/predicted peptidoglycan peptidase
Reads supporting (aligned to +/- strand):  ref base C (14/25);  new base T (3/5);  total (17/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 = 4.93e-01

GTGAGATGGGGATTAACCCGGAGACGTGGGAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTA  >  minE/2303790‑2303922
                                                                  |                                                                  
gTGAGATGGGGATTAACCCGGAGACGTGGGAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCt                                                                <  1:1562448/71‑1 (MQ=255)
 tGAGATGGGGATTAACCCGGAGACGTGGGAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCt                                                                >  1:2144752/1‑70 (MQ=255)
       ggggATTAACCCGGAGACGTGGGAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGATTCGCCTGATgcg                                                          >  1:3877934/1‑70 (MQ=255)
             tAACCCGGAGACGTGGGAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCT                                                   >  1:1229701‑M1/1‑71 (MQ=255)
             tAACCCGGAGACGTGGGAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCT                                                   >  1:2036374‑M1/1‑71 (MQ=255)
              aaCCCGGAGACGTGGGAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCT                                                   >  1:5656704‑M1/1‑70 (MQ=255)
               aCCCGGAGACGTGGGAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTA                                                 >  1:3037683‑M1/1‑71 (MQ=255)
               aCCCGGAGACGTGGGAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTT                                                  <  1:4297467‑M1/70‑1 (MQ=255)
                  cGGAGACGTGGGAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCA                                              <  1:3562871‑M1/71‑1 (MQ=255)
                      gaCGTGGGAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCt                                                        <  1:2344273/57‑1 (MQ=255)
                            ggAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGt                                     >  1:2796021/1‑70 (MQ=255)
                             gAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCa                                   <  1:308941/71‑1 (MQ=255)
                                tACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTA                                        >  1:5352514‑M1/1‑63 (MQ=255)
                                tACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTT                                                  <  1:207959‑M1/53‑1 (MQ=255)
                                   tAAAGCAAAAAAGACCGGCGCTGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCat                                  <  1:787971/66‑1 (MQ=255)
                                     aaGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCa                            <  1:2753124/70‑1 (MQ=255)
                                       gCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAAtt                         >  1:2785279/1‑71 (MQ=255)
                                         aaaaaaGAGCGGCGGGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTa                       <  1:1439509/71‑1 (MQ=255)
                                          aaaaaGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTA                                        >  1:1297954‑M1/1‑53 (MQ=255)
                                           aaaaGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTAtt                     >  1:560992/1‑71 (MQ=255)
                                            aaaGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTAtt                     >  1:446652/1‑70 (MQ=255)
                                            aaaGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTg                    <  1:456896/71‑1 (MQ=255)
                                              agagCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTg                    <  1:3538734/69‑1 (MQ=255)
                                              agagCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTg                    <  1:2476149/69‑1 (MQ=255)
                                               gagCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAAt                 <  1:3732742/71‑1 (MQ=255)
                                                       cgGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTa                       <  1:3937039/57‑1 (MQ=255)
                                                          gctggAATTGCCTGATGCGCTACGCTTATCAGGCCTACGagttcagtgctctacgtaggccggataaggcg      <  1:423301‑M1/69‑33 (MQ=255)
                                                           gTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCa                                   <  1:3706707/41‑1 (MQ=255)
                                                           gTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACtt     <  1:4268655/71‑1 (MQ=255)
                                                            tGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCat                                  <  1:182463/41‑1 (MQ=255)
                                                            tGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAAATTATTGAATTTGCACGAACTTg    <  1:197399/71‑1 (MQ=255)
                                                             ggAATTGCCTGATGCGCTACGCTTATCAGGCCTACGagttcagtgctctacgtaggccggataaggcgtt    >  1:1703141‑M1/1‑36 (MQ=255)
                                                             ggAATTGCCTGATGCGCTACGCTTATCAGGCCTACGagttcagtgctctacgtaggccggataaggcgtt    >  1:4305531‑M1/1‑36 (MQ=255)
                                                             ggAATTGCCTGATGCGCTACGCTTATCAGGCCTACGagttcagtgctctacgtaggccggataaggc       <  1:2998829‑M1/67‑32 (MQ=255)
                                                              gAATTGCCTGATGCGCTACGCTTATCAGGCCTACGagttcagtgctctacgtaggccggataaggcgttca  >  1:2826639‑M1/1‑35 (MQ=255)
                                                              gAATTGCCTGATGCGCTACGCTTATCAGGCCTACGagttcagtgctctacgtaggccggataaggcgtt    <  1:1040741‑M1/69‑35 (MQ=255)
                                                              gAATTGCCTGATGCGCTACGCTTATCAGGCCTACGagttcagtgctctacgtaggccggataaggcgtt    <  1:2241225‑M1/69‑35 (MQ=255)
                                                              gAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTg              <  1:1978944/59‑1 (MQ=255)
                                                              gAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGt   <  1:4594440/70‑1 (MQ=255)
                                                              gAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGt   <  1:1434517/70‑1 (MQ=255)
                                                              gAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGt   <  1:843407/70‑1 (MQ=255)
                                                              gAATCGCCTGATGCGCTACGCTTATCAGACCTACGTCATATTGCAATTTa                       >  1:674553/1‑50 (MQ=255)
                                                               aaTTGCCTGATGCGCTACGCTTATCAGGCCTACGagttcagtgctctacgtaggccggataaggcgtt    <  1:4719136‑M1/68‑35 (MQ=255)
                                                               aaTCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAAt                 <  1:4816554/55‑1 (MQ=255)
                                                               aaTCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTg    <  1:4278153/68‑1 (MQ=255)
                                                                 tCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGa                   <  1:1657780/51‑1 (MQ=255)
                                                                 tCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTACTGAATTTg              >  1:327091/1‑56 (MQ=255)
                                                                  |                                                                  
GTGAGATGGGGATTAACCCGGAGACGTGGGAATACTAAAGCAAAAAAGAGCGGCGCGGAGTGGAATCGCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTA  >  minE/2303790‑2303922

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: