Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,023,139 A→G 35.5% G256G (GGT→GGC ygiQ ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,023,1390AG35.5% 37.9 / 32.7 59G256G (GGT→GGCygiQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (17/21);  new base G (13/8);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

CGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACGGGTGGAATCCACGCCG  >  minE/2023072‑2023202
                                                                   |                                                               
cGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCAc                                                                <  1:1925539/69‑1 (MQ=255)
cGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCAc                                                                <  1:1513785/69‑1 (MQ=255)
cGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCAc                                                                <  1:437382/69‑1 (MQ=255)
cGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCAc                                                                <  1:551504/69‑1 (MQ=255)
  gTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCTGTTTGTTATACGCGCACGGAAAATCTTCAcc                                                               >  1:420874/1‑68 (MQ=255)
  gTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCa                                                              >  1:1383836/1‑69 (MQ=255)
  gTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCa                                                              >  1:2318451/1‑69 (MQ=255)
  gTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCa                                                              >  1:1779762/1‑69 (MQ=255)
  gTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCa                                                              >  1:1594830/1‑69 (MQ=255)
                         tttCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCggg                                                  >  1:2265330/1‑58 (MQ=255)
                         tttCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCggg                                                  >  1:189657/1‑58 (MQ=255)
                         tttCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCggg                                                  >  1:3307570/1‑58 (MQ=255)
                         tttCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCggg                                                  >  1:211189/1‑58 (MQ=255)
                         tttCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCggg                                                  <  1:707192/58‑1 (MQ=255)
                         tttCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCggg                                                  >  1:1555454/1‑58 (MQ=255)
                         tttCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCggg                                                  <  1:1409387/58‑1 (MQ=255)
                         tttCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCggg                                                  >  1:1309251/1‑58 (MQ=255)
                                   aCCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCggg                                                  <  1:1370021/48‑1 (MQ=255)
                                    ccGGTTTGTTATCCGCGCCCGGCAAATCTTCGCCATACGGATGCGGGATTGGGTCGAtttt                                    >  1:303452/1‑61 (MQ=255)
                                    ccGGTTTGTTATCCGCGCACGGCAAATCTTCGCCATa                                                            >  1:2967924/1‑37 (MQ=255)
                                    ccGGTTTGTTATCCGCGCACGGCAAATCTTCGCCATACGGATGCGGGATTGGGTCGAtttt                                    >  1:3137451/1‑61 (MQ=255)
                                    ccGGTTTGTTATCCGCGCACGGCAAATCTTCGCCATACGGATGCGGGATTGGGTCGAtttt                                    >  1:350383/1‑61 (MQ=255)
                                    ccGGTTTGTTATCCGCGCACGGCAAATCTTCGCCATACGGATGCGGGATTGGGTCGAtttt                                    >  1:3123423/1‑61 (MQ=255)
                                    ccGGTTTGTTATCCGCGCACGGCAAATCTTCGCCATACGGATGCGGGATTGGGTCGAtttt                                    >  1:662499/1‑61 (MQ=255)
                                    ccGGTTTGTTATCCGCGCACGGCAAATCTTCGCCATACGGATGCGGGATTGGGTCGAtttt                                    >  1:3120473/1‑61 (MQ=255)
                                    ccGGTTTGTTATCCGCGCACGGCAAATCTTCGCCATACGGATGCGGGATTGGGTCGAtttt                                    >  1:744122/1‑61 (MQ=255)
                                    ccGGTTTGTTATCCGCGCACGGCAAATCTTCGCCATACGGATGCGGGATTGGGTCGAtttt                                    >  1:3099763/1‑61 (MQ=255)
                                    ccGGTTTGTTATCCGCGCACGGCAAATCTTCGCCATACGGATGCGGGATTGGGTCGAtttt                                    >  1:95525/1‑61 (MQ=255)
                                    ccGGTTTGTTATCCGCGCACGGCAAATCTTCGCCATACGGATGCGGGATTGGGTCGAtttt                                    >  1:2696314/1‑61 (MQ=255)
                                    ccGGTTTGTTATCCGCGCACGGCAAATCTTCGCCATACGGATGCGGGATTGGGTCGAtttt                                    >  1:2389713/1‑61 (MQ=255)
                                    ccGGTTTGTTATCCGCGCACGGCAAATCTTCGCCATACGGATGCGGGATTGGGTCGAtttt                                    >  1:1655375/1‑61 (MQ=255)
                                          tGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCaa                      <  1:9185/69‑1 (MQ=255)
                                          tGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCaa                      <  1:3007426/69‑1 (MQ=255)
                                          tGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCaa                      <  1:2309829/69‑1 (MQ=255)
                                          tGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCaa                      <  1:483896/69‑1 (MQ=255)
                                          tGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCaa                      <  1:185998/69‑1 (MQ=255)
                                           gTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCaa                      <  1:238668/68‑1 (MQ=255)
                                              aTCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACg                  <  1:2764204/69‑1 (MQ=255)
                                              aTCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACg                  <  1:735356/69‑1 (MQ=255)
                                              aTCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACg                  <  1:2517064/69‑1 (MQ=255)
                                                       cGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACGGGTGGAATc         <  1:1790401/69‑1 (MQ=255)
                                                       cGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACGGGTGGAATc         <  1:1145868/69‑1 (MQ=255)
                                                        ggCATATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACGGGTGGAATc         <  1:1441901/68‑1 (MQ=255)
                                                          cAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCaa                      <  1:1074611/53‑1 (MQ=255)
                                                           aaaTCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACg                  <  1:1712904/56‑1 (MQ=255)
                                                            aaTCTTCGCCATACGGATGCGGGATTGGGTCGATTTTCCCTGGGGTATCAAGACGGGTGGAATCCAcac    <  1:3282424/69‑3 (MQ=255)
                                                            aaTCTTCGCCATACGGATGCGGGATTGGGTCGATTTTCCCTGGGGTATCAAGACGGGTGGAATCCAcac    <  1:1714302/69‑3 (MQ=255)
                                                            aaTCTTCGCCATACGGATGCGGGATTGGGTCGATTTTCCCTGGGGTATCAAGACGGGTGGAATCCAcac    <  1:1984730/69‑3 (MQ=255)
                                                            aaTCTTCGCCATACGGATGCGGGATTGGGTCGATTTTCCCTGGGGTATCAAGACGGGTGGAATCCAcac    <  1:488292/69‑3 (MQ=255)
                                                            aaTCTTCGCCATACGGATGCGGGATTGGGTCGATTTTCCCTGGGGTATCAAGACGGGTGGAATCCAcac    <  1:2847823/69‑3 (MQ=255)
                                                            aaTCTTCGCCATACGGATGCGGGATTGGGTCGATTTTCCCTGGGGTATCAAGACGGGTGGAATCCAcac    <  1:1591362/69‑3 (MQ=255)
                                                             aTCTTCGCCATACGGATGCGGGATTGGGTCGATTTTCCCTGGGGTATCAAGACGGGTGGAATCCAcac    <  1:2729898/68‑3 (MQ=255)
                                                             aTCTTCGCCATACGGATGCGGGATTGGGTCGATTTTCCCTGGGGTATCAAGACGGGTGGAATCCAcac    <  1:815131/68‑3 (MQ=255)
                                                              tCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGa                    >  1:3368074/1‑51 (MQ=255)
                                                              tCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACGGGTGGAATCCAcgccg  >  1:742739/1‑69 (MQ=255)
                                                              tCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACGGGTGGAATCCAcgccg  >  1:2849726/1‑69 (MQ=255)
                                                              tCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACGGGTGGAATCCAcgccg  >  1:156964/1‑69 (MQ=255)
                                                              tCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACGGGTGGAATCCAcgccg  >  1:580422/1‑69 (MQ=255)
                                                              tCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACGGGTGGAATCCAcgccg  >  1:2205513/1‑69 (MQ=255)
                                                                   |                                                               
CGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGGCAAATCTTCACCATACGGATGCGGGATTGGGTCGATTTTTCCAGGGGTATCAAGACGGGTGGAATCCACGCCG  >  minE/2023072‑2023202

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: