Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,623,335 T→C 50.7% A130A (GCT→GCC ftsY → fused Signal Recognition Particle (SRP) receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,623,3350TC50.7% ‑1.1 / 64.9 73A130A (GCT→GCCftsYfused Signal Recognition Particle (SRP) receptor
Reads supporting (aligned to +/- strand):  ref base T (13/23);  new base C (15/22);  total (28/45)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GTGAAGAGCTGCCGTTGCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCGGCTAAAGAAGAAA  >  minE/2623271‑2623393
                                                                |                                                          
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:2930981/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:2494099/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:2224001/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:2067477/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:2545367/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:2030346/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:3367364/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:1977305/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:607582/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:2944711/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:2653095/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:720852/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:1739485/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:788271/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:828096/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:2798502/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:27380/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:2453609/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:929604/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:2287278/68‑1 (MQ=255)
gTGAAGAACTGCCGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGACAGGCCGaa                                                         <  1:610215/68‑1 (MQ=255)
gTGAAGAACTGACGTTGTCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCCGaa                                                         <  1:2972576/68‑1 (MQ=255)
     gagCTGCCGTTGCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGa                                                    <  1:2779809/68‑1 (MQ=255)
     gagCTGCCGTTGCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGa                                                    <  1:750153/68‑1 (MQ=255)
       gCTGCCGTTGCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAAc                                                 >  1:2363530/1‑69 (MQ=255)
        cTGCCGTTGCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAAc                                                 >  1:1293118/1‑68 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGTAAACCg                                               <  1:2415232/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:2957658/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:3298514/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:2450169/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:2444283/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:3212998/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:2329260/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:2228621/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:213842/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:3311914/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:1905409/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:613070/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:713570/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:882606/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCg                                               <  1:998435/62‑1 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAgaga                                         >  1:862827/1‑68 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGAt                                        >  1:3285323/1‑69 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGAt                                        >  1:1089636/1‑69 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGAt                                        >  1:2217968/1‑69 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGAt                                        >  1:1957240/1‑69 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGAt                                        >  1:1937349/1‑69 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGAt                                        >  1:2002024/1‑69 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGAt                                        >  1:1876939/1‑69 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGAt                                        >  1:1730379/1‑69 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGAt                                        >  1:2247107/1‑69 (MQ=255)
                gCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGAt                                        >  1:2607974/1‑69 (MQ=255)
                                        ggTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaagaaga                   <  1:2752739/66‑1 (MQ=255)
                                         gTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCg               <  1:723024/69‑1 (MQ=255)
                                         gTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCg               <  1:285858/69‑1 (MQ=255)
                                         gTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCg               <  1:563528/69‑1 (MQ=255)
                                                  gaagaGTGGCAGGCTGAAGCGGAAAGCGTAGAGATTGTCTAAGCGGCGgaagaaga                   <  1:108021/56‑1 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAgc                               >  1:1433620/1‑44 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:170671/1‑52 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:1367695/1‑52 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:2717843/1‑52 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:740509/1‑52 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:1201031/1‑52 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:690103/1‑52 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:1899855/1‑52 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:3431439/1‑52 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:3350560/1‑52 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:3330709/1‑52 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:2467670/1‑52 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:2521185/1‑52 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:2709370/1‑52 (MQ=255)
                                                  gaagaGTGGCAGGCCGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGgaag                       >  1:2899657/1‑52 (MQ=255)
                                                    agagTGGCAGGCTGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCGGCTAAagaaga    <  1:983262/69‑1 (MQ=255)
                                                                tGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCGGCTAAAGAAGaaa  <  1:553564/59‑1 (MQ=255)
                                                                tGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCGGCTAAAGAAGaaa  <  1:699185/59‑1 (MQ=255)
                                                                tGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCGGCTAAAGAAGaaa  <  1:1892151/59‑1 (MQ=255)
                                                                tGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCGGCTAAAGAAGaaa  <  1:2056709/59‑1 (MQ=255)
                                                                tGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCGGCTAAAGAAGaaa  <  1:3309633/59‑1 (MQ=255)
                                                                tGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCGGCTAAAGAAGaaa  <  1:1747872/59‑1 (MQ=255)
                                                                tGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCGGCTAAAGAAGaaa  <  1:3301301/59‑1 (MQ=255)
                                                                tGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCGGCTAAAGAAGaaa  <  1:2405026/59‑1 (MQ=255)
                                                                |                                                          
GTGAAGAGCTGCCGTTGCCGGAAGACGTCAACGCCGAAGCGGTTTCGCCAGAAGAGTGGCAGGCTGAAGCGGAAACCGTAGAGATTGTCGAAGCGGCGGAAGAAGAAGCGGCTAAAGAAGAAA  >  minE/2623271‑2623393

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