Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,951,017 A→T 6.0% N172Y (AAT→TAT)  bglJ → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,951,0170AT 6.0% 340.7 / 10.6 117N172Y (AAT→TAT) bglJDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (62/47);  new base T (3/4);  total (66/51)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 125 total aligned reads displayed.

GGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATG  >  minE/2950948‑2951084
                                                                     |                                                                   
ggCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTaa                                                                    >  1:2369565/1‑71 (MQ=255)
ggCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTaa                                                                    >  1:959391/1‑71 (MQ=255)
ggCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTaa                                                                    >  1:2373773/1‑71 (MQ=255)
     cTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGa                                                               >  1:1549068/1‑71 (MQ=255)
      tCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGAt                                                              >  1:2804141/1‑71 (MQ=255)
        aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGt                                                            >  1:817898/1‑71 (MQ=255)
          tGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAAtgt                                                                 <  1:245845/64‑1 (MQ=255)
          tGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTc                                                           >  1:2493357/1‑70 (MQ=255)
           gACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCg                                                          >  1:3051355/1‑70 (MQ=255)
           gACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCg                                                          >  1:3055843/1‑70 (MQ=255)
            acacAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGa                                                               >  1:255619/1‑64 (MQ=255)
             cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGaaa                                                       <  1:3100491/71‑1 (MQ=255)
             cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGaaa                                                       >  1:2499140/1‑71 (MQ=255)
              acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAAc                                                      <  1:3317983/71‑1 (MQ=255)
              acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAAc                                                      >  1:345916/1‑71 (MQ=255)
              acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAAc                                                      >  1:1448833/1‑71 (MQ=255)
               caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACt                                                     <  1:1821489/71‑1 (MQ=255)
               caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACt                                                     >  1:254684/1‑71 (MQ=255)
               caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACt                                                     <  1:4525/71‑1 (MQ=255)
               caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACt                                                     <  1:707355/71‑1 (MQ=255)
                aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTg                                                    >  1:2143847/1‑71 (MQ=255)
                aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTg                                                    >  1:2148302/1‑71 (MQ=255)
                 aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTa                                                                     <  1:3518967/53‑1 (MQ=255)
                 aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTgg                                                   <  1:1569219/71‑1 (MQ=255)
                 aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTgg                                                   >  1:3003925/1‑71 (MQ=255)
                 aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTgg                                                   <  1:2763927/71‑1 (MQ=255)
                  aTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTa                                                                     <  1:2617405/52‑1 (MQ=255)
                      ccGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGt                                               >  1:3595935/1‑70 (MQ=255)
                      ccGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGt                                               >  1:3163739/1‑70 (MQ=255)
                          gcagcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTa                                                                     >  1:825325/1‑44 (MQ=255)
                          gcagcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGtt                                          >  1:1015467/1‑71 (MQ=255)
                          gcagcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGt                                           >  1:2695094/1‑70 (MQ=255)
                           cagcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGa                                                         <  1:1846873/55‑1 (MQ=255)
                           cagcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTc                                         >  1:2990725/1‑71 (MQ=255)
                            agcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTgg                                                   <  1:317946/60‑1 (MQ=255)
                             gcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCt                                        >  1:3017538/1‑70 (MQ=255)
                               ttAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGa                                                         <  1:427479/51‑1 (MQ=255)
                               ttAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGa                                                         <  1:2563989/51‑1 (MQ=255)
                               ttAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGAc                                     >  1:1272757/1‑71 (MQ=255)
                                tAAACGCAATATCAAAACGATCCGTGCACATAAATTTa                                                                     >  1:2654971/1‑38 (MQ=255)
                                tAAACGCAATATCAAAACGATCCGTGCACATAAATTTAAtgt                                                                 >  1:2237103/1‑42 (MQ=255)
                                tAAACGCAATATCAAAACGATCCGTGCACATAAATTTAAtgt                                                                 >  1:963073/1‑42 (MQ=255)
                                tAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACg                                    >  1:2155150/1‑71 (MQ=255)
                                     gCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAggg                               <  1:3189585/71‑1 (MQ=255)
                                     gCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAggg                               >  1:2495166/1‑71 (MQ=255)
                                     gCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAggg                               >  1:1254262/1‑71 (MQ=255)
                                     gCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAgg                                >  1:1925344/1‑70 (MQ=255)
                                       aaTATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGt                                               >  1:1359809/1‑53 (MQ=255)
                                       aaTATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCt                                        <  1:1354561/60‑1 (MQ=255)
                                       aaTATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCt                                        <  1:301359/60‑1 (MQ=255)
                                        atatCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttg                            <  1:3255054/71‑1 (MQ=255)
                                        atatCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttg                            <  1:1867583/71‑1 (MQ=255)
                                          atCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                          >  1:2939484/1‑71 (MQ=255)
                                           tCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                          <  1:2239390/70‑1 (MQ=255)
                                           tCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                         >  1:3600156/1‑71 (MQ=255)
                                           tCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                         >  1:2435603/1‑71 (MQ=255)
                                           tCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                         >  1:1283465/1‑71 (MQ=255)
                                           tCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                         <  1:877472/71‑1 (MQ=255)
                                             aaaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTgg                        <  1:2021589/70‑1 (MQ=255)
                                              aaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAg                      >  1:2603468/1‑71 (MQ=255)
                                               aaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGtt                                          <  1:3099319/50‑1 (MQ=255)
                                               aaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                     >  1:592060/1‑71 (MQ=255)
                                               aaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                     <  1:3626756/71‑1 (MQ=255)
                                               aaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                     <  1:2272111/71‑1 (MQ=255)
                                               aaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                     <  1:2029022/71‑1 (MQ=255)
                                                aCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGc                    >  1:3593364/1‑71 (MQ=255)
                                                 cGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCt                                        >  1:2955893/1‑50 (MQ=255)
                                                  gATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCg                  <  1:1929056/71‑1 (MQ=255)
                                                    tCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCg                  >  1:2355273/1‑69 (MQ=255)
                                                    tCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCg                  >  1:2615497/1‑69 (MQ=255)
                                                     ccGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCg                  >  1:1001391/1‑68 (MQ=255)
                                                     cagTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCg                  >  1:1146078/3‑68 (MQ=255)
                                                      cGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAg                                 <  1:1185474/52‑1 (MQ=255)
                                                      cGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCata              >  1:1411260/1‑69 (MQ=255)
                                                         tcacaTTAAATTTATGTGATGTCGAAACTGGGCGTCAGTTc                                         <  1:429708/40‑1 (MQ=21)
                                                         tcacaTTAAATTTATGTGATGTCGAAACTGGGCGTCAGTTc                                         <  1:3661211/40‑1 (MQ=21)
                                                         tcacaTTAAATTTATGTGATGTCGAAACTGGGCGTCAGTTc                                         <  1:1070937/40‑1 (MQ=21)
                                                         tcacaTTAAATTTATGTGATGTCGAAACTGGGCGTCAGTTc                                         <  1:3289853/40‑1 (MQ=21)
                                                         tcacaTTAAATTTATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGatat           >  1:2942933/2‑71 (MQ=255)
                                                         gCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGatat           >  1:2739620/1‑71 (MQ=255)
                                                         gCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGAACTGACg                                    <  1:916027/46‑1 (MQ=39)
                                                         gCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGAACTGACGcc                                  >  1:1330238/1‑47 (MQ=38)
                                                          cacaTAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGat             >  1:492199/1‑68 (MQ=255)
                                                           acaTTAAATTTATGTGATGTCGAAACTGGGCGTCAGTTCTg                                       >  1:2225819/1‑41 (MQ=25)
                                                           acaTTAAATTTATGTGATGTCGAAACTGGGCGTCAGTTCTg                                       >  1:2416462/1‑41 (MQ=25)
                                                            caTAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa                                      >  1:1543184/1‑41 (MQ=255)
                                                            caTAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGc                 <  1:2484211/62‑1 (MQ=255)
                                                            attAAATTTATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTg       >  1:1720524/3‑71 (MQ=255)
                                                             aTAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTg       <  1:1410261/71‑1 (MQ=255)
                                                               aaaTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                       <  1:742013/53‑1 (MQ=255)
                                                                aaTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                       <  1:964956/52‑1 (MQ=255)
                                                                aaTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTta    >  1:2968448/1‑71 (MQ=255)
                                                                 aTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                         >  1:2401507/1‑49 (MQ=255)
                                                                 aTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTc         <  1:1337176/65‑1 (MQ=255)
                                                                 aTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTc         >  1:188472/1‑65 (MQ=255)
                                                                 aTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTc         <  1:2131514/65‑1 (MQ=255)
                                                                 aTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtat   >  1:3365964/1‑71 (MQ=255)
                                                                 aTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtat   >  1:2797563/1‑71 (MQ=255)
                                                                  tttAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatg  <  1:1406613/71‑1 (MQ=255)
                                                                  tttAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatg  <  1:2537088/71‑1 (MQ=255)
                                                                     |                                                                   
GGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATG  >  minE/2950948‑2951084

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: