Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,624,094 A→C 42.6% T383T (ACA→ACC 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,0940AC42.6% 6.0 / 46.8 47T383T (ACA→ACCftsYfused Signal Recognition Particle (SRP) receptor
Reads supporting (aligned to +/- strand):  ref base A (19/8);  new base C (9/11);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TTCGACGCCATTCAGGCAGCTAAAGCGCGTAATATCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAA  >  minE/2624038‑2624159
                                                        |                                                                 
ttCGACGCCATTCAGGCAGCTAAAGCGCGTAATATCGACGTCCTGATTGCCGATACAGCCGGACGcc                                                         >  1:504102/1‑67 (MQ=255)
            cAGGCAGCTAAAGCGCGTAATATCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCg                                           >  1:436797/1‑69 (MQ=255)
            cAGGCAGCTAAAGCGCGTAATATCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCg                                           >  1:1575538/1‑69 (MQ=255)
            cAGGCAGCTAAAGCGCGTAATATCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCg                                           >  1:1600882/1‑69 (MQ=255)
            cAGGCAGCTAAAGCGCGTAATATCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCg                                           >  1:373587/1‑69 (MQ=255)
                                tatCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                         <  1:2715806/67‑1 (MQ=255)
                                  tggACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTgaag                      >  1:1787864/3‑68 (MQ=255)
                                  tCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTgaag                      >  1:3333068/1‑68 (MQ=255)
                                  tCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTgaag                      >  1:2810155/1‑68 (MQ=255)
                                  tCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTgaag                      >  1:170073/1‑68 (MQ=255)
                                  tCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGCAAGAGTTgaa                       >  1:2564438/1‑67 (MQ=255)
                                    gACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                         <  1:1854364/63‑1 (MQ=255)
                                    gACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                         <  1:414359/63‑1 (MQ=255)
                                    gACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                         <  1:2761901/63‑1 (MQ=255)
                                    gACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                         <  1:3428339/63‑1 (MQ=255)
                                    gACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                         <  1:2243010/63‑1 (MQ=255)
                                    gACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                         <  1:2356700/63‑1 (MQ=255)
                                     aCGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTg                         <  1:1495269/62‑1 (MQ=255)
                                       gTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATc                >  1:3123813/1‑69 (MQ=255)
                                       gTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATc                >  1:2273921/1‑69 (MQ=255)
                                       gTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATc                >  1:1887181/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:1482378/1‑69 (MQ=255)
                                               tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt        >  1:462058/1‑69 (MQ=255)
                                               tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt        >  1:573523/1‑69 (MQ=255)
                                               tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt        >  1:423563/1‑69 (MQ=255)
                                               tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt        >  1:797919/1‑69 (MQ=255)
                                               tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt        >  1:2576444/1‑69 (MQ=255)
                                               tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt        >  1:2213102/1‑69 (MQ=255)
                                               tGCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGt        >  1:2006328/1‑69 (MQ=255)
                                                gCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                    <  1:302868/56‑1 (MQ=255)
                                                gCCGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                    <  1:1592200/56‑1 (MQ=255)
                                                 ccGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                    <  1:2730018/55‑1 (MQ=255)
                                                 ccGATACCGCCGGACGCCTGCAAAACAAATCGCACCTGATGGAAGAGTTgaagaa                    <  1:686671/55‑1 (MQ=255)
                                                      aCAGCCGGCCGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaa  >  1:2557939/1‑68 (MQ=255)
                                                      aCAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaa  >  1:3285888/1‑68 (MQ=255)
                                                      aCAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaa  >  1:2903809/1‑68 (MQ=255)
                                                      aCAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaa  >  1:2597100/1‑68 (MQ=255)
                                                      aCAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaa  >  1:206600/1‑68 (MQ=255)
                                                      aCAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATgaa  >  1:13274/1‑68 (MQ=255)
                                                        |                                                                 
TTCGACGCCATTCAGGCAGCTAAAGCGCGTAATATCGACGTCCTGATTGCCGATACAGCCGGACGCCTGCAGAACAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAA  >  minE/2624038‑2624159

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: