Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,581,266 G→T 16.6% A44D (GCC→GAC)  ttgB ← efflux pump membrane protein TtgB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,581,2660GT16.6% 134.1 / 12.9 52A44D (GCC→GAC) ttgBefflux pump membrane protein TtgB
Reads supporting (aligned to +/- strand):  ref base G (19/24);  new base T (4/5);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGCGGTAATGGTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAGAACTTCGACATGGTTTAC  >  NC_002947/1581118‑1581403
                                                                                                                                                    |                                                                                                                                         
ggTGGCGGTAATGGTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACggcgatgg                                                                                                                                           >  4:1943/1‑149 (MQ=255)
    gCGGTAATGGTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACggcgatggcgat                                                                                                                                       <  5:328749/149‑1 (MQ=255)
          aTGGTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGAtggctgg                                                                                                                                 <  1:353024/149‑1 (MQ=255)
            ggTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTggcg                                                                                                                               >  4:347356/1‑149 (MQ=255)
              tCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTggcggc                                                                                                                             <  2:326698/149‑1 (MQ=255)
                aTGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGgcgc                                                                                                                           <  1:19015/149‑1 (MQ=255)
                aTGCTGCCGTCGGAGGTACTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGgcgc                                                                                                                           <  5:75334/149‑1 (MQ=255)
                   ctgcCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGa                                                                                                                         <  1:357147/148‑1 (MQ=255)
                    tgcCGTCGGAGTTACGTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATa                                                                                                                       >  5:215330/1‑149 (MQ=255)
                      ccGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGAtggct                                                                                                                                   <  6:196004/135‑1 (MQ=255)
                      ccGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGAtggct                                                                                                                                   >  5:196004/1‑135 (MQ=255)
                        gTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTgg                                                                                                                   >  7:146664/1‑149 (MQ=255)
                                 aCTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTgg                                                                                                           <  8:146664/148‑1 (MQ=255)
                                   tttCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGtt                                                                                                         >  6:395621/1‑148 (MQ=255)
                                       cgacgaCACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATc                                                                                                     <  3:1943/148‑1 (MQ=255)
                                        gacgacACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGGCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCgg                                                                                                   <  6:206937/149‑1 (MQ=255)
                                        gacgacACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCgg                                                                                                   >  6:53375/1‑149 (MQ=255)
                                              acaTAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCt                                                                                              >  3:145232/1‑148 (MQ=255)
                                              acaTAACGCAGGTTGTCGAAACCGTTGAGCTGCTGCTCGAGCACTTGCACCACGGTGTCCTGCACGTTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATTGCTGGCTGCGCGATACTGGGGTACTGGTTGATCGGCAGCtt                                                                                             <  5:20974/149‑1 (MQ=255)
                                               caTAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTc                                                                                            <  1:44699/149‑1 (MQ=255)
                                               caTAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTc                                                                                            <  3:143073/149‑1 (MQ=255)
                                                aTAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGgcgc                                                                                                                           >  4:194117/1‑117 (MQ=255)
                                                aTAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGAGCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGgcgc                                                                                                                           <  3:194117/117‑1 (MQ=255)
                                                    cGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGAt                                                                                       <  3:347356/149‑1 (MQ=255)
                                                         gTTGTCGATACGGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACa                                                                                  <  5:143145/149‑1 (MQ=255)
                                                                     gTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGAc                                                                        >  3:298588/1‑147 (MQ=255)
                                                                     gTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGAc                                                                        <  4:298588/147‑1 (MQ=255)
                                                                         aGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCAt                                                                  <  6:215330/149‑1 (MQ=255)
                                                                          gctgctgctCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGa                                                                                        >  8:376423/1‑126 (MQ=255)
                                                                          gctgctgctCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGa                                                                                        <  7:376423/126‑1 (MQ=255)
                                                                                    gATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGtt                                                                                                         >  2:377975/1‑99 (MQ=255)
                                                                                    gATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGtt                                                                                                         <  1:377975/99‑1 (MQ=255)
                                                                                          cTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCa                                                 >  4:167872/1‑149 (MQ=255)
                                                                                             caccacGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCAc                                                <  5:53375/147‑1 (MQ=255)
                                                                                                cacGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAgg                                           >  3:266177/1‑149 (MQ=255)
                                                                                                cacGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAgg                                           >  4:127889/1‑149 (MQ=255)
                                                                                                   cgtgtCCTGCACGGTCTGCGCCGACGCGCCCGGGTCGGTCACGGCGATGTCGCTGGCTGGCCGCGCGATACTGGGGTACTGGTTGACCGGCAGCTTCAGGATCGACAAGGCCCCGACCAGCATGATCCCCAAGGCGATCACCCAGGCGa                                        <  8:344021/148‑1 (MQ=255)
                                                                                                     tgtCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCAcc                                                           >  6:167014/1‑128 (MQ=255)
                                                                                                     tgtCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCAcc                                                           <  5:167014/128‑1 (MQ=255)
                                                                                                     tgtCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAg                                      <  3:127889/149‑1 (MQ=255)
                                                                                                      gtCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCTAAGCTATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCg                                         >  8:307278/1‑145 (MQ=255)
                                                                                                       tCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCAccc                                              >  6:237433/1‑139 (MQ=255)
                                                                                                       tCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCAccc                                              <  5:237433/139‑1 (MQ=255)
                                                                                                         cTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCg                                         >  7:156795/1‑142 (MQ=255)
                                                                                                            cACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGc                               >  7:48089/1‑149 (MQ=255)
                                                                                                               ggTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAACGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGcgat                            >  8:211444/1‑149 (MQ=255)
                                                                                                                       ccGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAg                    <  8:48089/149‑1 (MQ=255)
                                                                                                                              gcCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAGAACTTc              >  4:150612/1‑148 (MQ=255)
                                                                                                                               cccGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAGAACTTCGa            <  5:395621/149‑1 (MQ=255)
                                                                                                                                  tggTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAg                                                                                          >  7:266241/2‑68 (MQ=255)
                                                                                                                                  gggTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAg                                                                                          <  8:266241/68‑1 (MQ=255)
                                                                                                                                   ggTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAGAACTTCGACATg        <  3:150612/149‑1 (MQ=255)
                                                                                                                                       ggTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAGAACTTCGACATGGtt     >  8:7017/1‑148 (MQ=255)
                                                                                                                                         tCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAGAACTTCGACATGGTTTAc  <  3:167872/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
GGTGGCGGTAATGGTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACCGTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTTGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAGAACTTCGACATGGTTTAC  >  NC_002947/1581118‑1581403

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