Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,589,664 A→T 57.1% intergenic (+300/‑239) acrD → / → yffB aminoglycoside/multidrug efflux system/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,589,6640AT57.1% 15.0 / 64.0 49intergenic (+300/‑239)acrD/yffBaminoglycoside/multidrug efflux system/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (7/14);  new base T (10/18);  total (17/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

ATGACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAAAATAC  >  W3110S.gb/2589599‑2589728
                                                                 |                                                                
aTGACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTa                                                                <  1:1143362/68‑1 (MQ=255)
aTGACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTa                                                                <  1:1367977/68‑1 (MQ=255)
  gACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTacac                                                             >  1:2436101/1‑69 (MQ=255)
   aCTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGa                                                                  <  1:3643094/63‑1 (MQ=255)
       ccTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCaa                                                        >  1:1174056/1‑69 (MQ=255)
             ggggACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTc                                                  <  1:951118/69‑1 (MQ=255)
             ggggACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTc                                                  <  1:2679047/69‑1 (MQ=255)
             ggggACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTc                                                  <  1:2362663/69‑1 (MQ=255)
             ggggACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTc                                                  <  1:1001004/69‑1 (MQ=255)
                      tCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATAtt                                          >  1:275336/1‑68 (MQ=255)
                      tCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATAtt                                          >  1:2719647/1‑68 (MQ=255)
                        ttATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCAt                                       >  1:30250/1‑69 (MQ=255)
                          aTGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCa                                        <  1:427390/66‑1 (MQ=255)
                            gTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCAt                                       <  1:1543474/65‑1 (MQ=255)
                            gTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCAt                                       <  1:716829/65‑1 (MQ=255)
                            gTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCAt                                       <  1:212051/65‑1 (MQ=255)
                            gTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCAt                                       <  1:1250541/65‑1 (MQ=255)
                             tAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCAt                                       <  1:3185100/64‑1 (MQ=255)
                             tAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCAt                                       <  1:2763507/64‑1 (MQ=255)
                               aaTACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCa                                 <  1:1674256/68‑1 (MQ=255)
                                 tACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAAc                                                       <  1:2675626/44‑1 (MQ=255)
                                 tACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAAc                                                       <  1:2738532/44‑1 (MQ=255)
                                 tACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCt                              >  1:2785617/1‑69 (MQ=255)
                                     aTGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAAt                           <  1:2674145/68‑1 (MQ=255)
                                     aTGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCAAGATATTCATACGGCATCTAAt                           <  1:2368503/68‑1 (MQ=255)
                                     aTGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCAAGATATTCATACGGCATCTAAt                           <  1:2364410/68‑1 (MQ=255)
                                        gTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACt                        <  1:439032/68‑1 (MQ=255)
                                        gTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACt                        <  1:1045467/68‑1 (MQ=255)
                                        gTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACt                        <  1:2978835/68‑1 (MQ=255)
                                        gTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTg                       >  1:802342/1‑69 (MQ=255)
                                        gTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTg                       >  1:2928436/1‑69 (MQ=255)
                                        gTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGTATCTAATACTg                       >  1:1347878/1‑69 (MQ=255)
                                          ccTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGa                      <  1:1705650/68‑1 (MQ=255)
                                          ccTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGa                      <  1:3358342/68‑1 (MQ=255)
                                             gCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGAttt                   <  1:2540586/68‑1 (MQ=255)
                                                 tgatgaTTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTc              <  1:3192589/69‑1 (MQ=255)
                                                   atgatTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTg            >  1:997947/1‑69 (MQ=255)
                                                     gatTTGTATGTGTTACACAGCAACATTTCGAGATATTCa                                        <  1:3032667/39‑1 (MQ=255)
                                                     gatTTGTATGTGTTACACAGCAACATTTCGAGATATTCa                                        <  1:2998729/39‑1 (MQ=255)
                                                     gatTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTc              <  1:3343509/65‑1 (MQ=255)
                                                     gatTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAACTc              <  1:1674980/65‑1 (MQ=255)
                                                        ttGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTaa                            >  1:1622687/1‑48 (MQ=255)
                                                        ttGTATGTGATACACAGCAACATTTCGAGATATTCa                                        >  1:2609622/1‑36 (MQ=255)
                                                        ttGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTaa                            >  1:2708837/1‑48 (MQ=255)
                                                        ttGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTaa                            >  1:2572005/1‑48 (MQ=255)
                                                        ttGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTaa                            >  1:2059262/1‑48 (MQ=255)
                                                        ttGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTaa                            <  1:74619/48‑1 (MQ=255)
                                                          gTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTaaaa     >  1:746965/1‑69 (MQ=255)
                                                          gTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTaaaa     >  1:1102939/1‑69 (MQ=255)
                                                             tgtgTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAAAATac  <  1:3314534/69‑1 (MQ=255)
                                                                 |                                                                
ATGACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAAAATAC  >  W3110S.gb/2589599‑2589728

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: