Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,236,056 A→G 7.2% L491P (CTG→CCG)  cvrA ← putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,236,0560AG 7.2% 143.0 / 2.6 57L491P (CTG→CCG) cvrAputative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (25/26);  new base G (2/2);  total (28/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 = 3.48e-01

CGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGATGAAGCAACTGGTTATCACGAAACAGTGCCGCA  >  CP009273/1235911‑1236199
                                                                                                                                                 |                                                                                                                                                
cgcgCCTAACAACTGCTGGACAATTTCACCCAGCTTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCTCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAgcg                                                                                                                                               <  2:167187/149‑1 (MQ=255)
 gcgcCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCCGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCGgcgc                                                                                                                                              >  2:682139/1‑149 (MQ=255)
 gcgcCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCGgcgc                                                                                                                                              >  2:682125/1‑149 (MQ=255)
          aacTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACcggcgg                                                                                                                                     >  2:839569/1‑149 (MQ=255)
          aacTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACcggcgg                                                                                                                                     >  1:344677/1‑149 (MQ=255)
                    ccttttCACCCCGCGTTTGCTGCTTAGCACGAGATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTg                                                                                                                           <  1:430004/146‑1 (MQ=255)
                       tttCCCCCCGCGTTTGCTGCTTAGCACAATATTCTCGCCCGTTTTTCAGGCCATATAGCAGCGCAACCTCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAAc                                                                                                                        <  2:161462/149‑1 (MQ=255)
                       tttCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAAc                                                                                                                        <  2:783980/149‑1 (MQ=255)
                             ccAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTa                                                                                                                  >  2:91616/1‑149 (MQ=255)
                              cAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGACACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTAc                                                                                                                 >  1:711128/1‑149 (MQ=255)
                                 cGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGa                                                                                                              <  2:422512/149‑1 (MQ=255)
                                       ctgctTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCg                                                                                                        >  2:730782/1‑149 (MQ=255)
                                              acACGATATTCTTGCCCGTCTCCCAGGCCATATATCAGCGCAACATCGGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAGGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGTACAGTTTACCGAGCGCCGGAAGATc                                                                                                 <  1:40232/148‑1 (MQ=255)
                                                 cGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATcgcg                                                                                              >  1:660395/1‑149 (MQ=255)
                                                       tctcGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCg                                                                                        <  2:111091/149‑1 (MQ=255)
                                                        ctcGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGa                                                                                       <  2:344238/149‑1 (MQ=255)
                                                        ctcGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGa                                                                                       >  1:724186/1‑149 (MQ=255)
                                                        ctcGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGa                                                                                       >  2:652999/1‑149 (MQ=255)
                                                             ccGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAAt                                                                                  >  2:787097/1‑149 (MQ=255)
                                                             ccGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAAt                                                                                  >  2:447309/1‑149 (MQ=255)
                                                             ccGGCTCCCAGGCCAGATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAAt                                                                                  <  2:281099/149‑1 (MQ=255)
                                                              cGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAAtt                                                                                 >  2:58475/1‑149 (MQ=255)
                                                                   tCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTacaca                                                                            <  1:331622/149‑1 (MQ=255)
                                                                      aGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACacacaa                                                                         >  2:184118/1‑149 (MQ=255)
                                                                       ggCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACacacaac                                                                        >  1:770203/1‑149 (MQ=255)
                                                                        gCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACacacaaca                                                                       >  2:833101/1‑149 (MQ=255)
                                                                         ccATATATCGGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGGGCGACGGGCGGGGAGTGGCTGAACGGTTTACCGAGCGCCGGAAGATAGCGTTCCCGGCCGATTACgcacgaccc                                                                      >  1:406031/1‑147 (MQ=255)
                                                                               atCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGcc                                                               >  2:175565/1‑149 (MQ=255)
                                                                                 cAGCGCAACATCATCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAGTAAGCGTTGATCCTGCGCGACCGTCGTCGACTGGCGGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAGTTACACACATCGCATCGCCtt                                                             <  1:282615/149‑1 (MQ=255)
                                                                                    cgcAACATCAGCATATTTGGCGCTGGCTTCGAGCATGAAGTCACCAAAGGAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACAGAGCGCCGGAAGATCGCGTTCGCGACCAATTACACACAACACAACGCCTTcgc                                                          >  1:63592/1‑149 (MQ=255)
                                                                                    cgcAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACGGTTTCCCGAGCGCAGGAAGAACGCGGTCCCGACCGATTACACCAAACCCATCGCCTTcgc                                                          >  1:266297/1‑149 (MQ=255)
                                                                                    cgcAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGCCCAATTACACACAACACCACGCCTTcgc                                                          >  1:713269/1‑149 (MQ=255)
                                                                                    cgcAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTcgc                                                          >  1:860172/1‑149 (MQ=255)
                                                                                       aaCCTCAGCATAGTTGGCGCTGGCTTCGAGAATGAAGTCCCCACAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGCAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCa                                                       <  1:337606/149‑1 (MQ=255)
                                                                                          aTCAGCATATTAGGCGCGGGCTGCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTc                                                    <  2:31743/149‑1 (MQ=255)
                                                                                           tCAGCATATTTGGTGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCg                                                   >  1:212999/1‑149 (MQ=255)
                                                                                             aGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGTAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCggg                                                 <  1:723531/149‑1 (MQ=255)
                                                                                              gCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCGGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGt                                                <  1:682139/149‑1 (MQ=255)
                                                                                              gCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCGGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGt                                                <  1:682125/149‑1 (MQ=255)
                                                                                              gCAGATTGGGCGCTGGCTTCGAGAATGAAGTCACCAAAGACGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGt                                                <  1:120330/149‑1 (MQ=255)
                                                                                                      ggCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCgg                                        >  2:98838/1‑149 (MQ=255)
                                                                                                            gtcgtCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGAt                                  <  1:865495/145‑1 (MQ=255)
                                                                                                            ggCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGAt                                  <  2:776960/149‑1 (MQ=255)
                                                                                                             gCTTCGAGAATGAAGTCACCAAAGAAGCGTTGCTGCCGCGCGACGGGCGGCGACTGGGTGAACAGTTTACCGAGCGCCGGGAGGTCGGGTTCCCGACAAATTAAACACACACAATCACCTTCGCGCAGGCAGGGACTGCCGGTCGGATg                                 >  1:78531/1‑149 (MQ=255)
                                                                                                               ttCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGATGaa                               >  2:81491/1‑149 (MQ=255)
                                                                                                                      aTGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGATGATGCAACTg                        <  1:298062/149‑1 (MQ=255)
                                                                                                                      aTGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGATGAAGCAACTg                        <  2:691971/149‑1 (MQ=255)
                                                                                                                      aTGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGATGAAGCAACTg                        <  2:445913/149‑1 (MQ=255)
                                                                                                                       tGAAGGCACCAACGCCGCGTGGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCTCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGATGAAGCAACTgg                       <  2:112149/149‑1 (MQ=255)
                                                                                                                       ggAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGATGAAGCAACTgg                       <  2:347771/148‑1 (MQ=255)
                                                                                                                             cACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGATGAAGCAACTGGTTATCa                 <  2:711128/149‑1 (MQ=255)
                                                                                                                              aCCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGATGAAGCAACTGGTTATCAc                <  1:91616/149‑1 (MQ=255)
                                                                                                                                cAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGATGAAGCAACTGGTTATCACGa              <  1:455107/149‑1 (MQ=255)
                                                                                                                                          ttGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGATGAAGCAACTGGTTATCACGAAACAGTGCCg    <  2:344677/149‑1 (MQ=255)
                                                                                                                                            gctCCAGCGCGACCGGCTGCGACTGGGTTAACGGTTCACCGAGCGCCGCAACATCGCGTCCCCGACCAATTATACACGACACACCTCCCTTGCGCAGTCGGGTGCTGCCGGTAGGACGAAGCAACTGTTTATCACGGCACAGTGCCCCa  <  1:379605/147‑1 (MQ=255)
                                                                                                                                            gATCCAGCGCGACCGGCGGCGCCTGGCTGAACTGTTTCCCGCGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAGCACATCGCCTTCGCGCAGTGCGGGGCGGCCGGTAGGATGAAGCAACGGGTTATCACGACACCGAGCCGCa  >  1:7353/1‑149 (MQ=255)
                                                                                                                                            gATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGCTGAAGCAACTGGTTATCACGAAACAGTGCCGCa  >  1:663204/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGCGCCTAACAACTGCTGGACAATTTCACCCAGCGTTTGCTGCTTATCACGATATTCTCGCCCGTCTTCCAGGCCATATATCAGCGCAACATCAGCATATTTGGCGCTGGCTTCGAGAATGAAGTCACCAAAGAAGCGTTGATCCAGCGCGACCGGCGGCGACTGGCTGAACAGTTTACCGAGCGCCGGAAGATCGCGTTCCCGACCAATTACACACAACACATCGCCTTCGCGCAGTCGGGTGCTGCCGGTAGGATGAAGCAACTGGTTATCACGAAACAGTGCCGCA  >  CP009273/1235911‑1236199

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