Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 145,530 T→A 14.3% 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,5300TA14.3% 169.6 / 19.8 63T233S (ACA→TCA) sfsApredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (21/33);  new base A (4/5);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

GGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCT  >  minE/145466‑145597
                                                                |                                                                   
ggAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCgg                                                               >  1:3803016/1‑71 (MQ=255)
    aaGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGt                                                           >  1:275595/1‑71 (MQ=255)
    aaGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGt                                                           <  1:2727128/71‑1 (MQ=255)
      gCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGa                                                         <  1:4465598/71‑1 (MQ=255)
       cGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAt                                                        <  1:2279102/71‑1 (MQ=255)
        gTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGt                                                                     <  1:4468428/57‑1 (MQ=255)
        gTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGt                                                                     <  1:3849317/57‑1 (MQ=255)
        gTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGt                                                                     <  1:2777874/57‑1 (MQ=255)
        gTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGt                                                                     <  1:1785464/57‑1 (MQ=255)
        gTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGAtt                                                       <  1:2233031/71‑1 (MQ=255)
         tATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAtt                                                       <  1:4639385/70‑1 (MQ=255)
         tATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAtt                                                       <  1:1713358/70‑1 (MQ=255)
          aTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTa                                                                    <  1:1912536/56‑1 (MQ=255)
           tttGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTa                                                                    <  1:1162096/55‑1 (MQ=255)
             tGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAtttttt                                                   <  1:1243227/70‑1 (MQ=255)
              gcgcACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTaa                                                                   >  1:547909/1‑53 (MQ=255)
              gcgcACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAttt                                                      >  1:2068017/1‑66 (MQ=255)
              gcgcACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAttt                                                      >  1:4680531/1‑66 (MQ=255)
                gcACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAttttt                                                    <  1:883206/66‑1 (MQ=255)
                gcACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAttttt                                                    <  1:3003405/66‑1 (MQ=255)
                    gcgACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCa                                           >  1:1444521/1‑71 (MQ=255)
                    gcgACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCa                                           >  1:2966577/1‑71 (MQ=255)
                    gcgACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCa                                           >  1:4574/1‑71 (MQ=255)
                        ccAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGcc                                       >  1:1574409/1‑71 (MQ=255)
                         cAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCt                                      >  1:4472967/1‑71 (MQ=255)
                            aaTGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCgg                                                               <  1:4581925/43‑1 (MQ=255)
                            aaTGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCgg                                                               <  1:1271974/43‑1 (MQ=255)
                               gTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGcc                                            <  1:336693/59‑1 (MQ=255)
                                  aaTTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGa                              <  1:1323642/70‑1 (MQ=255)
                                    ttAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCGGAAAt                           <  1:554266/71‑1 (MQ=255)
                                    ttAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAAt                           <  1:3257191/71‑1 (MQ=255)
                                    ttAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAAt                           <  1:1233054/71‑1 (MQ=255)
                                    ttAACCAGTTACTTACTTTACTACAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAAt                           >  1:1870773/1‑71 (MQ=255)
                                      aaCCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAAtt                          <  1:4585183/70‑1 (MQ=255)
                                       aCCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGa                              <  1:2694799/65‑1 (MQ=255)
                                        ccAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTcc                       >  1:2405386/1‑71 (MQ=255)
                                           gTTACTTACTTTACTACAATGTAACCGTCAGTGAtttttt                                                   <  1:4428680/40‑1 (MQ=255)
                                           gTTACTTACTTTACTACAATGTAACCGGCAGTGAtttttt                                                   <  1:811866/40‑1 (MQ=255)
                                           gTTACTTACTTTACTACAATGTAACCGGCAGTGAttttt                                                    <  1:1380903/39‑1 (MQ=255)
                                           gTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCt                    <  1:3403509/71‑1 (MQ=255)
                                            ttacttacTTTACTACAATGTAACCGGCAGTGAtttttt                                                   <  1:3211955/39‑1 (MQ=255)
                                             tacttactTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCa                                           >  1:3641154/1‑46 (MQ=255)
                                             tacttactTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCa                                           >  1:117935/1‑46 (MQ=255)
                                             tacttactTTACTACAATGTAACCGGCAGTGATTATTTAAGAGcc                                            >  1:3367934/1‑45 (MQ=255)
                                              acttacttTACTACAATGTAACCGGCAGTGAtttttt                                                   <  1:2525582/37‑1 (MQ=255)
                                              acttacttTACTACAATGTAACCGGCAGTGATTTTTTaa                                                 >  1:1559530/1‑39 (MQ=255)
                                              acttacttTACTACAATGTAACCGGCAGTGATTTTTTaa                                                 >  1:1876472/1‑39 (MQ=255)
                                              acttacttTACTACAATGTAACCGGCAGTGATTTTTTaa                                                 >  1:2621057/1‑39 (MQ=255)
                                              acttacttTACTACAATGTAACCGGCAGTGATTTTTTaa                                                 >  1:408648/1‑39 (MQ=255)
                                               cttacttTACTACAATGTAACCGGCAGTTATTTTTTAAGAGCCATGCCTTCAGCAGAAAt                           <  1:2132136/60‑1 (MQ=255)
                                               cttacttTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAAt                           <  1:3828424/60‑1 (MQ=255)
                                                  actttactACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTa               <  1:1254699/69‑1 (MQ=255)
                                                  actttactACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTa               <  1:1200427/69‑1 (MQ=255)
                                                    tttactACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGc            >  1:3411024/1‑70 (MQ=255)
                                                      tactacAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAAt                           <  1:2876018/53‑1 (MQ=255)
                                                      tactacAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCttt                  >  1:238038/1‑62 (MQ=255)
                                                      tactacAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGc            >  1:1343809/1‑68 (MQ=255)
                                                      tactacAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGc            >  1:2116401/1‑68 (MQ=255)
                                                        ctacAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGt                <  1:200597/62‑1 (MQ=255)
                                                        ctacAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGt                <  1:291248/62‑1 (MQ=255)
                                                        ctacAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGt                <  1:82469/62‑1 (MQ=255)
                                                        ctacAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGt                <  1:1698372/62‑1 (MQ=255)
                                                         tacAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAAt      >  1:3547479/1‑71 (MQ=255)
                                                          acAATGAAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTaa              >  1:2147704/1‑62 (MQ=255)
                                                           cAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAAt      <  1:1036391/69‑1 (MQ=255)
                                                           cAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAAt      <  1:870925/69‑1 (MQ=255)
                                                             aTGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCt  <  1:1234141/71‑1 (MQ=255)
                                                                |                                                                   
GGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCT  >  minE/145466‑145597

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: