Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,624,178 G→T 47.2% A411A (GCG→GCT 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,1780GT47.2% 21.5 / 73.0 74A411A (GCG→GCTftsYfused Signal Recognition Particle (SRP) receptor
Reads supporting (aligned to +/- strand):  ref base G (9/30);  new base T (10/25);  total (19/55)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCGGGCAGAACGCGGTAAGCCAGGCCAAA  >  minE/2624112‑2624241
                                                                  |                                                               
cAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAgcgc                                                                >  1:1012086/1‑68 (MQ=255)
cAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAgcgc                                                                >  1:2362337/1‑68 (MQ=255)
cAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAgcgc                                                                >  1:2166325/1‑68 (MQ=255)
cAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAgcgc                                                                >  1:1616478/1‑68 (MQ=255)
cAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAgcgc                                                                >  1:1409643/1‑68 (MQ=255)
cAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAgcgc                                                                >  1:934424/1‑68 (MQ=255)
cAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAgcgc                                                                >  1:973896/1‑68 (MQ=255)
cAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAgcgc                                                                >  1:935385/1‑68 (MQ=255)
         ccTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGaatt                                                      <  1:3159177/69‑3 (MQ=255)
         ccTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGt                                                      <  1:2381798/69‑1 (MQ=255)
           tGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAATAAACTCGACGTTGAAGCGCCGCATGAAGTTa                                                    <  1:1892437/69‑1 (MQ=255)
           tGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTa                                                    <  1:475861/69‑1 (MQ=255)
             aTGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGt                                                      <  1:1119124/65‑1 (MQ=255)
               ggAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCt                                                <  1:2662333/69‑1 (MQ=255)
                gAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTg                                               >  1:2397924/1‑69 (MQ=255)
                  agagTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCt                                                <  1:2198144/66‑1 (MQ=255)
                  agagTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCt                                                <  1:2249454/66‑1 (MQ=255)
                  agagTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCt                                                <  1:2598778/66‑1 (MQ=255)
                  agagTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCt                                                <  1:2951572/66‑1 (MQ=255)
                  agagTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGAc                                             <  1:979257/69‑1 (MQ=255)
                  agagTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGAc                                             <  1:626748/69‑1 (MQ=255)
                  agagTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGAc                                             <  1:2770874/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:100819/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:984780/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:281797/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:300509/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:3072883/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:3076895/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:3244936/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:3428901/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:387526/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:41331/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:547849/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:907263/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:1631048/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:1880986/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:1811203/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:1179361/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:1196557/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:2129340/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:1062063/69‑1 (MQ=255)
                     gTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:2301214/69‑1 (MQ=255)
                      ttGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:2634572/68‑1 (MQ=255)
                      ttGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTAtt                                         <  1:91716/69‑1 (MQ=255)
                      ttGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTAtt                                         <  1:1832368/69‑1 (MQ=255)
                       tGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTAtt                                         <  1:1460569/68‑1 (MQ=255)
                        gaagaaAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTAtt                                         <  1:1329390/67‑1 (MQ=255)
                          agaaAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGAt                                      >  1:1237225/1‑68 (MQ=255)
                          agaaAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGAt                                      >  1:576371/1‑68 (MQ=255)
                          agaaAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGAt                                      >  1:3073715/1‑68 (MQ=255)
                          agaaAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGAt                                      >  1:30552/1‑68 (MQ=255)
                          agaaAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGAt                                      >  1:1917018/1‑68 (MQ=255)
                          agaaAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGAt                                      >  1:2367645/1‑68 (MQ=255)
                          agaaAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGAt                                      >  1:944945/1‑68 (MQ=255)
                            aaaaTCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTATTGATGcc                                   >  1:1974087/1‑69 (MQ=255)
                            aaaaTCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTATTGATGcc                                   >  1:2622780/1‑69 (MQ=255)
                            aaaaTCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTATTGATGcc                                   >  1:1687968/1‑69 (MQ=255)
                            aaaaTCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTATTGATGcc                                   >  1:2789964/1‑69 (MQ=255)
                            aaaaTCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTATTGATGcc                                   >  1:2436310/1‑69 (MQ=255)
                            aaaaTCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTATTGATGcc                                   >  1:1139652/1‑69 (MQ=255)
                            aaaaTCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTATTGATGcc                                   >  1:2982666/1‑69 (MQ=255)
                            aaaaTCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTATTGATGcc                                   >  1:2261990/1‑69 (MQ=255)
                            aaaaTCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTATTGATGc                                    >  1:1191582/1‑68 (MQ=255)
                            aaaaTCGTCCGCGTGATGAAGAAACTCGACGTTAAGCTCCGCATGAAGTTATGCTGACTATTGATGcc                                   >  1:1541508/1‑68 (MQ=255)
                             aaaTCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTa                                                    <  1:2710905/51‑1 (MQ=255)
                              aaTCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAg                                 <  1:656411/69‑1 (MQ=255)
                                  gTCCGCGTGATGAAGAAACTCGACGTTGAAGCTCCGCTTGAAGTTATGCTGACTAt                                          <  1:569361/56‑1 (MQ=255)
                                     cgcgTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCggg                          <  1:914230/69‑1 (MQ=255)
                                     cgcgTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCggg                          <  1:233760/69‑1 (MQ=255)
                                     cgcgTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCggg                          <  1:2065859/69‑1 (MQ=255)
                                             gtagaaACTCGTCGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:15200/43‑1 (MQ=37)
                                             gaagaaACTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:891896/45‑1 (MQ=255)
                                                gaaACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCGGGCAGAACGCGGt               <  1:904466/69‑1 (MQ=255)
                                                    cTCGACGTTGAAGCTCCGCATGAAGTTATGCTGACTAt                                          <  1:1828252/38‑1 (MQ=255)
                                                    cTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCGGGCAGAACGCTGTAAGc           <  1:2377269/69‑1 (MQ=255)
                                                    cTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCGGGCAGAACGCGGTAAGc           <  1:1654670/69‑1 (MQ=255)
                                                    cTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCGGGCAGAACGCGGTAAGc           <  1:979402/69‑1 (MQ=255)
                                                      cGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCGGGCAGAACGCGGt               <  1:3434235/63‑1 (MQ=255)
                                                           ttGATGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCggg                          <  1:2013711/47‑1 (MQ=255)
                                                           ttGAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCGGGCAGAACGCGGTAAGCCAGGCCa    <  1:1987227/69‑1 (MQ=255)
                                                           ttGAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCGGGCAGAACGCGGTAAGCCAGGCCa    <  1:1033873/69‑1 (MQ=255)
                                                             gAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCGGGCAGAACGCGGTAAGCCAGGCCaaa  >  1:3092914/1‑69 (MQ=255)
                                                                  |                                                               
CAAATCGCACCTGATGGAAGAGTTGAAGAAAATCGTCCGCGTGATGAAGAAACTCGACGTTGAAGCGCCGCATGAAGTTATGCTGACTATTGATGCCAGCACCGGGCAGAACGCGGTAAGCCAGGCCAAA  >  minE/2624112‑2624241

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: