Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,468,574	T→A	7.8%	A344A (GCA→GCT)	sasA_16 ←	Adaptive‑response sensory‑kinase SasA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,468,574	0	T	A	7.8%	159.6 / 4.0	51	A344A (GCA→GCT)	sasA_16	Adaptive‑response sensory‑kinase SasA
Reads supporting (aligned to +/- strand):  ref base T (23/24);  new base A (2/2);  total (25/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 = 9.20e-01

GGTTTCCTAGCTGTCAACGGTAGGCGCTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGCTTGAACAGATAGGGCATCACATCGCTGGGAATCGGCTCTCCGGAATTGT  >  ACB122/4468440‑4468701                                                                                                                                       |                                                                                                                                gggttCCTAGCTGTCAACGGTAGGCGCTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCcgat                                                                                                                                 <  5:13251/132‑1 (MQ=255) ggTTTCCTAGCTGTCAACGGTAGGCGCTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCACCGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCcgat                                                                                                                                 <  7:344710/135‑1 (MQ=255) ggTTTCCTAGCTGTCAACGGTAGGCGCTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCcgat                                                                                                                                 <  3:49679/135‑1 (MQ=255) ggTTTCCTAGCTGGCAACGGTAGGCGCGTCGAGGCGGACTGAAGGCATCTGGACAGAATTAGGTGAAAGGTCGTCCCGCCTTCAACGCCAGACTCGTGCTCTATCCTGCCAGGATTTCTGGCGACTATTTCcgat                                                                                                                                 <  5:79623/135‑1 (MQ=255)       cTAGCTGTCAACGGTAGGCGCTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCcgatgcgat                                                                                                                            <  6:227321/134‑1 (MQ=255)       cTAGCTGTCAACGGTAGGCGCTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCcgatgcgat                                                                                                                            >  4:173045/1‑134 (MQ=255)       cTAGCTGTCAACGGTAGGCGCTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCcgatgcgat                                                                                                                            >  7:237874/1‑134 (MQ=255)        tAGCTGTCAACGGTAGGCGCTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCcgatgcgatg                                                                                                                           <  3:295140/134‑1 (MQ=255)           cTGTCAACGGTAGGCGCTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATa                                                                                                                       <  3:72270/135‑1 (MQ=255)               cAACGGTAGGCGCTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGCCCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGccc                                                                                                                   >  5:255473/1‑135 (MQ=255)               cAACGGTAGGCGCTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGc                                                                                                                     <  6:53554/133‑1 (MQ=255)               cAACGGTAGGCGCTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGc                                                                                                                     <  7:158077/133‑1 (MQ=255)                       ggCGCTTCCAGTCGGATTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTa                                                                                                           >  6:336933/1‑135 (MQ=255)                          gcTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAAcc                                                                                                        <  5:231168/135‑1 (MQ=255)                                  tCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCt                                                                                                <  3:173045/135‑1 (MQ=255)                                         gAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCt                                                                                         >  6:90012/1‑135 (MQ=255)                                              cATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTttct                                                                                    <  8:342584/135‑1 (MQ=255)                                                   ggACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCcgat                                                                                                                                 <  1:116338/85‑1 (MQ=255)                                                   ggACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCcgat                                                                                                                                 >  2:116338/1‑85 (MQ=255)                                                   ggACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGc                                                                               >  7:12188/1‑135 (MQ=255)                                                   ggACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGc                                                                               >  3:22144/1‑135 (MQ=255)                                                   ggACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGc                                                                               >  3:92576/1‑135 (MQ=255)                                                   ggACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGc                                                                               >  2:175402/1‑135 (MQ=255)                                                           tGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGc                                                                       <  7:371141/135‑1 (MQ=255)                                                            gACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGc                                                                                                                     >  3:100895/1‑88 (MQ=255)                                                            gACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGc                                                                                                                     <  4:100895/88‑1 (MQ=255)                                                            gACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCt                                                                      >  1:232358/1‑135 (MQ=255)                                                            gACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGAAGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCt                                                                      >  8:147140/1‑135 (MQ=255)                                                            gACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGATCCCCACTATTTCCGATGCGATGAAAAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCt                                                                      >  8:1210/1‑135 (MQ=255)                                                                tGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGaa                                                                  <  4:169159/135‑1 (MQ=255)                                                                 gAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAat                                                                 >  8:66303/1‑135 (MQ=255)                                                                      gtcgtcCCCCCTTCAACGCCAGACTCGCCCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCCACCCTAACCCAGCGGCTGCGCCCTGTTCTTCTGCCGCGTAGCGTGAATATCTg                                                            >  1:168604/1‑135 (MQ=255)                                                                       tcgtcCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGc                                                           <  1:175402/135‑1 (MQ=255)                                                                       tcgtcCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGc                                                           <  7:257730/135‑1 (MQ=255)                                                                          tcCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGcc                                                                              <  3:234863/113‑1 (MQ=255)                                                                          tcCCCCCTGCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGcc                                                                              >  4:234863/1‑113 (MQ=255)                                                                            cccccTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCCACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCctct                                                      >  1:73230/1‑135 (MQ=255)                                                                            cccccTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCctct                                                      >  8:308718/1‑135 (MQ=255)                                                                            cccccTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCctct                                                      >  4:180591/1‑135 (MQ=255)                                                                            cccccTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGACCCAACCCTAACCCAGCGGCTGAGCCCTTATCTTCTGCCGCGTAGCTTCAAAATCTGCCctca                                                      >  3:331445/1‑134 (MQ=255)                                                                               ccTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGc                                                   <  7:20443/135‑1 (MQ=255)                                                                               ccTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGc                                                   <  7:308718/135‑1 (MQ=255)                                                                                 ttCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGAAGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGCtt                                                 <  7:147140/135‑1 (MQ=255)                                                                                     aCGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGCTTGAAc                                             <  5:90012/135‑1 (MQ=255)                                                                                     aCGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGCTTGAAc                                             <  7:66303/135‑1 (MQ=255)                                                                                           gACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGCTTGAACAGATAg                                       <  2:134758/135‑1 (MQ=255)                                                                                            aCTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGCGCCCTTTTCTTCTGCCGCGTAGCTTGAATATATGCCCGCTGGCGTTAACAGATAgg                                      >  5:292584/1‑135 (MQ=255)                                                                                                 aCCTCTATCCTGCCACCATGGCTGTCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGCTTGAACAGATAGGGCATc                                 <  6:182255/135‑1 (MQ=255)                                                                                                 aCCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGCTTGAACAGATAGGGCATc                                 <  8:237874/135‑1 (MQ=255)                                                                                                    tctATCCTGCCACCATGGCTGGCGACTATTTCCGAAGCGATGAATAGc                                                                                                                     >  7:41451/1‑48 (MQ=255)                                                                                                    tctATCCTGCCACCATGGCTGGCGACTATTTCCGAAGCGATGAATAGc                                                                                                                     <  8:41451/48‑1 (MQ=255)                                                                                                       aTCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCTACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGCTTGAACAGATAGGGCATCACAt                             <  4:96666/133‑1 (MQ=255)                                                                                                       aTCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGCTTGAACAGATAGGGCATCACAt                             <  3:372834/133‑1 (MQ=255)                                                                                                          cTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGCTTGAACAGATAGGGCATCACATCGCTg                        >  4:180759/1‑135 (MQ=255)                                                                                                               accaTGGCTGGCGACTATTTCCGATGCGATGAATAGc                                                                                                                     >  2:301877/1‑37 (MQ=255)                                                                                                               accaTGGCTGGCGACTATTTCCGATGCGATGAATAGc                                                                                                                     <  1:301877/37‑1 (MQ=255)                                                                                                                     gctggcGACTATTTCCGATGCGATGAATAGCCCCAACCCTAAcc                                                                                                        <  7:305252/44‑1 (MQ=255)                                                                                                                     gctggcGACTATTTCCGATGCGATGAATAGCCCCAACCCTAAcc                                                                                                        >  8:305252/1‑44 (MQ=255)                                                                                                                     gctggcGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGc                                                                       <  1:94107/77‑1 (MQ=255)                                                                                                                     gctggcGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGc                                                                       >  2:94107/1‑77 (MQ=255)                                                                                                                                 ttCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGCTTGAAAAGATAGCGCATCACATCGCTGGGAATCGGCTCTCCGGAATtgt  >  7:302117/1‑134 (MQ=255)                                                                                                                                       |                                                                                                                                GGTTTCCTAGCTGTCAACGGTAGGCGCTTCCAGTCGGACTGAAGCCATCAGGACAGAATGACCTGAAAGGTCGTCCCCCCTTCAACGCCAGACTCGACCTCTATCCTGCCACCATGGCTGGCGACTATTTCCGATGCGATGAATAGCCCCAACCCTAACCCAGCGGCTGAGCCCTTTTCTTCTGCCGCGTAGCTTGAATATCTGCCCTCTGGCTTGAACAGATAGGGCATCACATCGCTGGGAATCGGCTCTCCGGAATTGT  >  ACB122/4468440‑4468701

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑