Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,488,439 T→A K116I (AAA→ATA)  emrK ← EmrKY‑TolC multidrug resistance efflux pump, membrane fusion protein component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,488,4390TA93.0% 289.6 / 11.9 100K116I (AAA→ATA) emrKEmrKY‑TolC multidrug resistance efflux pump, membrane fusion protein component
Reads supporting (aligned to +/- strand):  ref base T (5/2);  new base A (68/25);  total (73/27)
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.47e-01

Only 100 of 105 total aligned reads displayed.

ATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTTTATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTAT  >  W3110S.gb/2488381‑2488502
                                                          |                                                               
aTCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      >  1:1480349/1‑70 (MQ=255)
aTCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      >  1:3188256/1‑70 (MQ=255)
aTCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      >  1:2444163/1‑70 (MQ=255)
aTCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      >  1:1649780/1‑70 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGGGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:2142826/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:10079/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:220924/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:3969807/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:2364587/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:468327/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:2048537/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:2032389/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:1930064/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:3502935/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:806948/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:2594124/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:936541/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:869790/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:1236528/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:2718838/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATACGAGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:2793013/1‑69 (MQ=255)
   ttctAAAGATTGTTGATACTGAATAAGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGta                                                    >  1:2119591/1‑69 (MQ=255)
              gttgATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGtt                                       >  1:3421745/1‑71 (MQ=255)
              gttgATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGtt                                       >  1:1834077/1‑71 (MQ=255)
              gttgATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGtt                                       >  1:1219183/1‑71 (MQ=255)
                 gATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGtata                                                               >  1:3981535/1‑44 (MQ=255)
                 gATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGtata                                                               >  1:2310532/1‑44 (MQ=255)
                 gATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGtata                                                               >  1:1283252/1‑44 (MQ=255)
                 gATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGtata                                                               >  1:2415410/1‑44 (MQ=255)
                 gATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGc                                    >  1:3654432/1‑71 (MQ=255)
                 gATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGc                                    >  1:2983094/1‑71 (MQ=255)
                 gATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGc                                    >  1:2779724/1‑71 (MQ=255)
                 gATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGc                                    >  1:3466156/1‑71 (MQ=255)
                 gATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGc                                    >  1:1227718/1‑71 (MQ=255)
                 gATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGc                                    >  1:905767/1‑71 (MQ=255)
                 gATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGc                                    >  1:1507276/1‑71 (MQ=255)
                 gATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGc                                    >  1:1512260/1‑71 (MQ=255)
                 gATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGc                                    >  1:1673039/1‑71 (MQ=255)
                     cTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTAtt                                           >  1:2302383/1‑60 (MQ=255)
                     cTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTAtt                                           >  1:2731669/1‑60 (MQ=255)
                     cTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTAtt                                           >  1:2681311/1‑60 (MQ=255)
                     cTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTAtt                                           >  1:2950957/1‑60 (MQ=255)
                     cTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTAtt                                           >  1:640127/1‑60 (MQ=255)
                     cTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTAtt                                           >  1:1811825/1‑60 (MQ=255)
                     cTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTAtt                                           >  1:2136378/1‑60 (MQ=255)
                     cTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTAtt                                           >  1:2172541/1‑60 (MQ=255)
                          tACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAata                           <  1:1745716/71‑1 (MQ=255)
                          tACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAata                           >  1:3381882/1‑71 (MQ=255)
                          tACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAata                           >  1:2544891/1‑71 (MQ=255)
                          tACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAata                           <  1:606925/71‑1 (MQ=255)
                          tACGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAata                           >  1:507187/1‑71 (MQ=255)
                            cGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATAtt                         >  1:822102/1‑71 (MQ=255)
                            cGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATAtt                         >  1:3721572/1‑71 (MQ=255)
                            cGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATAtt                         >  1:4093641/1‑71 (MQ=255)
                            cGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATAtt                         >  1:3964129/1‑71 (MQ=255)
                            cGTGCTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATAtt                         >  1:1013557/1‑71 (MQ=255)
                                cTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGc                      <  1:2643043/70‑1 (MQ=255)
                                cTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGc                      <  1:2400992/70‑1 (MQ=255)
                                cTGAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGc                      <  1:2641087/70‑1 (MQ=255)
                                  gAAGCGACTTCGGCACTGTATTGTTTATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGc                      >  1:278969/1‑68 (MQ=255)
                                  gAAGCGACTTCGGCACTGTATTGTTTATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGc                      >  1:912854/1‑68 (MQ=255)
                                  gAAGCGACTTCGGCACTGTATTGTTTATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGc                      >  1:2817187/1‑68 (MQ=255)
                                  gAAGCGACTTCGGCACTGTATTGTTTATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGc                      >  1:651430/1‑68 (MQ=255)
                                  gAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCg                                  <  1:3551663/56‑1 (MQ=255)
                                  gAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCg                                  <  1:3851810/56‑1 (MQ=255)
                                  gAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCg                                  <  1:2273431/56‑1 (MQ=255)
                                  gAAGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCg                                  <  1:1595506/56‑1 (MQ=255)
                                   aaGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      <  1:2106225/35‑1 (MQ=255)
                                   aaGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      <  1:441159/35‑1 (MQ=255)
                                   aaGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      <  1:3659232/35‑1 (MQ=255)
                                   aaGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      <  1:2406640/35‑1 (MQ=255)
                                   aaGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      <  1:3624793/35‑1 (MQ=255)
                                   aaGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      <  1:3677868/35‑1 (MQ=255)
                                   aaGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      <  1:3227244/35‑1 (MQ=255)
                                   aaGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      <  1:2974410/35‑1 (MQ=255)
                                   aaGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      <  1:3182596/35‑1 (MQ=255)
                                   aaGCGACTTCGGCACTGTATTGTATATCCTGTAAg                                                      <  1:299304/35‑1 (MQ=255)
                                    aGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGaa                                >  1:2760164/1‑56 (MQ=255)
                                    aGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGaa                                >  1:2330035/1‑56 (MQ=255)
                                    aGCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGaa                                >  1:3171062/1‑56 (MQ=255)
                                     gCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGaa                 >  1:2701/1‑69 (MQ=255)
                                     gCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGaa                 >  1:2281912/1‑69 (MQ=255)
                                     gCGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGaa                 >  1:1633498/1‑69 (MQ=255)
                                      cGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTg                       >  1:2093874/1‑63 (MQ=255)
                                      cGACTTCGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGatta               <  1:2436273/71‑1 (MQ=255)
                                            cGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTAttttt          <  1:2453777/70‑1 (MQ=255)
                                            cGGCACTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTAttttt          <  1:3331741/70‑1 (MQ=255)
                                                aCTGTATTGTTTATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCtt     <  1:1553939/71‑1 (MQ=255)
                                                aCTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCtt     <  1:1148032/71‑1 (MQ=255)
                                                aCTGTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCtt     <  1:415807/71‑1 (MQ=255)
                                                   gTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAg        >  1:50184/1‑65 (MQ=255)
                                                   gTATTGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTAt  >  1:1385031/1‑71 (MQ=255)
                                                      ttGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTa   >  1:2019953/1‑67 (MQ=255)
                                                      ttGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTa   >  1:1988752/1‑67 (MQ=255)
                                                      ttGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTa   >  1:309090/1‑67 (MQ=255)
                                                      ttGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTa   >  1:3969413/1‑67 (MQ=255)
                                                      ttGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTa   >  1:2601940/1‑67 (MQ=255)
                                                      ttGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTa   >  1:1259869/1‑67 (MQ=255)
                                                      ttGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTa   >  1:2468677/1‑67 (MQ=255)
                                                      ttGTATATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTa   >  1:2425356/1‑67 (MQ=255)
                                                          |                                                               
ATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTATTGTTTATCCTGTAAGTATAGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTAT  >  W3110S.gb/2488381‑2488502

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