Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,624,109 G→A 32.7% Q388Q (CAG→CAA ftsY → fused Signal Recognition Particle (SRP) receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,624,1090GA32.7% 49.6 / 29.2 61Q388Q (CAG→CAAftsYfused Signal Recognition Particle (SRP) receptor
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base A (9/11);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CAGGCAGCTAAAGCGCGTAATATCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGC  >  minE/2624050‑2624177
                                                           |                                                                    
cAGGCAGCTAAAGCGCGTAATATCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCg                                                             >  1:436797/1‑69 (MQ=255)
cAGGCAGCTAAAGCGCGTAATATCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCg                                                             >  1:373587/1‑69 (MQ=255)
cAGGCAGCTAAAGCGCGTAATATCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCg                                                             >  1:1575538/1‑69 (MQ=255)
cAGGCAGCTAAAGCGCGTAATATCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCg                                                             >  1:1600882/1‑69 (MQ=255)
                    tatCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                                           <  1:2715806/67‑1 (MQ=255)
                      tggACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTgaag                                        >  1:1787864/3‑68 (MQ=255)
                      tCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTgaag                                        >  1:170073/1‑68 (MQ=255)
                      tCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTgaag                                        >  1:3333068/1‑68 (MQ=255)
                      tCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTgaag                                        >  1:2810155/1‑68 (MQ=255)
                      tCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGCAAGAGTTgaa                                         >  1:2564438/1‑67 (MQ=255)
                        gACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                                           <  1:2243010/63‑1 (MQ=255)
                        gACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                                           <  1:2761901/63‑1 (MQ=255)
                        gACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                                           <  1:2356700/63‑1 (MQ=255)
                        gACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                                           <  1:1854364/63‑1 (MQ=255)
                        gACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                                           <  1:3428339/63‑1 (MQ=255)
                        gACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                                           <  1:414359/63‑1 (MQ=255)
                         aCGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                                           <  1:1495269/62‑1 (MQ=255)
                           gTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATc                                  >  1:2273921/1‑69 (MQ=255)
                           gTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATc                                  >  1:1887181/1‑69 (MQ=255)
                           gTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATc                                  >  1:3123813/1‑69 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                                      <  1:2508005/57‑1 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                                      <  1:2147157/57‑1 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                                      <  1:615186/57‑1 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                                      <  1:1698645/57‑1 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                                      <  1:1645272/57‑1 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                                      <  1:876213/57‑1 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                                      <  1:975965/57‑1 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATc                                  >  1:3236790/1‑61 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt                          >  1:423563/1‑69 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt                          >  1:2213102/1‑69 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt                          >  1:2576444/1‑69 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt                          >  1:462058/1‑69 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt                          >  1:573523/1‑69 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt                          >  1:2006328/1‑69 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt                          >  1:797919/1‑69 (MQ=255)
                                   tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt                          >  1:1482378/1‑69 (MQ=255)
                                    gCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                                      <  1:1592200/56‑1 (MQ=255)
                                    gCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                                      <  1:302868/56‑1 (MQ=255)
                                     ccGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                                      <  1:686671/55‑1 (MQ=255)
                                     ccGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                                      <  1:2730018/55‑1 (MQ=255)
                                          accgCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaa                    >  1:726538/4‑68 (MQ=255)
                                          aCAGCCGGCCGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaa                    >  1:2557939/1‑68 (MQ=255)
                                          aCAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaa                    >  1:3285888/1‑68 (MQ=255)
                                          aCAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaa                    >  1:2597100/1‑68 (MQ=255)
                                          aCAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaa                    >  1:2903809/1‑68 (MQ=255)
                                          aCAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaa                    >  1:13274/1‑68 (MQ=255)
                                          aCAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaa                    >  1:206600/1‑68 (MQ=255)
                                             gCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaagaa                 <  1:3228752/68‑1 (MQ=255)
                                             gCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaagaa                 <  1:2648617/68‑1 (MQ=255)
                                             gCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaagaa                 <  1:707241/68‑1 (MQ=255)
                                             gCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaagaa                 <  1:2211580/68‑1 (MQ=255)
                                                ggACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTc             <  1:3408289/69‑1 (MQ=255)
                                                ggACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTc             <  1:1512828/69‑1 (MQ=255)
                                                ggACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTc             <  1:1731469/69‑1 (MQ=255)
                                                ggACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTc             <  1:1753197/69‑1 (MQ=255)
                                                ggACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTc             <  1:1849292/69‑1 (MQ=255)
                                                ggACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTc             <  1:208400/69‑1 (MQ=255)
                                                ggACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTc             <  1:246503/69‑1 (MQ=255)
                                                ggACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTc             <  1:2535894/69‑1 (MQ=255)
                                                ggACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTc             <  1:2627251/69‑1 (MQ=255)
                                                         cAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGaa    >  1:477383/1‑69 (MQ=255)
                                                           gAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTc             <  1:2708072/58‑1 (MQ=255)
                                                           gAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAgc  >  1:3313716/1‑69 (MQ=255)
                                                           gAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAgc  >  1:1146395/1‑69 (MQ=255)
                                                           |                                                                    
CAGGCAGCTAAAGCGCGTAATATCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGC  >  minE/2624050‑2624177

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: