Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 1,312,115 C→T 25.7% Q96* (CAA→TAA)  aaeA_2 → p‑hydroxybenzoic acid efflux pump subunit AaeA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1431,312,1150CT25.7% 102.8 / 38.7 59Q96* (CAA→TAA) aaeA_2p‑hydroxybenzoic acid efflux pump subunit AaeA
Reads supporting (aligned to +/- strand):  ref base C (22/21);  new base T (5/10);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-02

GGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATGCCGCCATCGCGCGTGGCGAGCTGCATGAAGCCGAGGC  >  ACB143/1311983‑1312245
                                                                                                                                    |                                                                                                                                  
ggCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAAc                                                                                                                                 >  1:214550/1‑136 (MQ=255)
       ttCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCtt                                                                                                                          <  2:69045/136‑1 (MQ=255)
          gccaaCTCAAGGTCCATGACAACAAGCAGGTCAAGGCTGGCGACTTGAGGATGAGAATTGACAGCGAGCGCAACCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACAGCGCTTccg                                                                                                                       >  6:285454/4‑136 (MQ=255)
                  aatgTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGt                                                                                                                <  6:251934/132‑1 (MQ=255)
                    ggTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGaa                                                                                                             <  6:13536/136‑1 (MQ=255)
                         aTGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAagcagc                                                                                                        <  2:161206/136‑1 (MQ=255)
                           gACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCTAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGccc                                                                                                      >  3:8085/1‑136 (MQ=255)
                             cAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCgg                                                                                                    >  8:36000/1‑136 (MQ=255)
                             cAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCgg                                                                                                    >  5:33299/1‑136 (MQ=255)
                             cAACCAGCAGGTCAAGGCTGGCGACGTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGCGCGTGAAGCAGCCCgg                                                                                                    >  5:10450/1‑136 (MQ=255)
                             cAACCAGCAGGTCAAAGATGGCGACTTCTTGATGAGTATTGACCGCGAGCGCATCCAGGCCGCGCTCCATAAGCCCAGCGCAGTGACCGAGACACGCGCCAAGAACCAGCGCTTGCGCGAGCGTGAAGCAGCCCgg                                                                                                    >  7:263143/1‑136 (MQ=255)
                              aaCCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGc                                                                                                   <  5:10573/136‑1 (MQ=255)
                              aaCCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGc                                                                                                   <  2:42552/136‑1 (MQ=255)
                              aaCCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGc                                                                                                   <  7:143111/136‑1 (MQ=255)
                              aaCCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGc                                                                                                   <  8:135949/136‑1 (MQ=255)
                                  agcagGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTAc                                                                                               >  5:220181/1‑136 (MQ=255)
                                   gcagGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACt                                                                                              >  3:148430/1‑136 (MQ=255)
                                   gcagGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACt                                                                                              >  3:9579/1‑136 (MQ=255)
                                   gcagGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCTCGAGCGTGAAGCAGCCCGGCGTACt                                                                                              >  8:171347/1‑136 (MQ=255)
                                   gcagGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACt                                                                                              >  7:168498/1‑136 (MQ=255)
                                   gcagGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACt                                                                                              >  4:252849/1‑136 (MQ=255)
                                   gcagGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACt                                                                                              >  8:112878/1‑136 (MQ=255)
                                     agGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGc                                                                                            <  2:214550/136‑1 (MQ=255)
                                          aaGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTgg                                                                                       <  3:274005/136‑1 (MQ=255)
                                              cTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGc                                                                                    >  3:115960/1‑135 (MQ=255)
                                                    aCTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                                                                             <  6:167701/136‑1 (MQ=255)
                                                    aCTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                                                                             <  4:6651/136‑1 (MQ=255)
                                                    aCTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                                                                             <  4:9579/136‑1 (MQ=255)
                                                    aCTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                                                                             <  3:176411/136‑1 (MQ=255)
                                                    aCTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                                                                             <  8:262031/136‑1 (MQ=255)
                                                    aCTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                                                                             >  3:41126/1‑136 (MQ=255)
                                                                        cgcgAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGGGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCTGAGGCGATCAGTGCCGAGTTGCGTg                                                         <  1:239111/136‑1 (MQ=255)
                                                                        cgcgAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTg                                                         <  3:255401/136‑1 (MQ=255)
                                                                        cgcgAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTg                                                         <  3:252849/136‑1 (MQ=255)
                                                                        cgcgAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTg                                                         <  3:167000/136‑1 (MQ=255)
                                                                        cgcgAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGGGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTg                                                         <  6:287220/136‑1 (MQ=255)
                                                                        cgcgAGAGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGTCCGAGACACGCGCCCAGCAACTGCGCTTGCTCGTGCGTGAAGAAGCCGGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGtt                                                         <  2:96977/136‑2 (MQ=255)
                                                                         gcgAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTg                                                         >  7:248423/1‑135 (MQ=255)
                                                                         gcgAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTg                                                         >  7:248427/1‑135 (MQ=255)
                                                                         gcgAGCGCTGCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTg                                                         >  1:192231/1‑135 (MQ=255)
                                                                           gagaGCTTCCAGGCCGCGCTCGATCAGGCCCGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGCGGCGATCAGTGCCGAGTTGCGTGaa                                                       >  8:91318/1‑135 (MQ=255)
                                                                            aGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGaaaa                                                     >  7:241359/1‑136 (MQ=255)
                                                                             gcgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAAc                                                    <  6:200756/136‑1 (MQ=255)
                                                                             gcgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAAc                                                    <  4:219893/136‑1 (MQ=255)
                                                                             gcgcTCCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAAc                                                    <  2:192231/136‑1 (MQ=255)
                                                                              cgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACg                                                   >  3:157154/1‑136 (MQ=255)
                                                                                              tCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATgacgc                                   >  5:162203/1‑136 (MQ=255)
                                                                                              tCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGACCCGAGTTGCGTGAAAACGCCCAGATCAAAgccgc                                   >  4:89501/1‑136 (MQ=255)
                                                                                              tCGATAAAGACACCGCGATGACATATACACGCGCCCAGCAACTGCGCTTGCGAGATCTTGAAGAAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATgccgc                                   >  2:256340/1‑136 (MQ=255)
                                                                                                gATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGACAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATGCCGCCa                                 >  3:70733/1‑136 (MQ=255)
                                                                                                           cgcgGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATGCCGCCATCGCGCGTGGc                      <  5:285269/136‑1 (MQ=255)
                                                                                                           cgcgGTGACCGAGACACGCGCCCAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATGCCGCCATCGCGCGTGGc                      <  3:282900/136‑1 (MQ=255)
                                                                                                           cgcgGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATGCCGCCATCGCGCGTGGc                      <  6:11106/136‑1 (MQ=255)
                                                                                                           cgcgGGGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGACGAGTTGCGTGAAAACGCCCAGATCAATGCCGCCATCGCGCGTGGc                      <  7:195098/136‑1 (MQ=255)
                                                                                                                           cgcgCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATGCCGCCATCGCGCGTGGCGAGCTGCATGAAGcc       >  3:155214/1‑135 (MQ=255)
                                                                                                                             cgcCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCGCTGGGGCCAGAGGCGATCAGTGCCTAGTTGCGTGAAAACGCCCAGATCAATGCCGCCATCGCGCGTGGCGAGCTGCATGAAGCCGAg    <  5:76992/136‑1 (MQ=255)
                                                                                                                             cgcCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATGCCGCCATCGCGCGTGGCGAGCTGCATGAAGCCGAg    <  4:174969/136‑1 (MQ=255)
                                                                                                                             cgcCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGCCGTACGGCCCTGGGGCCAGAGTCGATCAGTGCCTACTTGCGTGAAAACGCCCAGATCAATGCCGCCCTCGCGCGTGGCGAGCTGCATGACGCCGAg    <  8:705/136‑1 (MQ=255)
                                                                                                                               cccAGTAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATGCCGCCATCGCGCGTGGCGAGCTGCATGAAGCCGAGgc  >  2:231890/1‑136 (MQ=255)
                                                                                                                                    |                                                                                                                                  
GGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATGCCGCCATCGCGCGTGGCGAGCTGCATGAAGCCGAGGC  >  ACB143/1311983‑1312245

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