Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,298,551 A→T 8.5% Q256H (CAA→CAT)  pvdQ → acyl‑homoserine lactone acylase PvdQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,298,5510AT 8.5% 136.6 / 2.6 47Q256H (CAA→CAT) pvdQacyl‑homoserine lactone acylase PvdQ
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base T (2/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CCGGCAGGAGCGGTTCGCCGCGGAGCGCGGCAGCAATGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGACACCTCCAAGCACTTCACGCTCTACCGCCTGCAACTCGA  >  NC_002947/3298413‑3298694
                                                                                                                                          |                                                                                                                                                
ccGGCAGGAGCGGTTCGCCGCGGAGCGCGGCAGCAATGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGAcc                                                                                                                                            >  4:6361/1‑145 (MQ=255)
ccGGCAGGAGCGGTTCGCCGCGGAGCGCGGCAGCAATGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGAcc                                                                                                                                            <  3:6361/145‑1 (MQ=255)
   gCAGGAGCGGTTCGCCGCGGAGCGCGGCAGCAATGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGTGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑g                                                                                                                                    <  4:234316/149‑1 (MQ=255)
   gCAGGAGCGGTTCGCCGCGGAGCGCGGCAGCAATGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑g                                                                                                                                    <  2:244276/149‑1 (MQ=255)
              tCGCCGCGGAGCGCGGCAGCAATGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTtgat                                                                                                                         >  1:205306/1‑149 (MQ=255)
               cgccgcGGAGCGCGGCAGCAATGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCTTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTtgat                                                                                                                         <  2:34919/148‑1 (MQ=255)
                gccgcGGAGCGCGGCAGCAATGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTtgatgt                                                                                                                       <  5:53866/149‑1 (MQ=255)
                gccgcGGAGCGCGGCAGCAATGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTtgatgt                                                                                                                       <  2:230965/149‑1 (MQ=255)
                  cgcgGAGCGCGGCAGCAATGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTtgatgtga                                                                                                                     <  1:370030/149‑1 (MQ=255)
                  cgcgGAGCGCGGCAGCAATGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTtgatgtga                                                                                                                     <  2:205306/149‑1 (MQ=255)
                          gcgGCAGCAATGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGggcgcg                                                                                                             >  6:365756/1‑149 (MQ=255)
                                  aaTGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTAcc                                                                                                     >  4:186609/1‑149 (MQ=255)
                                      cGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGc                                                                                                 <  3:186609/149‑1 (MQ=255)
                                           gCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGcc                                                                                            >  4:220133/1‑149 (MQ=255)
                                              gTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCggt                                                                                         <  3:66338/149‑1 (MQ=255)
                                                    gcgcAGCGTTAAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGTCGCGGCGCTACCGGCCATGCCGGTGGTCaa                                                                                   >  8:340161/1‑149 (MQ=255)
                                                               aGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTc                                                                         <  1:247261/148‑1 (MQ=255)
                                                               aGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTc                                                                         >  3:7150/1‑148 (MQ=255)
                                                                 ccAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGAGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACa                                                                                 >  8:158564/1‑138 (MQ=255)
                                                                 ccAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGAGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACa                                                                                 <  7:158564/138‑1 (MQ=255)
                                                                                tGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTg                                                        >  6:150901/1‑148 (MQ=255)
                                                                                      tGGCCAAGCGGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAAGTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGGCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGCCGcac                                                 <  7:340161/149‑1 (MQ=255)
                                                                                      tGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGca                                                  >  5:103905/1‑148 (MQ=255)
                                                                                          cAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCggtg                                                                                        <  2:136865/106‑1 (MQ=255)
                                                                                          cAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCggtg                                                                                        >  1:136865/1‑106 (MQ=255)
                                                                                          cAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAAc                                                                      >  5:119895/1‑124 (MQ=255)
                                                                                          cAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAAc                                                                      <  6:119895/124‑1 (MQ=255)
                                                                                           aaCCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCATCTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGt                                            >  2:299340/1‑149 (MQ=255)
                                                                                             cccGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCt                                                                                               >  8:206820/1‑96 (MQ=255)
                                                                                             cccGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCt                                                                                               <  7:206820/96‑1 (MQ=255)
                                                                                              ccGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGa                                         <  3:220133/149‑1 (MQ=255)
                                                                                                gCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGaca                                       >  4:95906/1‑149 (MQ=255)
                                                                                                    ttCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGACACCTc                                   >  3:48517/1‑149 (MQ=255)
                                                                                                      cccGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGACACCTcc                                  <  5:365756/148‑1 (MQ=255)
                                                                                                        cGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGACACCTCCAAg                               >  6:314562/1‑149 (MQ=255)
                                                                                                            gATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCGAGCTTGATGTGATGGGCGCGCCGCTACCTGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTTGACGCACACCGTCGACACCTCCAAGCACt                           >  5:88160/1‑149 (MQ=255)
                                                                                                                ggcggcATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGACACCTCCAAGCACTTCAc                       >  6:131614/1‑149 (MQ=255)
                                                                                                                   ggcATGCGCTTCTACCAGATGCAACTGACCATTCCCTGGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGACACCTCCAAGCACTTCACGc                     >  1:60600/1‑149 (MQ=255)
                                                                                                                        gcgcTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCggtg                                                                                        <  7:320017/76‑1 (MQ=255)
                                                                                                                        gcgcTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCggtg                                                                                        >  8:320017/1‑76 (MQ=255)
                                                                                                                        cccgatCTACCAGATGCATCTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGACACCTCCAAGCACTTCACGCTCTAc                <  3:11294/143‑1 (MQ=255)
                                                                                                                            ttCTACCAGATGCAACTGACCATTCCC‑TGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGAAGCACACCGGCGACACCTCCAAGCACTTAACGCCATACAGCCt           >  4:228371/1‑149 (MQ=255)
                                                                                                                            ttCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGACACCTCCAAGCACTTCACGctc                   >  1:334592/1‑141 (MQ=255)
                                                                                                                            ttCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGACACCTCCAAGCACTTCACGctc                   <  2:334592/141‑1 (MQ=255)
                                                                                                                                 ccAGATGCATCTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGACACCTCCAAGCACTTCACGCTCTAc                >  4:11294/1‑139 (MQ=255)
                                                                                                                                   aGATGCATCTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGACACCTCCAAGCACTTCACGCTCTACCGCCTGCAACTc    <  1:299340/149‑1 (MQ=255)
                                                                                                                                     aTGCAACTTACCATTCCC‑GTCCAGCTTTATGTCATGGGCGAGGCGCTACCTGGCCTGCCGGGGGTCAACATCTGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGACACCTCCAAGCACTTCACGCTCTACCGCCTGCAACTCGa  <  6:365840/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                
CCGGCAGGAGCGGTTCGCCGCGGAGCGCGGCAGCAATGCGGTGGCAGTAGGCGCGCAGCGTTCAGCCAATGGCCGCGGCCTGTTGCTGGCCAACCCGCACTTCCCGTGGATGGGCGGCATGCGCTTCTACCAGATGCAACTGACCATTCCC‑GGCCAGCTTGATGTGATGGGCGCGGCGCTACCGGGCCTGCCGGTGGTCAACATCGGTTTCAACCAGCACCTGGCGTGGACGCACACCGTCGACACCTCCAAGCACTTCACGCTCTACCGCCTGCAACTCGA  >  NC_002947/3298413‑3298694

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