Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,299,010	G→A	7.4%	T409T (ACG→ACA)	pvdQ →	acyl‑homoserine lactone acylase PvdQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,299,010	0	G	A	7.4%	158.2 / 2.3	54	T409T (ACG→ACA)	pvdQ	acyl‑homoserine lactone acylase PvdQ
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base A (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

CGCAGGCGGCCTACAGCGTGCGAGACGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCTGCAGCAACCCGCAGGCGCAAGGGCGCCTGGTGGTGCTGGATGGCTCGCGTAGCGCATGCCAGTGG  >  NC_002947/3298868‑3299157                                                                                                                                               |                                                                                                                                                    cgcAGGCGGCCTACAGCGTGCGAGACGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGc                                                                                                                                                <  8:270657/148‑1 (MQ=255)    aGGCGGCCTACAGCGTGCGAGACGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTg                                                                                                                                            >  3:33053/1‑149 (MQ=255)        ggcCTACAGCGTGCGAGACGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGAt                                                                                                                                         >  6:228625/1‑148 (MQ=255)           cTACGGCGTGCGAGACGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCagg                                                                                                                                     >  5:129315/1‑149 (MQ=255)              cAGCGTGCGAGACGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCaggg                                                                                                                                    <  4:33053/147‑1 (MQ=255)              cAGCGTGCGAGACGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCagg                                                                                                                                     >  5:234805/1‑146 (MQ=255)              cAGCGTGCGAGACGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCagg                                                                                                                                     <  6:234805/146‑1 (MQ=255)                gCGTGCGAGACGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCa                                                                                                                                       >  5:92204/1‑142 (MQ=255)                gCGTGCGAGACGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCa                                                                                                                                       <  6:92204/142‑1 (MQ=255)                   tGCGAGACGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAgggcgg                                                                                                                             >  7:51864/1‑149 (MQ=255)                      gagaCGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGcc                                                                                                                          >  4:126455/1‑149 (MQ=255)                          cGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGt                                                                                                                      <  4:154491/149‑1 (MQ=255)                                ccTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGa                                                                                                                <  5:188279/149‑1 (MQ=255)                                ccTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGa                                                                                                                <  7:3260/149‑1 (MQ=255)                                   cGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGc                                                                                                                                                    <  7:218168/110‑1 (MQ=255)                                   agagaACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCCGCATTGAAGCGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGc                                                                                                                                                    >  8:218168/2‑110 (MQ=255)                                    gagaACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAAc                                                                                                              >  8:206361/1‑147 (MQ=255)                                    gagaACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAAc                                                                                                              <  5:228625/147‑1 (MQ=255)                                    gagaACTCACGGGGGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAAc                                                                                                              <  7:159291/147‑1 (MQ=255)                                       aaCTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTc                                                                                                         >  4:141355/1‑149 (MQ=255)                                            aCGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGaa                                                                                                               <  6:19291/138‑1 (MQ=255)                                            aCGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGaa                                                                                                               >  5:19291/1‑138 (MQ=255)                                            aCGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAg                                                                                                    >  1:258527/1‑149 (MQ=255)                                                   ttGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTa                                                                                             <  3:126455/149‑1 (MQ=255)                                                      cagcagTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGt                                                                                          >  4:364407/1‑149 (MQ=255)                                                         cagTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGa                                                                                       <  8:51864/149‑1 (MQ=255)                                                          agTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGAc                                                                                      <  2:258527/149‑1 (MQ=255)                                                                     aTCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGtt                                                                           >  4:72109/1‑149 (MQ=255)                                                                     aTCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGtt                                                                           <  3:364407/149‑1 (MQ=255)                                                                      ccAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGGCGAGCAGTTGCAAGGTACCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTg                                                                          <  8:58524/148‑1 (MQ=255)                                                                         aCCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGt                                                                                                          <  4:331336/114‑1 (MQ=255)                                                                         aCCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGt                                                                                                          >  3:331336/1‑114 (MQ=255)                                                                         aCCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGAt                                                                       <  3:141355/149‑1 (MQ=255)                                                                                    aGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGTCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCTGCAGCaa                                                            <  7:206361/149‑1 (MQ=255)                                                                                         ggcggcATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCTGCAGCAACCCGc                                                       <  2:347895/149‑1 (MQ=255)                                                                                           cggcATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCTGCAGCAACCCGCAg                                                     >  1:198646/1‑149 (MQ=255)                                                                                                      ggATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCt                                                                   >  1:354786/1‑124 (MQ=255)                                                                                                      ggATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCt                                                                   <  2:354786/124‑1 (MQ=255)                                                                                                                aGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATAc                                                                     >  8:274195/1‑112 (MQ=255)                                                                                                                aGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATAc                                                                     <  7:274195/112‑1 (MQ=255)                                                                                                                  cAGTTGCAAGGTATCCCTTGGGTCAACACACTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTgg                                                                         >  5:234672/1‑106 (MQ=255)                                                                                                                  cAGTTGCAAGGTATCCCTTGGGTCAACACACTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTgg                                                                         <  6:234672/106‑1 (MQ=255)                                                                                                                  cAGTTGCAAGGTATCCCTTGGGTCAACACACTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTgg                                                                         >  5:289495/1‑106 (MQ=255)                                                                                                                  cAGTTGCAAGGTATCCCTTGGGTCAACACACTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTgg                                                                         <  6:289495/106‑1 (MQ=255)                                                                                                                    gTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCTGCAGCAACCCGCAGGCGCAAGGGCGCCTTGTGGTGCTg                             >  4:15457/1‑148 (MQ=255)                                                                                                                        cAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGtt                                                                           <  3:184085/98‑1 (MQ=255)                                                                                                                        cAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGtt                                                                           >  4:184085/1‑98 (MQ=255)                                                                                                                        cAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCTGCAGCAACCCGCAGGCGCAAGGGCGCCTGGTGGTGCTGGATgg                        >  5:299417/1‑149 (MQ=255)                                                                                                                              aTCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTgg                                                                         >  2:257492/1‑94 (MQ=255)                                                                                                                              aTCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTgg                                                                         <  1:257492/94‑1 (MQ=255)                                                                                                                                  cTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCTGCAGCAACCCGCAGGCGCAAGGGCGCCTGGTGGTGCTGGATGGCTCGCGTAgc              >  5:188988/1‑149 (MQ=255)                                                                                                                                         cAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCTGCAGCAACCCGCAGGCGCAAGGGCGCCTGGTGGTGCTGGATGGCTCGCGTAGCGCATGc        >  7:212999/1‑148 (MQ=255)                                                                                                                                           acacGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCTGCAGCAACCCGCAGGCGCAAGGGCGCCTGGTGGTGCTGGATGGCTCGCGTAGCGCATGCCAg     <  8:212999/149‑1 (MQ=255)                                                                                                                                              cGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCTGCAGCAACCCGCAGGCGCAAGGGCGCCTGGTGGTGCTGGATGGCTCGCGTAGCGCATGCCAGTgg  <  2:198646/149‑1 (MQ=255)                                                                                                                                               gCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCTGCAGCAACCCGCAGGCGCAAGGGCGCCTGGTGGTGCTGGATGGCTCGCGTAGCGCATGCCAGTgg  >  6:361720/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                                    CGCAGGCGGCCTACAGCGTGCGAGACGCCAACCTCGAGAACTCACGGGTGTTGCAGCAGTGGTACCAGATCAACCGAGCCGACAGCCTGGCGGCATTGAAGGGATCGGTCGAGCAGTTGCAAGGTATCCCTTGGGTCAACACGCTGGCCGTGGATCAGGGAGGGCGGGCCTTGTACCTGAACCAGTCGGTAGTGCCCTATGTCGACCAGCAACTGTTGGATACCTGCAGCAACCCGCAGGCGCAAGGGCGCCTGGTGGTGCTGGATGGCTCGCGTAGCGCATGCCAGTGG  >  NC_002947/3298868‑3299157

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