Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 558,981 G→C 10.6% S76S (TCG→TCC cmr → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE558,9810GC10.6% 125.5 / 5.8 47S76S (TCG→TCCcmrmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (17/25);  new base C (2/3);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

GGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATA  >  minE/558912‑559050
                                                                     |                                                                     
ggTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCg                                                                       >  1:4834167/1‑70 (MQ=255)
 gTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGa                                                                     <  1:5145216/71‑1 (MQ=255)
     cTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGTATCGt                                                                 <  1:2509527/71‑1 (MQ=255)
     cTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCg                                                                  <  1:3370055/70‑1 (MQ=255)
     cTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCg                                                                  <  1:4622537/70‑1 (MQ=255)
     cTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGt                                                                 <  1:1168012/71‑1 (MQ=255)
           cGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTa                                                                >  1:2171554/1‑66 (MQ=255)
           cGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTg                                                             <  1:4481084/69‑1 (MQ=255)
           cGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTg                                                             <  1:2587475/69‑1 (MQ=255)
           cGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGt                                                           >  1:983826/1‑71 (MQ=255)
           cGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGt                                                           >  1:648688/1‑71 (MQ=255)
                aCCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGcc                                                       <  1:4582602/70‑1 (MQ=255)
                  cgcgTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTc                                                    >  1:2157784/1‑71 (MQ=255)
                    tctATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCg                                                  >  1:1792202/3‑71 (MQ=255)
                    cgTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTcc                                                   <  1:2524799/70‑1 (MQ=255)
                     gTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCgg                                                 <  1:5559143/71‑1 (MQ=255)
                     gTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCgg                                                 <  1:3187092/71‑1 (MQ=255)
                        tCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCg                                                         <  1:4474671/60‑1 (MQ=255)
                        tCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCgg                                                 <  1:2676914/68‑1 (MQ=255)
                        tCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCgg                                                 <  1:3590787/68‑1 (MQ=255)
                                gcgggATGTTTTTACAATGGCTGCTGGGGCCGCTGTCgg                                                                      <  1:5157700/39‑1 (MQ=255)
                                    gATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAg                                  <  1:4317998/71‑1 (MQ=255)
                                       gTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCg                                                         >  1:3151194/1‑45 (MQ=255)
                                       gTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCg                                                         >  1:367359/1‑45 (MQ=255)
                                        tttttACAATGGCTGCTGGGGCCGCTGTCGGATCg                                                                  <  1:2850000/35‑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:1160401/43‑1 (MQ=255)
                                                              cGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCggg                                    <  1:5176543/43‑1 (MQ=255)
                                                                   tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAgtg                                <  1:3973549/42‑1 (MQ=255)
                                                                   tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAgtg                                <  1:5502867/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:2810463/1‑70 (MQ=255)
                                                                     ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAata  >  1:2849882/1‑70 (MQ=255)
                                                                     ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAata  >  1:505052/1‑70 (MQ=255)
                                                                     ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAata  >  1:2033159/1‑70 (MQ=255)
                                                                     ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAata  >  1:3764290/1‑70 (MQ=255)
                                                                     |                                                                     
GGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATA  >  minE/558912‑559050

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: