Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,228,356 T→A 30.5% P120P (CCT→CCA yhdZ → predicted amino‑acid transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,228,3560TA30.5% 96.8 / 54.2 82P120P (CCT→CCAyhdZpredicted amino‑acid transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (40/17);  new base A (11/14);  total (51/31)
Fisher's exact test for biased strand distribution p-value = 2.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTCTTCCCTCATCTGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTGAGAATTGCCGAACATGCG  >  minE/2228290‑2228419
                                                                  |                                                               
tctcTTCCCTCATCTGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTa                                                                >  1:2441862/1‑68 (MQ=255)
tctcTTCCCTCATCTGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTa                                                                >  1:3376214/1‑68 (MQ=255)
tctcTTCCCTCATCTGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTa                                                                >  1:485890/1‑68 (MQ=255)
tctcTTCCCTCATCTGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTa                                                                >  1:2775423/1‑68 (MQ=255)
tctcTTCCCTCATCTGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTa                                                                >  1:336075/1‑68 (MQ=255)
tctcTTCCCTCATCTGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTa                                                                >  1:3158301/1‑68 (MQ=255)
 ctctTCCCTCATCTGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTaa                                                               >  1:800257/1‑68 (MQ=255)
              tGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAgg                                                       >  1:682746/1‑63 (MQ=255)
              tGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGc                                                      >  1:2284956/1‑64 (MQ=255)
              tGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGc                                                      >  1:2730974/1‑64 (MQ=255)
              tGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGc                                                      >  1:1386983/1‑64 (MQ=255)
              tGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGc                                                      >  1:2089788/1‑64 (MQ=255)
              tGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGc                                                      >  1:2341564/1‑64 (MQ=255)
              tGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGc                                                      >  1:95341/1‑64 (MQ=255)
              tGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGc                                                      >  1:1936502/1‑64 (MQ=255)
              tGACAGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTaagaaaga                                                         >  1:1688980/1‑61 (MQ=255)
                   gTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTg                                            <  1:1845244/69‑1 (MQ=255)
                   gTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTg                                            <  1:2647626/69‑1 (MQ=255)
                    ttttACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTg                                            <  1:2347472/68‑1 (MQ=255)
                          aGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTg                                      <  1:2294239/68‑1 (MQ=255)
                          aGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTg                                      <  1:861704/68‑1 (MQ=255)
                          aGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTg                                      <  1:951520/68‑1 (MQ=255)
                          aGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTg                                      <  1:2206767/68‑1 (MQ=255)
                           gAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTg                                      <  1:2878610/67‑1 (MQ=255)
                             aCTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:1686684/1‑69 (MQ=255)
                             aCTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:1924668/1‑69 (MQ=255)
                             aCTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:1078743/1‑69 (MQ=255)
                                      tGGCACCGATTTGGGTACGCAAGATGCCAAAGAAAGAGGCTGAAgct                                               >  1:2547114/1‑45 (MQ=255)
                                      tGGCACCGATTTGGGTACGCAAGATGCCAAAGAAAGAGGCTGAAgct                                               >  1:2633750/1‑45 (MQ=255)
                                      tGGCACCGATTTGGGTACGCAAGATGCCAAAGAAAGAGGCTGAAgct                                               >  1:2623277/1‑45 (MQ=255)
                                      tGGCACCGATTTGGGTACGCAAGATGCCAAAGAAAGAGGCTGAAgct                                               >  1:2214022/1‑45 (MQ=255)
                                      tGGCACCGATTTGGGTACGCAAGATGCCAAAGAAAGAGGCTGAAgct                                               >  1:3036557/1‑45 (MQ=255)
                                      tGGCACCGATTTGGGTACGCAAGATGCCAAAGAAAGAGGCTGAAgct                                               >  1:3275081/1‑45 (MQ=255)
                                      tGGCACCGATTTGGGTACGCAAGATGCCAAAGAAAGAGGCTGAAgct                                               >  1:1758294/1‑45 (MQ=255)
                                      tGGCACCGATTTGGGTACGCAAGATGCCAAAGAAAGAGGCTGAAgct                                               >  1:545442/1‑45 (MQ=255)
                                      tGGCACCGATTTGGGTACGCAAGATGCCAAAGAAAGAGGCTGAAgct                                               >  1:1403309/1‑45 (MQ=255)
                                      tGGCACCGATTTGGGTACGCAAGATGCCAAAGAAAGAGGCTGAAgct                                               >  1:1312107/1‑45 (MQ=255)
                                      tGGCACCGATTTGGGTACGCAAGATGCCAAAGAAAGAGGCTGAAgct                                               >  1:1285205/1‑45 (MQ=255)
                                              aTTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGTGGTGCATTACCTAGAGCGGGTgaga                 <  1:1000544/69‑1 (MQ=255)
                                              aTTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTgaga                 <  1:743582/69‑1 (MQ=255)
                                              aTTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTgaga                 <  1:957485/69‑1 (MQ=255)
                                              aTTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTgaga                 <  1:1031483/69‑1 (MQ=255)
                                              aTTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTgaga                 <  1:1807893/69‑1 (MQ=255)
                                              aTTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTgaga                 <  1:1273898/69‑1 (MQ=255)
                                              aTTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTgaga                 <  1:1913922/69‑1 (MQ=255)
                                              aTTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTgaga                 <  1:1670333/69‑1 (MQ=255)
                                               tttGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTgaga                 <  1:3258382/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCGGGTGAGAATTGc            <  1:2041646/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCGGGTGAGAATTGc            <  1:798057/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCGGGTGAGAATTGc            <  1:32578/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCGGGTGAGAATTGc            <  1:754133/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCGGGTGAGAATTGc            <  1:67291/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCGGGTGAGAATTGc            <  1:89354/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCGGGTGAGAATTGc            <  1:1166599/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCGGGTGAGAATTGc            <  1:1049471/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCGGGTGAGAATTGc            <  1:2274344/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCGGGTGAGAATTGc            <  1:2361697/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCGGGTGAGAATTGc            <  1:1036753/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCGGGTGAGAATTGc            <  1:960024/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCGGGTGAGAATTGc            <  1:967604/68‑1 (MQ=255)
                                                    gTACGCAAGATGCCAAAGAAAGAGGCTGAAGCTCTGGCGATGCATTACCTTGAGCCGGTGAGAATTGc            <  1:3081543/68‑1 (MQ=255)
                                                         cAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTGAGAATTGCCGAACa      >  1:3094428/1‑69 (MQ=255)
                                                             aTGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:2535593/1‑37 (MQ=255)
                                                             aTGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:855385/1‑37 (MQ=255)
                                                             aTGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:517665/1‑37 (MQ=255)
                                                             aTGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:482453/1‑37 (MQ=255)
                                                             aTGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:2601813/1‑37 (MQ=255)
                                                             aTGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:253421/1‑37 (MQ=255)
                                                             aTGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:2478674/1‑37 (MQ=255)
                                                             aTGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:2866763/1‑37 (MQ=255)
                                                             aTGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:2224396/1‑37 (MQ=255)
                                                             aTGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:1757231/1‑37 (MQ=255)
                                                             aTGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:1716904/1‑37 (MQ=255)
                                                             aTGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCAt                                   >  1:3265285/1‑36 (MQ=255)
                                                             aTGCATAAGAAAGAGGCTGAAGATCTGGCGGTGCAtt                                  >  1:1458493/1‑37 (MQ=255)
                                                              tGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTGAGAATTGCCGAACATgcg  >  1:1859215/1‑68 (MQ=255)
                                                              tGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTGAGAATTGCCGAACATgcg  >  1:802298/1‑68 (MQ=255)
                                                              tGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTGAGAATTGCCGAACATgcg  >  1:305751/1‑68 (MQ=255)
                                                              tGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTGAGAATTGCCGAACATgcg  >  1:123577/1‑68 (MQ=255)
                                                              tGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTGAGAATTGCCGAACATgcg  >  1:1317290/1‑68 (MQ=255)
                                                              tGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTGAGAATTGCCGAACATgcg  >  1:2798271/1‑68 (MQ=255)
                                                              tGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTGAGAATTGCCGAACATgcg  >  1:3385526/1‑68 (MQ=255)
                                                                  |                                                               
TCTCTTCCCTCATCTGACCGTTTTACAGAACTGTACCCTGGCACCGATTTGGGTACGCAAGATGCCTAAGAAAGAGGCTGAAGATCTGGCGGTGCATTACCTAGAGCGGGTGAGAATTGCCGAACATGCG  >  minE/2228290‑2228419

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: