Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 558,999 T→G 7.8% R82R (CGT→CGG cmr → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported558,9990TG 7.8% 124.6 / 4.8 51R82R (CGT→CGGcmrmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (23/24);  new base G (2/2);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATT  >  Exported/558930‑559068
                                                                     |                                                                     
cgcgTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTc                                                                      >  1:2157784/1‑71 (MQ=255)
  tctATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCg                                                                    >  1:1792202/3‑71 (MQ=255)
  cgTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTcc                                                                     <  1:2524799/70‑1 (MQ=255)
   gTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCgg                                                                   <  1:3187092/71‑1 (MQ=255)
   gTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCgg                                                                   <  1:5559143/71‑1 (MQ=255)
      tCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCgg                                                                   <  1:3590787/68‑1 (MQ=255)
      tCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCgg                                                                   <  1:2676914/68‑1 (MQ=255)
                  gATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAg                                                    <  1:4317998/71‑1 (MQ=255)
                      tttttACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTcc                                                                     <  1:1351681/50‑1 (MQ=255)
                        tttACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTg                                                                 >  1:2286765/1‑52 (MQ=255)
                        tttACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGa                                                                >  1:776816/1‑53 (MQ=255)
                           aCAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCg                                                                    >  1:3574839/1‑46 (MQ=255)
                                 gctgctGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATc                                      <  1:1111768/70‑1 (MQ=255)
                                 gctgctGGGGCCGCAGTCGGATCGTATTGGTCGCCGTCCGGTGATGc                                                             >  1:150117/1‑47 (MQ=255)
                                   tgctgGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAgtggtgtggt                                           >  1:507946/1‑63 (MQ=255)
                                    gctgGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAgtggtgtg                                             >  1:3568688/1‑60 (MQ=255)
                                    gctgGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCa                                  >  1:2452505/1‑71 (MQ=255)
                                     ctgGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTc                                   <  1:4031155/69‑1 (MQ=255)
                                       ggggCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAgtggtgtgg                                            >  1:5454846/1‑58 (MQ=255)
                                            cGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCggg                                                      <  1:5176543/43‑1 (MQ=255)
                                            cGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCggg                                                      <  1:1160401/43‑1 (MQ=255)
                                                 tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAgtg                                                  <  1:5502867/42‑1 (MQ=255)
                                                 tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAgtg                                                  <  1:3973549/42‑1 (MQ=255)
                                                 tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGtt                                          >  1:4896074/1‑50 (MQ=255)
                                                 tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTa                                        >  1:3870143/1‑52 (MQ=255)
                                                 tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATc                                      <  1:5011852/54‑1 (MQ=255)
                                                 tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCaa                      <  1:4420962/70‑1 (MQ=255)
                                                 tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCaa                      <  1:3473048/70‑1 (MQ=255)
                                                 tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCaa                      <  1:1857765/70‑1 (MQ=255)
                                                 tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAat                     <  1:3471494/71‑1 (MQ=255)
                                                 tCGGATAGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTa                                        >  1:4902140/1‑52 (MQ=255)
                                                 tCGGACCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTa                                        >  1:1146864/1‑52 (MQ=255)
                                                   ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTg                              >  1:3027490/1‑60 (MQ=255)
                                                   ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCaa                      <  1:3168184/68‑1 (MQ=255)
                                                   ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAata                    >  1:3764290/1‑70 (MQ=255)
                                                   ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAata                    >  1:2033159/1‑70 (MQ=255)
                                                   ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAata                    >  1:2810463/1‑70 (MQ=255)
                                                   ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAata                    >  1:505052/1‑70 (MQ=255)
                                                   ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAata                    >  1:2849882/1‑70 (MQ=255)
                                                    gATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGttt                                         <  1:1103271/48‑1 (MQ=255)
                                                      tCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATc                                      >  1:3124240/1‑49 (MQ=255)
                                                      tCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGc                >  1:3998887/1‑71 (MQ=255)
                                                             ggTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAat                     <  1:2869052/59‑1 (MQ=255)
                                                             ggTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAat                     <  1:3398795/59‑1 (MQ=255)
                                                             ggTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGgcgc          >  1:5061011/1‑70 (MQ=255)
                                                             ggTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGgcgc          <  1:2907472/70‑1 (MQ=255)
                                                             ggTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGgcgc          <  1:4396756/70‑1 (MQ=255)
                                                                   cGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGCTgg             >  1:2693211/1‑61 (MQ=255)
                                                                    gTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTTGCaa                      >  1:4451761/1‑51 (MQ=255)
                                                                    gTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCaa                      <  1:4211763/51‑1 (MQ=255)
                                                                    gTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACAtt  <  1:1903323/71‑1 (MQ=255)
                                                                     tCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCaa                      <  1:790432/50‑1 (MQ=255)
                                                                     tCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGGCGCaaaa      <  1:4960116/66‑1 (MQ=255)
                                                                     |                                                                     
CGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATT  >  Exported/558930‑559068

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: