Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,026,414 G→A 28.2% intergenic (‑173/+61) plsC ← / ← parC 1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/DNA topoisomerase IV, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,026,4140GA28.2% 68.3 / 29.6 46intergenic (‑173/+61)plsC/parC1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/DNA topoisomerase IV, subunit A
Reads supporting (aligned to +/- strand):  ref base G (10/23);  new base A (6/7);  total (16/30)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

CATGAATAAACAACGGTTTCAAGGGAATTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGCCCTCCCCGGCATAT  >  minE/2026347‑2026475
                                                                   |                                                             
cATGAATAAACAACGGTTTCAAGGGAATTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACAACt                                                            <  1:2410535/71‑1 (MQ=255)
  tGAATAAACAACGGTTTCAAGGGAATTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTa                                                          <  1:2697975/71‑1 (MQ=255)
    aaTAAACAACGGTTTCAAGGGAATTTTTATCGTTTCCGGGTGTTGAATACGCTGCCGGATGACAACTTaa                                                         <  1:1071133/70‑1 (MQ=255)
     aTAAACAACGGTTTCAAGGGAATTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACg                                                       <  1:2752891/71‑1 (MQ=255)
     aTAAACAACGGTTTCAAGGGAATTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACg                                                       >  1:3457167/1‑71 (MQ=255)
         acaacGGTTTCAAGGGAATTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTc                                                   <  1:1048795/71‑1 (MQ=255)
          caacGGTTTCAAGGGAATTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCa                                                  <  1:2172273/71‑1 (MQ=255)
              ggTTTCAAGGGAATTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCg                                              <  1:55748/71‑1 (MQ=255)
                    aaGGGAATTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACAACTTAACGTTTCATCCGGCGTTc                                        >  1:3376405/1‑71 (MQ=255)
                    aaGGGAATTTTTATCGTTTCC‑‑GGGTTGAATACGCTGCCGGATGACAACTTAACGTTTCATCCGGCGTTcc                                       <  1:712801/70‑1 (MQ=255)
                     aGGGAATTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTcc                                       <  1:2173621/71‑1 (MQ=255)
                        gAATTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTg                                    >  1:1633426/1‑71 (MQ=255)
                         aaTTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTg                                    >  1:1241884/1‑70 (MQ=255)
                          aTTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGt                                                      >  1:1848779/1‑51 (MQ=255)
                          aTTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGt                                                      >  1:1118640/1‑51 (MQ=255)
                           tttttATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACAACTTAACGTTTCATCCGGCGTTCCTTGc                                   >  1:738873/1‑69 (MQ=255)
                             tttATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAg                                <  1:3878152/70‑1 (MQ=255)
                             tttATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAg                                <  1:333794/70‑1 (MQ=255)
                                 tCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCggg                            <  1:571258/70‑1 (MQ=255)
                                 tCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCggg                            <  1:3331750/70‑1 (MQ=255)
                                   gTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGt                                          >  1:2424388/1‑54 (MQ=255)
                                   gTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGggag                          <  1:1695605/70‑1 (MQ=255)
                                   gTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGggag                          <  1:3401091/70‑1 (MQ=255)
                                    tttCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGggag                          <  1:1999039/69‑1 (MQ=255)
                                    tttCCGGGGGTTGAATACGCTGCCGGATGACAACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGggagga                        <  1:2390541/71‑1 (MQ=255)
                                      tCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGaaa                      >  1:3035065/1‑71 (MQ=255)
                                        cGGGGGTTGAATACGCTGCCGGATGACAACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGaaa                      <  1:196524/69‑1 (MQ=255)
                                         gggggTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAg                   <  1:1680763/71‑1 (MQ=255)
                                          ggggTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAg                   >  1:2659607/1‑70 (MQ=255)
                                          ggggTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAg                   <  1:2390752/70‑1 (MQ=255)
                                              ttGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGccc              <  1:139113/71‑1 (MQ=255)
                                                  aTACGCTGCCGGATGACAACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGCCCTccc          >  1:834232/1‑71 (MQ=255)
                                                   tACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGggagga                        <  1:356957/56‑1 (MQ=255)
                                                   tACGCTGCCGGATGACAACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGCCCTcccc         >  1:1109005/1‑71 (MQ=255)
                                                    aCGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGCCCTCCCCg        >  1:2887142/1‑71 (MQ=255)
                                                     cGCTGCCGGATGACAACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGCCCTCCCCg        <  1:997647/70‑1 (MQ=255)
                                                        tGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGggagga                        <  1:3138873/51‑1 (MQ=255)
                                                         gCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGCCCTCCCCGGCata   >  1:456757/1‑71 (MQ=255)
                                                         gCCGGATGACAACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGCCCTCCCCGGCata   >  1:827867/1‑71 (MQ=255)
                                                          ccGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGgga                           <  1:3185301/46‑1 (MQ=255)
                                                          ccGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGCCCTCCCCGGCatat  <  1:362092/71‑1 (MQ=255)
                                                          ccGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGCCCTCCCCGGCatat  <  1:748149/71‑1 (MQ=255)
                                                          ccGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGCCCTCCCCGGCatat  <  1:752761/71‑1 (MQ=255)
                                                          ccGGATGACAACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGCCCTCCCCGGCatat  >  1:3854940/1‑71 (MQ=255)
                                                           cGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGCCCTCCCCgg       <  1:2628094/65‑1 (MQ=255)
                                                               tGACAACTTAACGTTTCATCCGGCGTTCCTTGCAAg                                <  1:3617025/36‑1 (MQ=255)
                                                                   |                                                             
CATGAATAAACAACGGTTTCAAGGGAATTTTTATCGTTTCCGGGGGTTGAATACGCTGCCGGATGACGACTTAACGTTTCATCCGGCGTTCCTTGCAAGCGGGAGGAAACAGCGCCCTCCCCGGCATAT  >  minE/2026347‑2026475

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: