Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 558,985 C→A 10.4% R78S (CGT→AGT) ‡ cmr → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported558,9850CA10.4% 170.0 / 4.2 56R78S (CGT→AGT) ‡cmrmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (24/26);  new base A (3/3);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGC  >  Exported/558917‑559054
                                                                    |                                                                     
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)
                               gATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAg                                      <  1:4317998/71‑1 (MQ=255)
                                  gTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCg                                                             >  1:367359/1‑45 (MQ=255)
                                  gTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCg                                                             >  1:3151194/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: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:3473048/70‑1 (MQ=255)
                                                              tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCaa        <  1:1857765/70‑1 (MQ=255)
                                                              tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCaa        <  1:4420962/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:2033159/1‑70 (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:3764290/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)
                                                                    |                                                                     
CTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGC  >  Exported/558917‑559054

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: