Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 506,394 C→T 7.0% A494V (GCG→GTG)  ushA → 5'‑nucleotidase/UDP‑sugar hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913506,3940CT 7.0% 142.8 / 2.5 57A494V (GCG→GTG) ushA5'‑nucleotidase/UDP‑sugar hydrolase
Reads supporting (aligned to +/- strand):  ref base C (29/24);  new base T (2/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

CTGACCGCCGTCGCGCAGATGAAGCCAGATTCAGGTGCCTACCCGCAATTTGCCAACGTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGAAAAGCTCGCCGCTGGATGTGAGTGTTTAT  >  NC_000913/506252‑506536
                                                                                                                                              |                                                                                                                                              
cTGACCGCCGTCGCGCAGATGAAGCCAGATTCAGGTGCCTACCCGCAATTTGCCAACGTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAACGGCGAACCGGTCGATCCGTCGAAACCTTACCGTCTGGc                                                                                                                                                >  2:573980/1‑143 (MQ=255)
    ccgccgTCGCGCAGATGAAGCCGGAGTCAGGGGCCTACCCGCGATGTGCCAACGTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACa                                                                                                                                            <  1:143095/143‑1 (MQ=255)
     cgccgTCGCGCAGATGAAGCCAGATTCAGGTGCCTACCCGCAATTTGCCAACGTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACAt                                                                                                                                           >  1:376353/1‑143 (MQ=255)
             cgcAGATGAAGCCAGATTCAGGTGCCTACCCGCAATTTGCCAACGTTAGCTCTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACtcc                                                                                                                                   >  2:588615/1‑141 (MQ=255)
                aGATGAAGCCAGATTCAGGTGCCTACCCGCAATTTGCCAACGTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGACCCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAAt                                                                                                                                >  1:877140/1‑143 (MQ=255)
                          aGATTCAGGTGCCTACCCGCAATTTGCCAACGTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACcggcg                                                                                                                      >  1:884092/1‑143 (MQ=255)
                                      cTACCCGCAATTTGCCAACGTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGTACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATc                                                                                                          >  2:994874/1‑143 (MQ=255)
                                        aCCCGCAATTTGCCAACGTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACACTACACTTCAATGCCACCGGCGGTGAAGGATATCcg                                                                                                        >  1:34101/1‑143 (MQ=255)
                                           cGCAATTTGCCAACGTTAGCTTTGTGGCGAAAGACGGTAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCgcgc                                                                                                     >  1:303912/1‑143 (MQ=255)
                                           cGCAATTTGCCAACGTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCgcgc                                                                                                     >  2:565087/1‑143 (MQ=255)
                                           cGCAATTTGCCAACGTTAGCTGTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGCCGGTGCTGGGTATgcgcgc                                                                                                     >  2:520635/1‑143 (MQ=255)
                                           cGCAAATTGCCAACGTTAGCATAGTGGCGAAAGACGGCAAACTCAACGACCTTAAAATCAAAGGCGAACCGGGCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCgcgc                                                                                                     <  1:25666/143‑1 (MQ=255)
                                               aTTTGCCAACGTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTg                                                                                                 <  1:224868/143‑1 (MQ=255)
                                                tttGCCAACGTTAGCTGTGTGGCGAAAGGCGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGa                                                                                                <  1:219150/143‑1 (MQ=255)
                                                 gtGCCAACGTGAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGAt                                                                                               <  1:194768/142‑1 (MQ=255)
                                                   gcaaaCGTGAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAATTCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATaa                                                                                             <  1:518135/140‑1 (MQ=255)
                                                   gCCAACGTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATaa                                                                                             >  2:248290/1‑143 (MQ=255)
                                                       aCGTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACaaa                                                                                         >  1:367198/1‑143 (MQ=255)
                                                         gTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGTGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAAcc                                                                                       >  1:830941/1‑143 (MQ=255)
                                                         gTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATACGGCGAAAACTTACCGTATGGTGACATTAAACTTCAATGCCACCGGCGCTGATGCATATCCGCCCCTTCATAACAAccc                                                                                       >  1:830926/1‑143 (MQ=255)
                                                          ttAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCCCCGGCGGTGATGGATATCCGCGCCTTGCTAACAAACCg                                                                                      >  1:634680/1‑143 (MQ=255)
                                                          ttAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGCTATCCGCGCATTGATAAAAAACCg                                                                                      >  2:672663/1‑143 (MQ=255)
                                                            aGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTCAAGACAATGGTGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCggg                                                                                    <  1:336406/143‑1 (MQ=255)
                                                                  gtgGCGACAGACGGCAAGCTGAACGACCTTGACATCAGAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAATCTTCAATGCCACCGGCGGTTATGGATATCCGCGCCTTGATAACAAACCGGGCTAtgt                                                                              <  2:777782/143‑1 (MQ=255)
                                                                  gtgGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTAtgt                                                                              <  1:320479/143‑1 (MQ=255)
                                                                    ggCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGa                                                                            >  2:173958/1‑143 (MQ=255)
                                                                                cAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATAACGGCATTa                                                                >  2:80224/1‑143 (MQ=255)
                                                                                   gcTGAAGGACCGTAGAATCAAAGGCGAACCGGTCGATGCGGGGAAAACTTACCGTATGGCGACATTAAACTTGAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTg                                                             <  2:1042762/142‑1 (MQ=255)
                                                                                    cTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGa                                                            >  1:167548/1‑143 (MQ=255)
                                                                                       agcGACCTTAAAATCAAGGGCGAACCGGTCGATCCGGCGAAAACTTAGCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGc                                                         <  2:877140/141‑1 (MQ=255)
                                                                                         cGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCCATGTGAATACCGGCTTTATTGACGCCg                                                       >  2:230515/1‑143 (MQ=255)
                                                                                          gACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAGCCGGGCTATGTGAATACCGGCTTTATTGATGCCGa                                                      >  2:847453/1‑143 (MQ=255)
                                                                                            tcTTATAATCAAAGGCGAACCGGTCGTTCCGGCGAAAACTTACCTTATGGTGACATTAAACTTCAATGACACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAg                                                    <  2:830926/142‑1 (MQ=255)
                                                                                            ccTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGTGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAg                                                    <  2:830941/143‑1 (MQ=255)
                                                                                             cTTAAAATCAAGGGCGCACCTGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGt                                                   <  2:140871/143‑1 (MQ=255)
                                                                                              ttAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTAGGTGAATACCGGCTTTATTGATGCCGACGTg                                                  >  1:35636/1‑143 (MQ=255)
                                                                                                 aaaCCAAAGGCGCACCGGTCGATCCGGCGAGAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTg                                               <  2:726717/143‑1 (MQ=255)
                                                                                                     cagaGGCGAACCGGTCGAGCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAg                                           <  2:315197/140‑1 (MQ=255)
                                                                                                     cAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAg                                           >  1:882404/1‑143 (MQ=255)
                                                                                                      aaaGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGTCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGGAGTGGTGaaatg                                          >  2:987230/1‑141 (MQ=255)
                                                                                                            gAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCCTTATTGATGCCGAAGTGCTGAAAGCGTatat                                    >  2:815898/1‑143 (MQ=255)
                                                                                                             aaCCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATc                                   <  2:367198/143‑1 (MQ=255)
                                                                                                                 ggTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGa                               <  1:432669/143‑1 (MQ=255)
                                                                                                                    caaaccGACGAAACCTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGATATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGaaaa                            <  1:87851/139‑1 (MQ=255)
                                                                                                                      aTCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCGACGCCTTGCTAAAAAACCGGGCTATATAAATACAGGCTATAATACTGACCAAGCGCGGAAAGCGTACATCCAGAAAAGc                          >  2:266285/1‑143 (MQ=255)
                                                                                                                      aTCCGGCGAAAACTTACCGTATGGCGACACTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGaaaacc                          >  2:304057/1‑141 (MQ=255)
                                                                                                                       tgcGGCGAAAGCGTACCGTAGGGCGGCAGTAAGCTGCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGAAAAGCt                         <  2:462211/141‑1 (MQ=255)
                                                                                                                         ccgCGAAAACGTACCGTATGGCGACATTAATCTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGAAAAGCTcg                       <  2:651235/141‑1 (MQ=255)
                                                                                                                         cGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGAAAAGCTcg                       <  1:994874/143‑1 (MQ=255)
                                                                                                                           gCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGAAAAGCTcgcc                     >  1:279815/1‑143 (MQ=255)
                                                                                                                                  cTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGAAAAGCTCGCCGCTGGAt              <  2:167548/143‑1 (MQ=255)
                                                                                                                                   ttACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGCAGTGCTGAAAGCGTATATCCAGAAAAGCTCGCCGCCGGagg             >  2:152710/1‑141 (MQ=255)
                                                                                                                                      ccGTATGGCTACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGAAAAGCTCGCCGCTGGATGTGa          <  2:382768/143‑1 (MQ=255)
                                                                                                                                       cGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCAGAAAGCGGATAACCAGAAAAGCTCGCCGCTGGATGTGAg         >  2:963419/1‑143 (MQ=255)
                                                                                                                                       cGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCGATGTGAATACCGGCGTTATTGAGGCCGAAGTGCTGAAAGCGGATAACAGGAAAAGCTGGCCGCCGGATGTGAg         >  2:964061/1‑143 (MQ=255)
                                                                                                                                           tGGCGACATCAGTCTTCAATTCCACCGGCGGTGATGGATAGCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGAAAAGCTCGCCGCTGGATGTGAGTGtt     <  2:165268/143‑1 (MQ=255)
                                                                                                                                             gCGTCATTAGACTTTGATGCCACCGCCCGTGATGTCTATTCTCGTCTTCAGAACAAACCGGGCTATGTGAAGACCGGCTGTATTCATGGCGATGTGCTGAAAGCGTATATCCAGAAAAGCTCGCCGCGGGATGTGAGTGTTTa   <  2:310510/143‑1 (MQ=255)
                                                                                                                                             gCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGAAAAGCTCGCCGCTGGATGTGAGTGTTTa   >  1:299533/1‑143 (MQ=255)
                                                                                                                                              cGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGACGTGCTGAAAGCGTGTATACAGGGATACGCGCCGCTGGATGTGAGTGTTTAt  >  1:937779/1‑143 (MQ=255)
                                                                                                                                              cGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGAAAAGCTCGCCGCTGGATGTGAGTGGTTAt  >  1:715524/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                              
CTGACCGCCGTCGCGCAGATGAAGCCAGATTCAGGTGCCTACCCGCAATTTGCCAACGTTAGCTTTGTGGCGAAAGACGGCAAACTGAACGACCTTAAAATCAAAGGCGAACCGGTCGATCCGGCGAAAACTTACCGTATGGCGACATTAAACTTCAATGCCACCGGCGGTGATGGATATCCGCGCCTTGATAACAAACCGGGCTATGTGAATACCGGCTTTATTGATGCCGAAGTGCTGAAAGCGTATATCCAGAAAAGCTCGCCGCTGGATGTGAGTGTTTAT  >  NC_000913/506252‑506536

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: ‑