Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 51,837 G→A 50.0% P182S (CCG→TCG)  surA ← peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE51,8370GA50.0% ‑4.7 / 54.4 38P182S (CCG→TCG) surApeptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base G (9/10);  new base A (12/7);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

AATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGC  >  minE/51775‑51900
                                                              |                                                               
aaTGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGa                                                         <  1:377081/71‑1 (MQ=255)
  tGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGAt                                                              >  1:3536380/1‑64 (MQ=255)
  tGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGa                                                         <  1:3884953/69‑1 (MQ=255)
   ggCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGAtg                                                       >  1:191272/1‑70 (MQ=255)
   ggCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGAGATCAGGAtgt                                                      >  1:302153/1‑71 (MQ=255)
             tGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGAGATCAGGATGTGGCTCAGGtt                                            >  1:1868770/1‑71 (MQ=255)
              ggCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGtt                                            >  1:2621379/1‑70 (MQ=255)
               gCTTTCCGCTTCGTTCACCTGAGCAGAGGTCGGGTTTTCCGGCAGCGAGATCAGGa                                                         <  1:3639870/56‑1 (MQ=255)
                   tCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCt                                       <  1:4046938/70‑1 (MQ=255)
                      gCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTc                                      >  1:4083605/1‑68 (MQ=255)
                      gCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAg                                    >  1:3325109/1‑70 (MQ=255)
                      gCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGAGATCAg                                                           <  1:100089/47‑1 (MQ=255)
                          cGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTg                               >  1:2714662/1‑71 (MQ=255)
                           gTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTc                                                 <  1:211713/52‑1 (MQ=255)
                              cACCTGATCAGAGGTCGGGTTTTCCGGCAGCGAGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgt                           >  1:2792661/1‑71 (MQ=255)
                                  tGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtt                       >  1:2732729/1‑71 (MQ=255)
                                      cAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTgg                   >  1:3664184/1‑71 (MQ=255)
                                       agagGTCGGGTTTTCCGGCAGCGAGATCAGGATGTGGCTCAGGTTCa                                          <  1:1317068/47‑1 (MQ=255)
                                       agagGTCGGGTTTTCCGGCAGCGAGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGt                  >  1:530868/1‑71 (MQ=255)
                                         agGTCGGGTTTTACGGCAGCGAGATCAGGATGTGGCTCAGGTTCAGc                                        >  1:2500684/1‑47 (MQ=255)
                                          ggTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTAc               <  1:182681/71‑1 (MQ=255)
                                          ggTCGGGTTTTCCGGCAGCGAGATCAGGATGTGGCTCAg                                               >  1:2337330/1‑39 (MQ=255)
                                           gTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtttt                     <  1:1769662/64‑1 (MQ=255)
                                           gTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtttt                     <  1:3784662/64‑1 (MQ=255)
                                           gTCGGGTTTTCCGGCAGCGAGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgtcg                         <  1:2826156/60‑1 (MQ=255)
                                            tCGGGTTTTCCGGCAGCGAGATCAGGATGTGGCTCAGGTTCa                                          <  1:1885946/42‑1 (MQ=255)
                                             cGGGTTTTCCGGCAGCGAGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCa            >  1:3629740/1‑71 (MQ=255)
                                              gggTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCAc           >  1:1619224/1‑71 (MQ=255)
                                              gggTTTTCCGGCAGCGAGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCAc           >  1:2528306/1‑71 (MQ=255)
                                                gTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCt                                       <  1:2047272/41‑1 (MQ=255)
                                                gTTTTCCGGCAGCGAGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACct         <  1:1327769/71‑1 (MQ=255)
                                                 ttttCCGGCAGCGAGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtttt                     >  1:2757001/1‑58 (MQ=255)
                                                 ttttCCGGCAGCGAGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGttt                      >  1:2997923/1‑57 (MQ=255)
                                                   ttCCGGCAGCGAGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACctgct      >  1:447974/1‑71 (MQ=255)
                                                       ggCATCGAGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTgcgc  <  1:2513238/71‑1 (MQ=255)
                                                        gCAGCGAGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTgcgc  >  1:2515398/1‑70 (MQ=255)
                                                           gCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtttt                     <  1:323768/48‑1 (MQ=255)
                                                           gCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtttt                     <  1:1228423/48‑1 (MQ=255)
                                                              |                                                               
AATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGC  >  minE/51775‑51900

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: