Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,951,014 T→A 5.7% F171I (TTT→ATT)  bglJ → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,951,0140TA 5.7% 369.5 / 10.5 123F171I (TTT→ATT) bglJDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (68/48);  new base A (3/4);  total (71/52)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 127 total aligned reads displayed.

GTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATG  >  minE/2950946‑2951084
                                                                    |                                                                      
gTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAAttt                                                                      <  1:2817799/71‑1 (MQ=255)
 tGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAAttt                                                                      >  1:581073/1‑70 (MQ=255)
  ggCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTaa                                                                    >  1:959391/1‑71 (MQ=255)
  ggCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTaa                                                                    >  1:2369565/1‑71 (MQ=255)
   gCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTaa                                                                    <  1:2728189/70‑1 (MQ=255)
       cTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAAtgt                                                                 <  1:3258305/69‑1 (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)
               cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGaaa                                                       >  1:2499140/1‑71 (MQ=255)
                acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAAc                                                      <  1:3317983/71‑1 (MQ=255)
                acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAAc                                                      >  1:345916/1‑71 (MQ=255)
                 caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACt                                                     <  1:4525/71‑1 (MQ=255)
                 caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACt                                                     <  1:2315471/71‑1 (MQ=255)
                 caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACt                                                     <  1:707355/71‑1 (MQ=255)
                 caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACt                                                     <  1:1821489/71‑1 (MQ=255)
                 caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACt                                                     >  1:254684/1‑71 (MQ=255)
                 caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACt                                                     <  1:3510728/71‑1 (MQ=255)
                  aaaTTGCCTAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTg                                                    >  1:3623302/1‑71 (MQ=255)
                  aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTg                                                    <  1:2978378/71‑1 (MQ=255)
                  aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTg                                                    >  1:2143847/1‑71 (MQ=255)
                   aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTa                                                                     <  1:3518967/53‑1 (MQ=255)
                   aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTgg                                                   >  1:1456267/1‑71 (MQ=255)
                   aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTgg                                                   >  1:3003925/1‑71 (MQ=255)
                   aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTgg                                                   <  1:2763927/71‑1 (MQ=255)
                   aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTgg                                                   <  1:1569219/71‑1 (MQ=255)
                    aTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTa                                                                     <  1:2617405/52‑1 (MQ=255)
                        ccGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGt                                               >  1:3163739/1‑70 (MQ=255)
                           agcagcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGt                                           >  1:3026605/1‑71 (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:2412261/55‑1 (MQ=255)
                             cagcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGa                                                         <  1:1846873/55‑1 (MQ=255)
                             cagcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTc                                         >  1:2990725/1‑71 (MQ=255)
                              agcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTgg                                                   <  1:317946/60‑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:963073/1‑42 (MQ=255)
                                  tAAACGCAATATCAAAACGATCCGTGCACATAAATTTAAtgt                                                                 >  1:2237103/1‑42 (MQ=255)
                                  tAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACg                                    >  1:2155150/1‑71 (MQ=255)
                                   aaaCGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGAc                                     <  1:442014/69‑1 (MQ=255)
                                    aaCGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGaaa                                                       >  1:2545701/1‑50 (MQ=255)
                                       gCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAggg                               >  1:2490706/1‑71 (MQ=255)
                                       gCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAggg                               <  1:3189585/71‑1 (MQ=255)
                                       gCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAggg                               >  1:1254262/1‑71 (MQ=255)
                                       gCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAgg                                >  1:1925344/1‑70 (MQ=255)
                                         aaTATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTTGGCGTCAGTTCt                                        <  1:2223777/60‑1 (MQ=255)
                                         aaTATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGt                                               >  1:1359809/1‑53 (MQ=255)
                                         aaTATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCt                                        <  1:1354561/60‑1 (MQ=255)
                                         aaTATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCt                                        <  1:301359/60‑1 (MQ=255)
                                         aaTATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAggg                               >  1:1890011/1‑69 (MQ=255)
                                          atatCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttg                            <  1:1867583/71‑1 (MQ=255)
                                            atCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                          >  1:732793/1‑71 (MQ=255)
                                            atCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                          >  1:2939484/1‑71 (MQ=255)
                                             tCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                          <  1:3103913/70‑1 (MQ=255)
                                             tCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                         <  1:877472/71‑1 (MQ=255)
                                             tCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                         >  1:3600156/1‑71 (MQ=255)
                                             tCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                         >  1:1283465/1‑71 (MQ=255)
                                                aaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAg                      >  1:2603468/1‑71 (MQ=255)
                                                 aaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGtt                                          <  1:3099319/50‑1 (MQ=255)
                                                 aaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                     <  1:2029022/71‑1 (MQ=255)
                                                 aaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                     >  1:2607760/1‑71 (MQ=255)
                                                 aaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                     <  1:3626756/71‑1 (MQ=255)
                                                 aaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                     <  1:2272111/71‑1 (MQ=255)
                                                 aaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                     >  1:592060/1‑71 (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)
                                                       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:3289853/40‑1 (MQ=21)
                                                           tcacaTTAAATTTATGTGATGTCGAAACTGGGCGTCAGTTc                                         <  1:429708/40‑1 (MQ=21)
                                                           tcacaTTAAATTTATGTGATGTCGAAACTGGGCGTCAGTTc                                         <  1:3661211/40‑1 (MQ=21)
                                                           tcacaTTAAATTTATGTGATGTCGAAACTGGGCGTCAGTTc                                         <  1:1070937/40‑1 (MQ=21)
                                                           tcacaTTAAATTTATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGatat           >  1:2942933/2‑71 (MQ=255)
                                                           gCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGAACTGACGcc                                  >  1:1330238/1‑47 (MQ=38)
                                                             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)
                                                               aTAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTg       <  1:1410261/71‑1 (MQ=255)
                                                                  aaTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                       <  1:1739748/52‑1 (MQ=255)
                                                                  aaTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                       <  1:964956/52‑1 (MQ=255)
                                                                   aTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                         >  1:2401507/1‑49 (MQ=255)
                                                                   aTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                         >  1:888935/1‑49 (MQ=255)
                                                                   aTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATAtt          <  1:2094601/64‑1 (MQ=255)
                                                                   aTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTc         <  1:2131514/65‑1 (MQ=255)
                                                                   aTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTc         <  1:1337176/65‑1 (MQ=255)
                                                                   aTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtat   >  1:2797563/1‑71 (MQ=255)
                                                                    tttAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatg  <  1:1406613/71‑1 (MQ=255)
                                                                    tttAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatg  <  1:2537088/71‑1 (MQ=255)
                                                                    |                                                                      
GTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATG  >  minE/2950946‑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: