Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	1,626,923	C→T	6.9%	S272S (AGC→AGT)	lepB →	Signal peptidase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	1,626,923	0	C	T	6.9%	187.2 / 3.5	58	S272S (AGC→AGT)	lepB	Signal peptidase I
Reads supporting (aligned to +/- strand):  ref base C (22/32);  new base T (2/2);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

CAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATC  >  KT2440/1626794‑1627053                                                                                                                                  |                                                                                                                                   ctaCTCCAACTTTAGCCGTTGCTGGTAGGCCCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTg                                                                                                                              <  7:387902/134‑1 (MQ=255)     tCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGccc                                                                                                                           <  4:322535/135‑1 (MQ=255)       cAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCaa                                                                                                                         >  4:391227/1‑135 (MQ=255)          cGATCGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCt                                                                                                                      <  5:345325/135‑1 (MQ=255)          cGATAGCCGTTGCTGGGATGACCCGACCATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCt                                                                                                                      <  5:425415/135‑1 (MQ=255)                 cGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCgg                                                                                                               >  5:167085/1‑135 (MQ=255)                  gTTTCTGGGATGACCCGAACATTCACAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGTCAAGGCCTTCGCGGTGTGGATGATCTGGACAGAGCCCAAACTCAGCCACTTGCCCAAACTGTCGCggg                                                                                                              >  2:386693/1‑135 (MQ=255)                     tCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTg                                                                                                           <  8:10882/135‑1 (MQ=255)                         ggATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCt                                                                                                       <  1:231239/135‑1 (MQ=255)                          gATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGcc                                                                                                        <  5:194853/133‑1 (MQ=255)                          gATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCt                                                                                                       <  2:136075/134‑1 (MQ=255)                          gATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCt                                                                                                       <  4:178566/134‑1 (MQ=255)                               cccGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCa                                                                                                 <  2:369625/135‑1 (MQ=255)                                ccGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCAt                                                                                                <  4:47377/135‑1 (MQ=255)                                 cGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTg                                                                                              >  7:86059/1‑136 (MQ=255)                                 cGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTg                                                                                              >  8:66912/1‑136 (MQ=255)                                  gAACATTCCCAATGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTg                                                                                              <  3:389566/135‑1 (MQ=255)                                     cATTCCCACGGAACTGCTCGGCATGTTTCCGGACCGGATCATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTTATCCATTGAt                                                                                            <  4:291736/134‑1 (MQ=255)                                      aTTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATAcc                                                                                         <  7:79146/136‑1 (MQ=255)                                       ttCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATAcc                                                                                         <  4:131573/135‑1 (MQ=255)                                           cAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGGTCCATTTAAACCCAta                                                                                     >  2:408447/1‑134 (MQ=255)                                           cAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCg                                                                                    <  8:60514/136‑1 (MQ=255)                                           cAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGGTGGCCTGCTCCATTGATACCCATc                                                                                     >  2:427631/1‑135 (MQ=255)                                            taGGAAGTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCg                                                                                    <  1:255632/134‑1 (MQ=255)                                            aaGGAACTGCACGGCTTGGTTCCGGACGGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCg                                                                                    <  3:225482/135‑1 (MQ=255)                                            aaGGAACTGCACGGCATTGTTCCGGACAGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATct                                                                                    <  5:134695/135‑2 (MQ=255)                                            aaGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCg                                                                                    <  5:15355/135‑1 (MQ=255)                                            aaGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCg                                                                                    <  6:370223/135‑1 (MQ=255)                                              ggAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGgc                                                                                  <  2:159628/135‑1 (MQ=255)                                              ggAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGgc                                                                                  <  6:113610/135‑1 (MQ=255)                                                     cACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTGca                                                                          >  8:73145/1‑136 (MQ=255)                                                     cACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGGTGGGCGGCTCCATTGATTCCCATCGgcgcggggc                                                                           >  5:429454/1‑130 (MQ=255)                                                       cGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCGGTcccc                                                                         >  6:211509/1‑133 (MQ=255)                                                             ggTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCATCCTGTCGCGGGTTGGGCGGCTCCATTGCTACCCATCGGCGCTGTGCACCCGACAg                                                                  >  8:135780/1‑136 (MQ=255)                                                                tCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCTCACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATACAATGAAACCCCTCCGAGCGGTCCAGCCCACAgcg                                                                >  1:420169/1‑135 (MQ=255)                                                                       cggaACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGCTCCATTGTTACCCAACGGCCGTGGCCAGCCGCCACCGCCGaaa                                                          >  6:234031/1‑133 (MQ=255)                                                                       cggaACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCGGTGCAGCAGCCAGCGCCgcagg                                                         >  7:372885/1‑131 (MQ=255)                                                                       cggaACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGGTGGCCTGCTCCATTGATACCCATCGGCGCTGTCCCGCACCCAGCGCCgcaag                                                         >  2:273098/1‑131 (MQ=255)                                                                        ggaACATTGGCGGCACGGCCTGCGCGGTGGGGAGGAGCTGCCCAGAGCCCAAACTCAGCCACTGGCCCACCCTGTCGCGGGTGGGCCTGATCCATTGATACCCATCGCCGCTGCCCAGCAGACAGCGCCGAATGCa                                                       <  8:86059/136‑1 (MQ=255)                                                                         taaCATTGTCGGCAAGGCCTTCGCGGTGTGCATGAGCGGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGGCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCa                                                       <  6:167085/134‑1 (MQ=255)                                                                         gaACATTGTCGGCAAGGCCTTCGCGGGGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCa                                                       <  4:197476/135‑1 (MQ=255)                                                                            cATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGTCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTg                                                                                              <  4:60280/93‑1 (MQ=255)                                                                            cATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGTCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTg                                                                                              <  4:66257/93‑1 (MQ=255)                                                                            cATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGTCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTg                                                                                              >  3:66257/1‑93 (MQ=255)                                                                            cATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGTCACTTGCCCAACCTGTCGCGGGGTGGCCTGATCCATTg                                                                                              >  3:60280/1‑93 (MQ=255)                                                                            cATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGACCAGAGCCCAAACTAAGCCAATTCCCCAAACTGTCGCTGGTTGGCCTGCTCCATTGATACCAATCGGCGATGTCCACCAGACACCGACGCAAGCAttaa                                                   >  8:124042/1‑134 (MQ=255)                                                                                  cGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATAGGCGCTGTGCCGCAGCAAGCGCCGAATGCAtttattaa                                               >  2:364218/1‑129 (MQ=255)                                                                                   ggCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCACCCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATa                                            <  7:445671/136‑1 (MQ=255)                                                                                     cAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATCCCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAg                                           <  5:416880/135‑1 (MQ=255)                                                                                     cAAGGCCTTCGCGGTGTGGAGGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGAGCCATTGATACCCATCGGCGCTGGCCAGCAGACAGCGCCGAATGCATTTCTGACATAg                                           <  3:388623/135‑1 (MQ=255)                                                                                           cTTCGCGGGGTGGAGGCGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTgtgt                                    <  7:209817/136‑1 (MQ=255)                                                                                                       gATGAGCTGGCCAGAGCCCAACCTCAGCCACTTGCCCAACCTATCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGCAAGCCGTTCTGCGATAGGCTGGGTTCGCCGGGAtgg                        >  7:116205/1‑134 (MQ=255)                                                                                                       gATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGGCCAGCAGACCCCGCCGCATGCCTTTGTGACATAGGCTGTGTTCTGAGGGCTCg                        >  8:209817/1‑136 (MQ=255)                                                                                                                gCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGGCCAGCAGCCCGCGCCGCCTGCATTTCTGACATAGGCGGGGTTGTCAGGGGTCGGGAGCTTc                >  5:230906/1‑135 (MQ=255)                                                                                                                 ccacaccccAAACCCCGCAAGTGGCCCAACCTGGCGGGGGTGGGCCTGATCCTTTCAAACCCATCGGCGCTGTCCCCCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCg               <  1:427631/129‑1 (MQ=255)                                                                                                                  cAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGa                                                <  5:230546/101‑1 (MQ=255)                                                                                                                  cAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGa                                                >  6:230546/1‑101 (MQ=255)                                                                                                                             cTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCCCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATc  >  7:346813/1‑136 (MQ=255)                                                                                                                                  |                                                                                                                                   CAACTCCAACGATAGCCGTTTCTGGGATGACCCGAACATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATC  >  KT2440/1626794‑1627053

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