Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 891,052 A→G 42.3% A179A (GCT→GCC sapF ← predicted antimicrobial peptide transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE891,0520AG42.3% 24.6 / 64.3 85A179A (GCT→GCCsapFpredicted antimicrobial peptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (37/12);  new base G (16/20);  total (53/32)
Fisher's exact test for biased strand distribution p-value = 6.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCA  >  minE/890994‑891113
                                                          |                                                             
tttCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCtt                                                     >  1:3081656/1‑69 (MQ=255)
            ccAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCaata                                          >  1:2577758/1‑68 (MQ=255)
            ccAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCaata                                          >  1:2346031/1‑68 (MQ=255)
            ccAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCaata                                          >  1:334942/1‑68 (MQ=255)
            ccAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCaata                                          >  1:3348669/1‑68 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:98440/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:3397396/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:3399435/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:343277/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:2230817/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:386253/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:2053099/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:663614/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:701543/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:706193/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:1625565/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:1549573/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:749208/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:1231002/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:3005637/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:30227/69‑1 (MQ=255)
                cATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:908921/69‑1 (MQ=255)
                 aTCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:2193952/68‑1 (MQ=255)
                 aTCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:978086/68‑1 (MQ=255)
                             aaCTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCaa                         >  1:956579/1‑68 (MQ=255)
                             aaCTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCaa                         >  1:977799/1‑68 (MQ=255)
                             aaCTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAat                        >  1:2140289/1‑69 (MQ=255)
                             aaCTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAat                        >  1:3394177/1‑69 (MQ=255)
                             aaCTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAat                        >  1:1202019/1‑69 (MQ=255)
                             aaCTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAat                        >  1:1521534/1‑69 (MQ=255)
                             aaCTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAat                        >  1:237363/1‑69 (MQ=255)
                             aaCTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAat                        >  1:671984/1‑69 (MQ=255)
                             aaCTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAat                        >  1:1609214/1‑69 (MQ=255)
                             aaCTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAat                        >  1:218167/1‑69 (MQ=255)
                             aaCTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAat                        >  1:314209/1‑69 (MQ=255)
                             aaCTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAat                        >  1:618889/1‑69 (MQ=255)
                               cTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCaat                                           <  1:2402745/48‑1 (MQ=255)
                                 gCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCa                    >  1:2594605/1‑69 (MQ=255)
                                           cGACATATCCAGTGAGGCGAGGGCTTCATCGGCAATAATGAc                                     <  1:2395862/42‑1 (MQ=255)
                                            gACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgcgc          >  1:253653/1‑68 (MQ=255)
                                            gACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgc            >  1:2942419/1‑66 (MQ=255)
                                            gACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgc            >  1:892482/1‑66 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACtttt                                 >  1:2241929/1‑40 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACt                                    >  1:885367/1‑37 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCTCAATATCAgcgcgc              >  1:558388/1‑59 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgc              >  1:677883/1‑59 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgc              >  1:933083/1‑59 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgc              >  1:947880/1‑59 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgc              >  1:2850421/1‑59 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgc              >  1:2323520/1‑59 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgc              >  1:1755893/1‑59 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgc              >  1:1702485/1‑59 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgc              >  1:1504478/1‑59 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgc              >  1:2155296/1‑59 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgc              >  1:2532314/1‑59 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgc              >  1:1419672/1‑59 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgc              >  1:2870782/1‑59 (MQ=255)
                                                 atCCAGTGAGGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgc              >  1:2893934/1‑59 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgcgc          <  1:466663/62‑1 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgcgc          <  1:1160402/62‑1 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgcgc          <  1:1771264/62‑1 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgcgc          <  1:1227202/62‑1 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgcgc          <  1:1865331/62‑1 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgcgc          <  1:525488/62‑1 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgcgc          <  1:2409753/62‑1 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgcgc          <  1:1490556/62‑1 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgcgc          <  1:125478/62‑1 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgcg           >  1:1238085/1‑61 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGccagacc    >  1:2087009/1‑68 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGccaga      >  1:1294348/1‑66 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGccaga      >  1:847691/1‑66 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGccaga      >  1:727475/1‑66 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGccaga      >  1:1997120/1‑66 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGccaga      >  1:2899555/1‑66 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGccaga      >  1:2507511/1‑66 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGccaga      >  1:3291801/1‑66 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGccaga      >  1:3021934/1‑66 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGccaga      >  1:3013749/1‑66 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGcca        >  1:652924/1‑64 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGAGCGCGccaga      >  1:1231965/1‑66 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGAATATCAGCGCGCGCGccagaccca  >  1:1122534/1‑69 (MQ=255)
                                                  tCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGCGCGCAATATCAGCGCGCGCGccaga      >  1:385011/1‑66 (MQ=255)
                                                  tCCAGTCAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAgcgcgcgcgc          >  1:2683733/1‑62 (MQ=255)
                                                    cAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGccagaccca  <  1:2890338/68‑1 (MQ=255)
                                                     aGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGccagaccca  <  1:1300421/67‑1 (MQ=255)
                                                         taGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATc                     <  1:1802110/43‑1 (MQ=255)
                                                          |                                                             
TTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCA  >  minE/890994‑891113

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: