Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 287,902 T→C 18.9% S182S (TCT→TCC lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE287,9020TC18.9% 126.0 / 19.2 84S182S (TCT→TCClonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (30/38);  new base C (2/14);  total (32/52)
Fisher's exact test for biased strand distribution p-value = 2.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGA  >  minE/287834‑287970
                                                                    |                                                                    
gacgaTCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCt                                                                      >  1:2019158/1‑69 (MQ=255)
gacgaTCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCt                                                                      >  1:1356809/1‑69 (MQ=255)
gacgaTCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCt                                                                      >  1:1428721/1‑69 (MQ=255)
gacgaTCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCt                                                                      >  1:1531605/1‑69 (MQ=255)
  cgaTCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTg                                                                     <  1:924234/68‑1 (MQ=255)
         gcgcGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGga                                                              <  1:1957119/68‑1 (MQ=255)
         gcgcGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGga                                                              <  1:902900/68‑1 (MQ=255)
         gcgcGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGga                                                              <  1:1685573/68‑1 (MQ=255)
         gcgcGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGga                                                              <  1:2141674/68‑1 (MQ=255)
               cTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTcc                                                       <  1:2492669/69‑1 (MQ=255)
                         ccATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTaacg                                             >  1:245703/1‑69 (MQ=255)
                         ccATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTaacg                                             >  1:954617/1‑69 (MQ=255)
                         ccATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTaacg                                             >  1:3336197/1‑69 (MQ=255)
                         ccATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTaacg                                             >  1:2453272/1‑69 (MQ=255)
                          cATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCCGTTCTGGAGATGt                                                         >  1:698230/1‑56 (MQ=255)
                          cATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGttcgt                                              >  1:854781/1‑64 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAAGGAACGt                                        >  1:2244096/1‑69 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGt                                        >  1:2304169/1‑69 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGt                                        >  1:2692656/1‑69 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGt                                        >  1:766180/1‑69 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGt                                        >  1:512170/1‑69 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGt                                        >  1:192598/1‑69 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGt                                        >  1:1002531/1‑69 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGt                                        >  1:1894703/1‑69 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGt                                        >  1:1704537/1‑69 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGt                                        >  1:1386881/1‑69 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGt                                        >  1:3370343/1‑69 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGt                                        >  1:1253783/1‑69 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGt                                        >  1:402377/1‑69 (MQ=255)
                              gctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGATCTGGAGATGTCCGACGTTAACGAACGt                                        >  1:3139136/1‑69 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                <  1:3224732/69‑1 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                <  1:1080755/69‑1 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                <  1:3247558/69‑1 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                >  1:488023/1‑69 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                >  1:3309122/1‑69 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                <  1:711586/69‑1 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                >  1:3394466/1‑69 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                <  1:2434497/69‑1 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                >  1:116176/1‑69 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                <  1:970271/69‑1 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                >  1:1290222/1‑69 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                <  1:1809356/69‑1 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                <  1:1924479/69‑1 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                >  1:1833543/1‑69 (MQ=255)
                                      tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                <  1:823834/69‑1 (MQ=255)
                                      tatGCCGCTGAAACTGGATGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAata                                >  1:323157/1‑69 (MQ=255)
                                      tatGCAGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTaacgaac                                          >  1:255515/1‑59 (MQ=255)
                                        tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTc                                       <  1:2537668/60‑1 (MQ=255)
                                        tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:635738/61‑1 (MQ=255)
                                        tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:696867/61‑1 (MQ=255)
                                        tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:923005/61‑1 (MQ=255)
                                        tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:808630/61‑1 (MQ=255)
                                        tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:807906/61‑1 (MQ=255)
                                        tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:1284243/61‑1 (MQ=255)
                                        tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:474732/61‑1 (MQ=255)
                                        tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:1263482/61‑1 (MQ=255)
                                        tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:3389199/61‑1 (MQ=255)
                                        tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:786253/61‑1 (MQ=255)
                                        tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:425127/61‑1 (MQ=255)
                                        tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:426827/61‑1 (MQ=255)
                                          ccGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCt                             <  1:541261/68‑1 (MQ=255)
                                          ccGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCt                             <  1:2906980/68‑1 (MQ=255)
                                          ccGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCt                             <  1:2939634/68‑1 (MQ=255)
                                          ccGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCt                             <  1:2740479/68‑1 (MQ=255)
                                          ccGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCt                             <  1:2695313/68‑1 (MQ=255)
                                          ccGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCt                             <  1:2363813/68‑1 (MQ=255)
                                               gAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGc                       <  1:917391/69‑1 (MQ=255)
                                               gAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGc                       <  1:3227748/69‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACTAACGTCTGGAATATCTGATGGCaa                     <  1:986803/68‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGTCaa                     <  1:1863134/68‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCaa                     <  1:336526/68‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCaa                     <  1:1664194/68‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCaa                     <  1:1782652/68‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCaa                     <  1:83783/68‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCaa                     <  1:1987061/68‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCaa                     <  1:808374/68‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCaa                     <  1:2409528/68‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCaa                     <  1:752490/68‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCaa                     <  1:2554330/68‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCaa                     <  1:994930/68‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCaa                     <  1:523880/68‑1 (MQ=255)
                                                  aCTGGCTGACAAACAGTCCGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCaa                     <  1:50056/68‑1 (MQ=255)
                                                       cTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:1881023/46‑1 (MQ=255)
                                                           cAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCt                                      <  1:1290135/42‑1 (MQ=255)
                                                                  tCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATc    >  1:994970/1‑69 (MQ=255)
                                                                  tCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATc    >  1:3264669/1‑69 (MQ=255)
                                                                   cTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATc    >  1:2046927/1‑68 (MQ=255)
                                                                   cTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCg   >  1:1537883/1‑69 (MQ=255)
                                                                    tGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGa  <  1:1552037/69‑1 (MQ=255)
                                                                    tGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGa  <  1:2767624/69‑1 (MQ=255)
                                                                    |                                                                    
GACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGA  >  minE/287834‑287970

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: