Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 558,993 T→C 9.8% G80G (GGT→GGC cmr → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported558,9930TC 9.8% 131.1 / 4.7 51G80G (GGT→GGCcmrmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base C (2/3);  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 = 9.97e-01

CGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGGCGC  >  Exported/558923‑559060
                                                                      |                                                                   
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)
                             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: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:3473048/70‑1 (MQ=255)
                                                        tCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCaa              <  1:4420962/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:2033159/1‑70 (MQ=255)
                                                          ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAata            >  1:2849882/1‑70 (MQ=255)
                                                          ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAata            >  1:3764290/1‑70 (MQ=255)
                                                          ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAata            >  1:505052/1‑70 (MQ=255)
                                                          ggATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAata            >  1:2810463/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:4396756/70‑1 (MQ=255)
                                                                    ggTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGgcgc  <  1:2907472/70‑1 (MQ=255)
                                                                      |                                                                   
CGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGGCGC  >  Exported/558923‑559060

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: