Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,811,594 G→A 43.6% R384R (CGC→CGT uvrA ← ATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,811,5940GA43.6% 13.8 / 39.8 55R384R (CGC→CGTuvrAATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC
Reads supporting (aligned to +/- strand):  ref base G (16/15);  new base A (14/10);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

GATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTT  >  minE/2811536‑2811662
                                                          |                                                                    
gATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCa                                                            <  1:2191493/69‑1 (MQ=255)
   tACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCatat                                                         <  1:1456921/69‑1 (MQ=255)
   tACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCatat                                                         <  1:2517117/69‑1 (MQ=255)
      tGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCAtatta                                                       >  1:1120772/1‑68 (MQ=255)
      tGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCAtatta                                                       <  1:766151/68‑1 (MQ=255)
         tAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATgc                                                    >  1:2775006/1‑68 (MQ=255)
         tAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATgc                                                    >  1:1757157/1‑68 (MQ=255)
                  cTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCt                                           <  1:1041072/68‑1 (MQ=255)
                  cTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCt                                           <  1:115925/68‑1 (MQ=255)
                  cTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCt                                           <  1:1326607/68‑1 (MQ=255)
                  cTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCt                                           <  1:320834/68‑1 (MQ=255)
                  cTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCt                                           <  1:3047428/68‑1 (MQ=255)
                  cTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCt                                           <  1:2928990/68‑1 (MQ=255)
                  cTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCt                                           <  1:2502816/68‑1 (MQ=255)
                    aaTTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCa                                                            >  1:982113/1‑49 (MQ=255)
                    aaTTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCa                                                            >  1:1039935/1‑49 (MQ=255)
                                   gcgcTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGacgacg                          <  1:2371332/68‑1 (MQ=255)
                                   gcgcTGGATTCCATCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGacgacg                          <  1:2767532/68‑1 (MQ=255)
                                         gATTCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCt                                           <  1:2498738/45‑1 (MQ=255)
                                           ttCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAgg                  >  1:1776710/1‑68 (MQ=255)
                                           ttCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAgg                  <  1:778445/68‑1 (MQ=255)
                                           ttCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGt                 >  1:423855/1‑69 (MQ=255)
                                           ttCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGt                 >  1:554396/1‑69 (MQ=255)
                                             ccGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTAt               <  1:2785796/69‑1 (MQ=255)
                                             ccGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTAt               <  1:3189333/69‑1 (MQ=255)
                                             ccGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTAt               <  1:498283/69‑1 (MQ=255)
                                             ccGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTAt               <  1:2764151/69‑1 (MQ=255)
                                             ccGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTAt               <  1:2739021/69‑1 (MQ=255)
                                             ccGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTAt               <  1:2606777/69‑1 (MQ=255)
                                             ccGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTAt               <  1:2567996/69‑1 (MQ=255)
                                             ccGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTAt               <  1:3330843/69‑1 (MQ=255)
                                             ccGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTAt               <  1:2418045/69‑1 (MQ=255)
                                             ccGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTAt               <  1:87461/69‑1 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTAt               >  1:3237771/1‑68 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:2768214/1‑69 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:610443/1‑69 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:1481900/1‑69 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:487589/1‑69 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:454597/1‑69 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:1680359/1‑69 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:2913689/1‑69 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:734767/1‑69 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:2788922/1‑69 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:842424/1‑69 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:2460450/1‑69 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:1914988/1‑69 (MQ=255)
                                              cGTCTCTTTATAACGGCGCTCCATATTGTGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:210686/1‑69 (MQ=255)
                                               gTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:615135/1‑68 (MQ=255)
                                               gTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:1454459/1‑68 (MQ=255)
                                               gTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:1316166/1‑68 (MQ=255)
                                               gTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:836754/1‑68 (MQ=255)
                                               gTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:3341511/1‑68 (MQ=255)
                                               gTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:3046514/1‑68 (MQ=255)
                                               gTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:2519227/1‑68 (MQ=255)
                                               gTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATc              >  1:2462184/1‑68 (MQ=255)
                                                          gcggcgCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGtt  <  1:3094997/69‑1 (MQ=255)
                                                          gcggcgCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGtt  <  1:1619544/69‑1 (MQ=255)
                                                          gcggcgCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGtt  >  1:2930813/1‑69 (MQ=255)
                                                          gcggcgCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGtt  >  1:2862130/1‑69 (MQ=255)
                                                          gcggcgCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGtt  >  1:1712104/1‑69 (MQ=255)
                                                          gcggcgCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGtt  >  1:799291/1‑69 (MQ=255)
                                                          gcggcgCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGtt  <  1:1784853/69‑1 (MQ=255)
                                                          gcggcgCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGtt  >  1:2280773/1‑69 (MQ=255)
                                                          gcggcgCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGtt  >  1:2036912/1‑69 (MQ=255)
                                                          |                                                                    
GATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTT  >  minE/2811536‑2811662

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: