Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 51,307 A→G 36.5% T358T (ACT→ACC surA ← peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE51,3070AG36.5% 27.6 / 30.3 52T358T (ACT→ACCsurApeptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base A (10/23);  new base G (6/13);  total (16/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAAATGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCGCCGCC  >  minE/51248‑51369
                                                           |                                                              
aTTTAATGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCACGGc                                                        <  1:346165/68‑1 (MQ=255)
aTTAAATGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGc                                                        <  1:3176606/68‑1 (MQ=255)
aTTAAATGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGc                                                        <  1:2249307/68‑1 (MQ=255)
aTTAAATGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGc                                                        <  1:950494/68‑1 (MQ=255)
aTTAAATGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGc                                                        <  1:1700268/68‑1 (MQ=255)
 ttAAACGCCAGCCTAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGCGt                                                      <  1:2923529/69‑1 (MQ=255)
 ttAAACGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGCGt                                                      <  1:2410362/69‑1 (MQ=255)
 ttAAACGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGCGt                                                      <  1:349530/69‑1 (MQ=255)
 ttAAACGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGCGt                                                      <  1:3341368/69‑1 (MQ=255)
 ttAAACGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGCGt                                                      <  1:2283313/69‑1 (MQ=255)
 ttAAACGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGCGt                                                      <  1:1310079/69‑1 (MQ=255)
        ccAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGa                                                                <  1:2532850/52‑1 (MQ=255)
        ccAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGa                                                                <  1:2465828/52‑1 (MQ=255)
             cGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCATGCGGGTCAGGGCGTCACGGAAGGCCg                                          <  1:150444/69‑1 (MQ=255)
             cGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCg                                          <  1:1749296/69‑1 (MQ=255)
             cGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCg                                          <  1:1088702/69‑1 (MQ=255)
             cGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCg                                          <  1:1108767/69‑1 (MQ=255)
             cGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCg                                          <  1:595157/69‑1 (MQ=255)
             cGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCg                                          <  1:3000698/69‑1 (MQ=255)
             cGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCg                                          <  1:873022/69‑1 (MQ=255)
             cGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCg                                          <  1:2614864/69‑1 (MQ=255)
             cGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCg                                          <  1:135447/69‑1 (MQ=255)
             cGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCg                                          <  1:1983830/69‑1 (MQ=255)
             cGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCg                                          <  1:2086397/69‑1 (MQ=255)
              gattGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCg                                          <  1:3350097/65‑1 (MQ=255)
              gAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCg                                          <  1:74170/68‑1 (MQ=255)
              gAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGc                                                        <  1:2693302/54‑1 (MQ=255)
                        tGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAata                               >  1:1148921/1‑69 (MQ=255)
                                    cTCATTTGACCTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAg                     >  1:3137814/1‑67 (MQ=255)
                                      cATTTGACCTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCa                 <  1:2258386/69‑1 (MQ=255)
                                       aTTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCa                 >  1:2652041/1‑68 (MQ=255)
                                       aTTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAg                >  1:2310804/1‑69 (MQ=255)
                                       aTTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAg                >  1:967366/1‑69 (MQ=255)
                                       aTTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAg                >  1:1104507/1‑69 (MQ=255)
                                       aTTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAg                >  1:1417095/1‑69 (MQ=255)
                                       aTTTGACCTTTGTTCAGGCGGGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAg                >  1:902969/1‑69 (MQ=255)
                                          tGACCTTTGTTCAGGCGAGTCAGGGCGTCACGGAAg                                              <  1:78144/36‑1 (MQ=255)
                                          tGACCTTTGTTCAGGCGAGTCAGGGCGTCACGGAAg                                              <  1:2186769/36‑1 (MQ=255)
                                          tGACCTTTGTTCAGGCGAGTCAGGGCGTCACGGAAg                                              <  1:1024608/36‑1 (MQ=255)
                                             ccTTTGTTCAGGCGAGTCATGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCgag           <  1:3402908/68‑1 (MQ=255)
                                             ccTTTGTTCAGGCGAGTCAGTGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCgag           <  1:947504/68‑1 (MQ=255)
                                             ccTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCgag           <  1:981926/68‑1 (MQ=255)
                                             ccTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCgag           <  1:3044227/68‑1 (MQ=255)
                                             ccTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCgag           <  1:2772039/68‑1 (MQ=255)
                                             ccTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCgag           <  1:2224945/68‑1 (MQ=255)
                                             ccTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCgag           <  1:1213616/68‑1 (MQ=255)
                                                ttGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATcg       >  1:852594/1‑69 (MQ=255)
                                                     cAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCgccgcc  >  1:3388272/1‑69 (MQ=255)
                                                     cAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCgccgcc  >  1:3238553/1‑69 (MQ=255)
                                                     cAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCgccgcc  >  1:653440/1‑69 (MQ=255)
                                                     cAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCgccgcc  >  1:3055937/1‑69 (MQ=255)
                                                     cAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCgccgcc  >  1:2874611/1‑69 (MQ=255)
                                                     cAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCgccgcc  >  1:2814099/1‑69 (MQ=255)
                                                     cAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCgccgcc  >  1:191954/1‑69 (MQ=255)
                                                           |                                                              
ATTAAATGCCAGCCGAATGAAGAGTGAACCGGTGCACTCATTTGACCTTTGTTCAGGCGAGTCAGGGCGTCACGGAAGGCCGGATCGAAAATATCTGGTGTAGCCCAGCCGAGATCGCCGCC  >  minE/51248‑51369

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: