Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 310,823 G→A 7.4% L684L (CTG→TTG) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE310,8230GA 7.4% 193.9 / 4.2 67L684L (CTG→TTG) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base A (3/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

TGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCG  >  minE/310755‑310891
                                                                    |                                                                    
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)
     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:2247475/71‑1 (MQ=255)
         aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaa                                                           <  1:1283897/71‑1 (MQ=255)
          agcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcga                                                           >  1:1523076/1‑68 (MQ=255)
           gcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcgaaa                                                         >  1:1711755/1‑70 (MQ=255)
           gcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcgaaa                                                         >  1:4088067/1‑70 (MQ=255)
           gcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag                                                         >  1:425769/1‑71 (MQ=255)
            caaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGc                                                        >  1:4575506/1‑71 (MQ=255)
            caaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGc                                                        >  1:2708401/1‑71 (MQ=255)
             aaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                                                       <  1:5341281/71‑1 (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:3974914/71‑1 (MQ=255)
                          ccTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTAccc                                          <  1:577185/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:2744478/48‑1 (MQ=255)
                                    tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                                                       <  1:3020273/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)
                                                      aGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                                                  <  1:3554108/35‑1 (MQ=255)
                                                         cTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTACCACGAtt                              >  1:1809172/1‑52 (MQ=255)
                                                            gTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGaa          <  1:790889/69‑1 (MQ=255)
                                                            gTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACACAGTTCCACGATTGCGGGCAGGTTAAAGGCGaaaa        >  1:1608180/1‑71 (MQ=255)
                                                               aaTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGaaaa        >  1:1431125/1‑68 (MQ=255)
                                                               aaTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTCCCACGATTGCGGGCAGGTTAAAGGCGaaaa        >  1:4461151/1‑68 (MQ=255)
                                                                aTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCAt    <  1:3107834/71‑1 (MQ=255)
                                                                aTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCAt    <  1:3114386/71‑1 (MQ=255)
                                                                 tCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCAt    <  1:5438985/70‑1 (MQ=255)
                                                                 tCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATc   <  1:2605012/71‑1 (MQ=255)
                                                                 tCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATc   >  1:1425093/1‑71 (MQ=255)
                                                                   aGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGc                            >  1:2873737/1‑44 (MQ=255)
                                                                   aGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCggg                         >  1:5109824/1‑47 (MQ=255)
                                                                   aGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCg  <  1:483712/70‑1 (MQ=255)
                                                                    |                                                                    
TGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCG  >  minE/310755‑310891

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: