Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	6,084,764	A→C	22.5%	N178T (AAC→ACC)	paiB →	Protease synthase and sporulation protein PAI 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	6,084,764	0	A	C	22.5%	102.8 / 3.6	59	N178T (AAC→ACC)	paiB	Protease synthase and sporulation protein PAI 2
Reads supporting (aligned to +/- strand):  ref base A (33/11);  new base C (2/13);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 3.90e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGAACAGGGTCGCAGCGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTG  >  ACB122/6084632‑6084897                                                                                                                                      |                                                                                                                                      cctaccAGGGTCTCAGCGAACATTGTAAAGTCACCGACGCCCCTGCCGACCACTTCGAGGGCTTGCTCGGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCATGGTG‑CGCGCAAGCTCAGCCAGAAcc                                                                                                                                     <  1:109173/130‑1 (MQ=255) aCGAACAGGGTCGCAGCGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAAcc                                                                                                                                     <  7:17343/135‑1 (MQ=255)    aaCAGGGTCGCAGCGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCt                                                                                                                                  >  7:178513/1‑135 (MQ=255)      cAGGGTCGCAGCGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTcc                                                                                                                                <  1:87871/135‑1 (MQ=255)       aGGGTCGCAGCTAACCATGTAAAGTCCCCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCg                                                                                                                               <  7:219330/135‑1 (MQ=255)            cGCAGCGAACCATGGAAAGTCACCGACGCCCCTGCCGCCTAAATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCACGGTG‑CGCGCCAGCTCAGCCAGAACCGCACCggg                                                                                                                             >  7:57548/1‑132 (MQ=255)            cGCAGCGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCCAGGTG‑AGCGCAAGCTCATCCAGCACCGCGCCggg                                                                                                                             >  5:26993/1‑132 (MQ=255)                gCGAACCATGGAAAGGCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACa                                                                                                                      <  8:59680/135‑1 (MQ=255)                 cGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCCAGCTCAGCCAGCACCGCTCCGGGCAAGACAt                                                                                                                     >  6:22958/1‑135 (MQ=255)                 cGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACAt                                                                                                                     >  1:157331/1‑135 (MQ=255)                 cGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCCTCTGCCCATCGCGCGCCAGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACAt                                                                                                                     >  3:322159/1‑135 (MQ=255)                           aaaGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGtgtg                                                                                                           >  7:27328/1‑135 (MQ=255)                           aaaGTCACCGACGCCCCTGCCGACAACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCAGCAAGCTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGtgtg                                                                                                           >  1:306787/1‑135 (MQ=255)                            aaGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTgc                                                                                                          <  4:25821/135‑1 (MQ=255)                                  ccGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAg                                                                                                    >  2:407634/1‑135 (MQ=255)                                  ccGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTCGTGTGCGCGacg                                                                                                    >  2:39370/1‑133 (MQ=255)                                    gccGCCCCGGCCTACTACTGCGCGGGCTTGCGCGGCGCCTTGGTCGGTAGCCGTCTGCCCTTGGCGCGCCGGCAAGGTG‑CGGGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCGGGTGGGCGCGAAGGc                                                                                                  <  4:280623/133‑1 (MQ=255)                                     aCGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGCCATCGCGGGTGTGCGCGCAGGcc                                                                                                 >  8:328056/1‑135 (MQ=255)                                         cacagCCTACTACATGGCGGGACTGCCCCGCGCCTGCGGCGGTCTCGGTCTGCCCCCGGCCCGACTGCATGGCG‑CGCGCAAGCTCAGCCAGACCCCCGCGGGGCAAGACAGCGCTGGTGTGCGCGACGCCCTgcc                                                                                             <  6:368527/131‑3 (MQ=255)                                            tGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAgcgc                                                                                          >  8:213010/1‑135 (MQ=255)                                                 aCTACATCGCGGGCAGCCGCCGCGCCATCGTCGGTAGCCGTCGCCCCCGCCCGCGCCGGCACGGGG‑CGCGCAAGCCCAGCCAGACCCGCGCCGGGAAAGACATCGGTGGGGTGCGCGAAGGCCTGCGCGCCAGcc                                                                                     <  8:203011/135‑1 (MQ=255)                                                           gggCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACt                                                                           >  6:347175/1‑135 (MQ=255)                                                           gggCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACt                                                                           >  2:217861/1‑135 (MQ=255)                                                             gCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTgg                                                                         >  7:156808/1‑135 (MQ=255)                                                               aTGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGACCCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGa                                                                        >  5:90780/1‑134 (MQ=255)                                                               aTGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGa                                                                        >  8:434641/1‑134 (MQ=255)                                                                   tCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAAt                                                                    >  4:87859/1‑134 (MQ=255)                                                                          cATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGACGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTcgc                                                            >  1:42223/1‑135 (MQ=255)                                                                          cATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTcgc                                                            >  4:431498/1‑135 (MQ=255)                                                                          cATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGAGCAAGCTCAGCCAGAACCGCACCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCCGCCCCGCTCAACAGGACAAACAACTcgc                                                            >  2:94254/1‑135 (MQ=255)                                                                          cATCGTCGGTATCCGTCTGCCCATCGCGCGCCAGCAAGGTG‑CGCGAAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTcgc                                                            >  1:141093/1‑135 (MQ=255)                                                                               tcgGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCACCTGGACAATCAACTCGCGgcgc                                                       >  7:144946/1‑135 (MQ=255)                                                                                cgGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCACCTCGCGGCGca                                                      >  1:455761/1‑135 (MQ=255)                                                                                            gCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGATCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAAct                                          >  6:321203/1‑135 (MQ=255)                                                                                                aTCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAActct                                        >  4:228010/1‑133 (MQ=255)                                                                                                  cgcgcgCC‑GCAACGGGCCGCCCAAGCCCCGCCAGACCGGCGCCGGGCAAGCCATGGCTGGGGGGCGCGACGGCCGCCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAg                                    <  8:395924/135‑1 (MQ=255)                                                                                                  agaccgcCTGCAAGGGG‑CGCGCACGCTCACCCAGACCCGCCCCGGGCACCCCCTCGCGGGGGGGCGCGACGGCCGGCCCCCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAg                                    <  6:335659/131‑1 (MQ=255)                                                                                                    cgcgcCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGCTCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGa                                  >  6:3296/1‑135 (MQ=255)                                                                                                      cgcCTGGACGGGG‑CGCGCAAGCGCAGCCAGAACCGCGCGGGGCAAGTCATCGCGGGGGGGCGCGACGGCCGGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAAt                                <  6:104799/135‑1 (MQ=255)                                                                                                       gctggCATGGTG‑GGGGCAAGAGCAGCCAGACCCGCTCCGGCCAATACATGGGTGGGGGGGGCGAAGGCCGGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATa                               <  2:373658/131‑1 (MQ=255)                                                                                                           tcAAGGGG‑CGCGAAATCTCAGCCAGAACCGCTCGGGGCAATACATCGGTGGGTGGCGCGACGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGccc                           <  8:224323/134‑1 (MQ=255)                                                                                                             tggggg‑cgcccAAGCCCCCCCACACCCGCTCCGGGCAAGACCGCGGGGGGGGGCGCGACGGCCGGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                         <  6:66047/130‑1 (MQ=255)                                                                                                             ccgggg‑cgcgcAAGCCCCGCCCGACCCGCGCCGGGCAAGACATGGGGGGGGGGCGGGTGGGCCCGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                         <  2:337096/130‑1 (MQ=255)                                                                                                             atggcg‑cgcgcAACCGCACCCAGACCCGCTCGGGGCAAGACTTGGCTGGTGTGCGCGACGGCCTGCCCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                         <  3:53362/130‑1 (MQ=255)                                                                                                             acgggg‑cgcgcAAGCTCAGCCAGACCCGCGGCGGGAAAGCCATCGCTGGGGTGCGGGCCGGCCTGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                         <  2:357798/130‑1 (MQ=255)                                                                                                             acggcg‑cgcgcAAGCTCAGCCAGCACCGCTCGGGGCAAGACTTCGCTGGGGGGCGCGAAGGCCCGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                         <  7:111939/130‑1 (MQ=255)                                                                                                             aaGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                         >  5:222459/1‑135 (MQ=255)                                                                                                              ccgggccgcgcCAGCCCACCCCGCCCCGCCCCGGGCAAGACCTCGCTGGCGCGCGCGCCGCCCGGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                         <  7:213010/129‑1 (MQ=255)                                                                                                              caggcccgcccAAGCTCACCCAGACCCGCTCCGGGCACGTCAGCGCGGGCGCGCGCCCCCCCCGCACCCCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                         <  6:222459/129‑1 (MQ=255)                                                                                                                   gctcgcAAGCTCAGCCAGAACCGCTCCGGCCAAGACATCGCGGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCa                    <  8:339822/132‑1 (MQ=255)                                                                                                                       gcAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCTCATGCGCCAACTCTGAAGGAATAGCCCATGCCCCAcc                 >  3:272654/1‑134 (MQ=255)                                                                                                                       gcAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGccctgtct                >  3:260562/1‑130 (MQ=255)                                                                                                                       gcAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGccctgtct                >  5:335659/1‑130 (MQ=255)                                                                                                                        cAAGCACAGCCTCACCCGCTCGGGGAAACACATCGCTGGGGGGGCCTAAGCCCGGACCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCtt               <  6:90780/135‑1 (MQ=255)                                                                                                                         aaGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTa              >  3:315025/1‑135 (MQ=255)                                                                                                                          aGCTCAGCCAGACCCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTAc             >  7:308664/1‑135 (MQ=255)                                                                                                                               cgCCAGACCCGCGCCGGGCACGCCCTCGCGGGGGGGCGCGCCGGCCGCCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTgc        <  7:126911/134‑1 (MQ=255)                                                                                                                               aGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTgc        >  8:266547/1‑135 (MQ=255)                                                                                                                               aGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTgc        >  2:295724/1‑135 (MQ=255)                                                                                                                                     aaCCGCTCCGGGCAAGACATCGGTGGTGGGCGCGAAGGCGGGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTg  <  6:129691/135‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      ACGAACAGGGTCGCAGCGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTG‑CGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTG  >  ACB122/6084632‑6084897

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