Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_tpiA 482,230 G→A P725L (CCG→CTG)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_tpiA482,2300GA98.2% 152.6 / NA 55P725L (CCG→CTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (30/24);  minor base C (1/0);  total (31/24)
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

GACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAA  >  NC_000913_3_vch_tpiA/482085‑482365
                                                                                                                                                 |                                                                                                                                       
gacatAAACTTTCTTCACACGACCGCGGTGGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTCATGTCGTTGATATAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCaggg                                                                                                                                      >  2:637527/1‑149 (MQ=255)
   atAACCTTTCTTAACAAGACCGCGGTCGATAAAGTCGGTCACATAGCTCCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCATGGTCCATATCAATCTTAAAATGCAAGGTAt                                                                                                                                   >  2:499393/1‑149 (MQ=255)
   atAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTTGTCTATATCAATCTTCAACTGCAGGGtgt                                                                                                                                   >  2:688732/1‑147 (MQ=255)
    tAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGACACACCCAGCGCCTGCGCTTTTTCCTGGTCCATATCAATCTTAACCTGCAGGGTATc                                                                                                                                  >  2:85661/1‑149 (MQ=255)
    tAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCGGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGTAGGGTTTc                                                                                                                                  >  2:87481/1‑149 (MQ=255)
       aCTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTc                                                                                                                               <  1:634779/149‑1 (MQ=255)
         tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGTTATCAATCTTACACTGCCGGGTATCTTCCa                                                                                                                             >  2:554564/1‑149 (MQ=255)
         tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTACACTGCAGGGTATCTTccc                                                                                                                             >  2:576264/1‑148 (MQ=255)
             ttcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGAcc                                                                                                                         >  2:407173/1‑149 (MQ=255)
                 cgcGAACGCTGTCGATAAAGGCGTTCACAGAGCTGCCGCCCCATGCAGCGCCCAGGGTGGTGGTAATGTCGTGGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAAGTGCAGGGTATCTTCCAGACCGttt                                                                                                                     <  2:18228/147‑1 (MQ=255)
                 cacGCCCGTGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGttt                                                                                                                     <  1:839210/149‑1 (MQ=255)
                    gACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGa                                                                                                                  <  1:324162/149‑1 (MQ=255)
                         cgGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAACCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTAc                                                                                                             >  2:304389/1‑149 (MQ=255)
                         cgGTCGATAAAGTCGTTCACATACCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCGGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTAc                                                                                                             <  1:711003/149‑1 (MQ=255)
                            tCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGTTATCTTCCAGa                                                                                                                           <  1:530560/132‑1 (MQ=255)
                            tCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGCTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGa                                                                                                                           >  2:530560/1‑132 (MQ=255)
                              gATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCATCGCCTGCGCTTTTTCCTGGTCGATAACAATACTAAACACAAGGGTATCTTCCACACCGCTTGGACCTACGCtac                                                                                                        >  2:16021/1‑147 (MQ=255)
                                tAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTc                                                                                                      >  2:128352/1‑149 (MQ=255)
                                  aaGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAACGCCTGCGCTTTTTCCTGGTCGATAGCAATCTTAAACATCAGGGTATCTTAAAGACCGTTTAGACGTACGCTgtcaaa                                                                                                    >  1:177656/1‑144 (MQ=255)
                                       gTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGTTATATACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACTGTTTGGACGTACGCTGGTCAACatat                                                                                               >  2:562232/1‑149 (MQ=255)
                                                     cgccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgct                                                                                 <  1:128352/149‑1 (MQ=255)
                                                               aGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTAGCTGCTTCTgcaa                                                                       <  2:313683/149‑1 (MQ=255)
                                                                gcgcCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaag                                                                      >  1:5289/1‑149 (MQ=255)
                                                                  gcCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaagca                                                                    >  1:522877/1‑149 (MQ=255)
                                                                  gcCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaagca                                                                    >  1:522920/1‑149 (MQ=255)
                                                                   cccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATTTCAATCTTAAACTGCAGGGTATCTTCCAGCCCGTTTGGACGTACGCTGGTCAACAGATCAGGGTGCTTCGCTGCTTCTgcaagcaa                                                                   <  1:41795/149‑1 (MQ=255)
                                                                    ccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGCTATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                  >  1:583378/1‑149 (MQ=255)
                                                                     cAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTCGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACt                                                                 <  1:452393/149‑1 (MQ=255)
                                                                       gagTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTgg                                                               <  1:278467/149‑1 (MQ=255)
                                                                         gtgttgTTAATGTCGGGGATAGAAACCCCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGtt                                                             <  2:14507/149‑1 (MQ=255)
                                                                            gtgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcg                                                          >  2:350494/1‑149 (MQ=255)
                                                                             tgtTAATGTCGGTGATAGAAATACCCAGCGCTTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTGGCTGCTTCTGCAAGCAACTGGTTAcgc                                                         <  2:280436/149‑1 (MQ=255)
                                                                             tgtTAAGGTCGTTGAGAGAGACACCCAGCGCCTGCGCTTTTTCCTGGTCGAGATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                         <  1:553520/149‑1 (MQ=255)
                                                                              gtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGGTATCAAACTTAAAGTGAAGGGTATCTTCCCGACTGTTCGGAAGCCCGCTGGGCATCATATCAGGCTGCTTCTCTGCTTCTTCAACTAACCGGTTAcgca                                                        >  2:469661/1‑148 (MQ=255)
                                                                               ttAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATATTCCAGACCGTTTGGACGTACGCGGGTAAACATATCAGGGTGCTTCGCTGCTTCTCCAACCAACTGGTTGcgcga                                                       >  2:774622/1‑148 (MQ=255)
                                                                                       gTTGATCGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                               <  1:157295/149‑1 (MQ=255)
                                                                                       gTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                               >  2:100938/1‑149 (MQ=255)
                                                                                            tAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtttt                                          >  1:539102/1‑149 (MQ=255)
                                                                                              gAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGTAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAATTtctc                                        <  2:526317/149‑1 (MQ=255)
                                                                                                  cgcccAGCGGCTGCGCTTTTTCCGGGCCGATATAAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGa                                    <  2:284172/147‑1 (MQ=255)
                                                                                                     ccAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCa                                 >  2:201071/1‑149 (MQ=255)
                                                                                                              gcgcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGcc                        <  1:221655/149‑1 (MQ=255)
                                                                                                               cgcTTTTTCCTGGTCGATATCAATATTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTAGTGACCAAGGCCAGCCt                       >  1:627313/1‑149 (MQ=255)
                                                                                                                gcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                      <  1:201071/149‑1 (MQ=255)
                                                                                                                gcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                      >  2:157114/1‑149 (MQ=255)
                                                                                                                gcTTTTTCCTGGTCGATACCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                      <  2:425162/149‑1 (MQ=255)
                                                                                                                  tttttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCTAGCCCGGt                    >  2:450253/1‑149 (MQ=255)
                                                                                                                  tttttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGt                    >  1:453619/1‑149 (MQ=255)
                                                                                                                       ccTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACCCTGTACAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATc               <  1:407173/149‑1 (MQ=255)
                                                                                                                         tGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTACATCac             >  2:165798/1‑148 (MQ=255)
                                                                                                                          ggTCGATATCAATCTTAAACTGCAGGGTAGCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGGCTATCAGc            <  2:566878/149‑1 (MQ=255)
                                                                                                                              gATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTACTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAACACCccaa        >  2:46115/1‑149 (MQ=255)
                                                                                                                               atatCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa       <  2:602144/149‑1 (MQ=255)
                                                                                                                                tatCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGCCGTACGCGGGTCAACATCTCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaag      >  1:120693/1‑149 (MQ=255)
                                                                                                                                   cAATCTTAAACGGCAGGATATCTTCAGGACCGGTTGGTCGTACGCTGGTCTACATAGCAGGGTGCTTAGGTGGTTCTGCAAGCAACGAGGTATGCGTCTGACTCAGTTTTGCATGACCCAGGCCAGCCTGGTCAATCAGCtcaatgtca   <  2:812801/149‑1 (MQ=255)
                                                                                                                                    aaTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaa  >  2:41995/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                       
GACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAA  >  NC_000913_3_vch_tpiA/482085‑482365

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: