Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 145,530 T→A 60.4% T233S (ACA→TCA)  sfsA ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE145,5300TA60.4% 36.3 / 61.6 53T233S (ACA→TCA) sfsApredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (14/7);  new base A (11/21);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 2.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTT  >  minE/145466‑145595
                                                                |                                                                 
ggAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCg                                                              <  1:3387311/70‑1 (MQ=255)
 gAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGc                                                            <  1:2388815/71‑1 (MQ=255)
  aaaaGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAAc                                                                <  1:12559/66‑1 (MQ=255)
  aaaaGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCa                                                           <  1:56221/71‑1 (MQ=255)
   aaaGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAg                                                          >  1:1168721/1‑71 (MQ=255)
    aaGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGt                                                         >  1:1178899/1‑71 (MQ=255)
     aGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGt                                                         >  1:2494882/1‑70 (MQ=255)
     aGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGt                                                         >  1:787401/1‑70 (MQ=255)
        gTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAAACGGCAGTGa                                                       <  1:1164141/69‑1 (MQ=255)
          aTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAtttt                                                   >  1:1568794/1‑71 (MQ=255)
          aTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAtttt                                                   >  1:1342106/1‑71 (MQ=255)
              gcgcACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTa                                                <  1:3091679/70‑1 (MQ=255)
              gcgcACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTa                                                <  1:592379/70‑1 (MQ=255)
              gcgcACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTa                                                <  1:3148872/70‑1 (MQ=255)
               cgcACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAag                                              <  1:1389476/71‑1 (MQ=255)
                   cgcgACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGc                                           >  1:791440/1‑70 (MQ=255)
                   cgcgACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGc                                           >  1:367214/1‑70 (MQ=255)
                   cgcgACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGc                                           >  1:1286266/1‑70 (MQ=255)
                         cAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCt                                    >  1:2865870/1‑71 (MQ=255)
                         cAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCt                                    <  1:2633421/71‑1 (MQ=255)
                         cAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCt                                    <  1:3586393/71‑1 (MQ=255)
                          aGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCtt                                   <  1:282085/71‑1 (MQ=255)
                          aGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCt                                    <  1:3143001/70‑1 (MQ=255)
                           gAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCtt                                   >  1:160892/1‑70 (MQ=255)
                           gAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTc                                  <  1:2557750/71‑1 (MQ=255)
                           gAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTc                                  <  1:2119233/71‑1 (MQ=255)
                            aaTGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCtt                                   >  1:2758616/1‑69 (MQ=255)
                            aaTGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTg                                                        <  1:1062383/48‑1 (MQ=255)
                            aaTGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTg                                                        <  1:1448097/48‑1 (MQ=255)
                             aTGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcag                                >  1:603165/1‑71 (MQ=255)
                             aTGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcag                                >  1:3188302/1‑71 (MQ=255)
                               gTAAATTAACCAGTTACTTACTTTACTACAATGaaa                                                                 <  1:985213/36‑1 (MQ=255)
                               gTAAATTAACCAGTTACTTACTTTACTACAATGaaa                                                                 <  1:2070916/36‑1 (MQ=255)
                               gTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTaa                                               >  1:3218697/1‑54 (MQ=255)
                               gTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTaa                                               >  1:2723101/1‑54 (MQ=255)
                               gTAAATTAACCAGTTACTTACTTTACCACAATGaaa                                                                 <  1:2394187/36‑1 (MQ=38)
                                tAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcag                                >  1:2180224/1‑68 (MQ=255)
                                 aaaTTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcagcag                             >  1:308459/1‑70 (MQ=255)
                                    ttAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATg                                       >  1:2258330/1‑57 (MQ=255)
                                       aCCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTc                      <  1:1201656/71‑1 (MQ=255)
                                       aCCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTc                      <  1:1661106/71‑1 (MQ=255)
                                           gTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCt                  <  1:2934103/71‑1 (MQ=255)
                                           gTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCt                  >  1:2697571/1‑71 (MQ=255)
                                           gTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCt                  >  1:1774255/1‑71 (MQ=255)
                                               cttacttTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCtt                                   <  1:182041/50‑1 (MQ=255)
                                               cttacttTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcag                                >  1:3508384/1‑53 (MQ=255)
                                               cttacttTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTg               <  1:972406/70‑1 (MQ=255)
                                                      tactacAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGc          >  1:857212/1‑68 (MQ=255)
                                                        ctacAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGaa     <  1:1399836/71‑1 (MQ=255)
                                                        ctacAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGa      <  1:2482506/70‑1 (MQ=255)
                                                        ctacAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGa      <  1:1581161/70‑1 (MQ=255)
                                                         tacAATGAAACCGGCAGTGATTTTTTAAGAGCCAt                                        <  1:2869760/35‑1 (MQ=255)
                                                           cAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAAttt  >  1:261225/1‑71 (MQ=255)
                                                           cAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAAtt   >  1:3214907/1‑70 (MQ=255)
                                                           cAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAAt    <  1:534119/69‑1 (MQ=255)
                                                           cAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAAt    <  1:2900543/69‑1 (MQ=255)
                                                                tAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcag                                >  1:3613531/1‑36 (MQ=255)
                                                                |                                                                 
GGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTT  >  minE/145466‑145595

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: