Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 558,978 G→T 12.8% L75L (CTG→CTT cmr → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported558,9780GT12.8% 112.8 / 4.0 39L75L (CTG→CTTcmrmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (15/19);  new base T (2/3);  total (17/22)
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.34e-01

TTGGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCA  >  Exported/558909‑559036
                                                                     |                                                          
ttgGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGt                                                           >  1:3592948/1‑71 (MQ=255)
   ggTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCg                                                         >  1:4834167/1‑70 (MQ=255)
    gTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGa                                                       <  1:5145216/71‑1 (MQ=255)
        cTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGTATCGt                                                   <  1:2509527/71‑1 (MQ=255)
        cTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCg                                                    <  1:4622537/70‑1 (MQ=255)
        cTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCg                                                    <  1:3370055/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:3187092/71‑1 (MQ=255)
                        gTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCgg                                   <  1:5559143/71‑1 (MQ=255)
                           tCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCg                                           <  1:4474671/60‑1 (MQ=255)
                           tCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCgg                                   <  1:3590787/68‑1 (MQ=255)
                           tCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCgg                                   <  1:2676914/68‑1 (MQ=255)
                                   gcgggATGTTTTTACAATGGCTGCTGGGGCCGCTGTCgg                                                        <  1:5157700/39‑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:5176543/43‑1 (MQ=255)
                                                                 cGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCggg                      <  1:1160401/43‑1 (MQ=255)
                                                                     |                                                          
TTGGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAGTGGTGTGGTTTATCGTCA  >  Exported/558909‑559036

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: