Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 310,806 G→C 7.9% G689G (GGC→GGG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE310,8060GC 7.9% 183.9 / 5.0 63G689G (GGC→GGGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base C (3/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAA  >  minE/310737‑310876
                                                                     |                                                                      
cATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGc                                                                       <  1:5482814/71‑1 (MQ=255)
cATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGc                                                                       <  1:83119/71‑1 (MQ=255)
  tatCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCa                                                                     <  1:194227/71‑1 (MQ=255)
            cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                                                           >  1:5678351/1‑71 (MQ=255)
            cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                                                           <  1:4248127/71‑1 (MQ=255)
            cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCa                                                            <  1:4569373/70‑1 (MQ=255)
            cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCa                                                            <  1:2166150/70‑1 (MQ=255)
               cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                                          >  1:2705732/1‑69 (MQ=255)
               cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                                          >  1:1785531/1‑69 (MQ=255)
               cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                                          >  1:1054075/1‑69 (MQ=255)
                  tgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                                                     >  1:3442614/1‑71 (MQ=255)
                   gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                                                     >  1:1410835/1‑70 (MQ=255)
                     ttCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                                                   >  1:4169904/1‑70 (MQ=255)
                     ttCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                                                   >  1:2339927/1‑70 (MQ=255)
                      tCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGcc                                                                  >  1:1967700/1‑54 (MQ=255)
                       cTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa                                                 <  1:1220595/70‑1 (MQ=255)
                       cTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa                                                  <  1:817166/69‑1 (MQ=255)
                         ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtc                                              <  1:1626494/70‑1 (MQ=255)
                           aagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcga                                            <  1:999838/71‑3 (MQ=255)
                           aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaa                                            <  1:1283897/71‑1 (MQ=255)
                           aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaa                                            <  1:2247475/71‑1 (MQ=255)
                            agcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcga                                            >  1:1523076/1‑68 (MQ=255)
                             gcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcgaaa                                          >  1:4088067/1‑70 (MQ=255)
                             gcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcgaaa                                          >  1:1711755/1‑70 (MQ=255)
                             gcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag                                          >  1:425769/1‑71 (MQ=255)
                              caaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGc                                         >  1:2708401/1‑71 (MQ=255)
                              caaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGc                                         >  1:4575506/1‑71 (MQ=255)
                               aaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                                        <  1:5341281/71‑1 (MQ=255)
                                 cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                                                                <  1:5583754/45‑1 (MQ=255)
                                 cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                                          >  1:436058/1‑51 (MQ=255)
                                    gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa                                                  <  1:2018286/56‑1 (MQ=255)
                                    gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa                                                  <  1:835662/56‑1 (MQ=255)
                                    gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                                   <  1:4797019/71‑1 (MQ=255)
                                    gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                                   <  1:4700272/71‑1 (MQ=255)
                                      tACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                                   >  1:4778008/1‑69 (MQ=255)
                                            ccTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTAccc                           <  1:577185/71‑1 (MQ=255)
                                            ccTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTAccc                           <  1:3974914/71‑1 (MQ=255)
                                                aGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGtt                       <  1:2198278/71‑1 (MQ=255)
                                                     gTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtc                                                    <  1:3457405/37‑1 (MQ=255)
                                                      tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                                        <  1:383746/48‑1 (MQ=255)
                                                      tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                                        <  1:4875957/48‑1 (MQ=255)
                                                      tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                                        <  1:3020273/48‑1 (MQ=255)
                                                      tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                                        <  1:2744478/48‑1 (MQ=255)
                                                      tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                                   <  1:1198118/53‑1 (MQ=255)
                                                      tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGc                                  <  1:1731220/54‑1 (MQ=255)
                                                       ttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGa                 >  1:2639938/1‑70 (MQ=255)
                                                        tttCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtca                                             >  1:1446223/1‑41 (MQ=255)
                                                        tttCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtca                                             >  1:2049266/1‑41 (MQ=255)
                                                        tttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCa                                 >  1:4645013/1‑53 (MQ=255)
                                                        tttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGt                        >  1:4281216/1‑62 (MQ=255)
                                                         ttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTTCCCAGTTCCACGATTg              >  1:3998606/1‑71 (MQ=255)
                                                         ttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGAtt               <  1:2025623/70‑1 (MQ=255)
                                                          tCGTGACCAAGGCCAGCCCGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTg              <  1:1264835/70‑1 (MQ=255)
                                                            gTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag                                          >  1:4482129/1‑40 (MQ=255)
                                                            gTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtaaaa                                           >  1:3451915/1‑39 (MQ=255)
                                                            gTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCgg           >  1:134305/1‑71 (MQ=255)
                                                            gTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCg            <  1:5109450/70‑1 (MQ=255)
                                                              gACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCggg          >  1:4076046/1‑70 (MQ=255)
                                                              gACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGAggg          >  1:1973690/1‑70 (MQ=255)
                                                               aCCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGc         <  1:3559267/70‑1 (MQ=255)
                                                                ccAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTcc                     >  1:2453408/1‑57 (MQ=255)
                                                                 cAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                                   <  1:76115/42‑1 (MQ=255)
                                                                   aGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGt     >  1:4766973/1‑70 (MQ=255)
                                                                     gCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCa                                 >  1:1194664/1‑40 (MQ=255)
                                                                     gCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTaa  <  1:2540955/71‑1 (MQ=255)
                                                                     gCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTa   >  1:4836064/1‑70 (MQ=255)
                                                                     |                                                                      
CATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAA  >  minE/310737‑310876

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